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13  ABSTRACT 

■Forces  and  motions  imposed  upon  a  moored  buoy  system  by  the  oceanic  environment  are 
investigated  using  numerical  methods  on  a  digital  computer.  Equations  of  motion  for  axisym- 
metric  buoys  are  developed  for  six  degrees  of  freedom  including  cross-coupled  hydrostatic 
and  hydrodynamic  forces.  Equations  of  motion  for  the  mooring  cable  are  developed  to  allow 
simultaneous  axial  and  transverse  wave  propagation.  A  unique  method  of  characteristics 
numerical  technique  is  shown  for  digital  computer  simulation  of  cable  dynamics.  The  equa¬ 
tions  of  motion  for  the  buoy  and  its  mooring  cable  are  coupled  and  are  excited  by  winds, 
currents,  and  a  quasi-random  wind  wave  model.  Simulated  buoy  system  parameter  response 
characteristics  are  compared  with  response  characteristics  observed  from  actual  oceano¬ 
graphic  buoy  system  moorings.  JJThis  comparison  indicates  that  steady  state  buoy  system 
forces  and  configurations  can  be  predicted  within  5  percent  and  that  dynamic  motions  can 
be  predicted  within  50  percent  usin^j  the  digital  computer  simulation. 
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ABSTRACT 


Motions  and  dynamic  forces  imposed  upon  a  moored  buoy  system  by  the 
oceanic  environment  are  of  vital  interest  to  the  user  of  the  system.  If  instru¬ 
mentation  for  monitoring  the  environment  is  motion  sensitive,  it  is  of  little 
value  if  its  response  to  platform  motion  is  greater  than  its  response  to  the 
changing  environment.  The  buoy  system  designer  is  also  concerned  with 
motions  and  forces  in  buoy  systems  in  order  to  design  for  the  highest  probabil¬ 
ity  of  system  survival  under  extreme  conditions. 

In  response  to  these  needs,  this  study  investigates  digital  computer 

s 

simulation  of  buoy  system  dynamics  for  simple  buoy  systems,  i.e. ,  a  surface 
buoy  moored  on  a  single  mooring  line.  The  buoy  system  can  be  excited  by 
winds,  waves,  and  currents.  Winds  can  act  from  any  compass  direction,  and 
currents  can  vary  in  strength  and  direction  as  a  function  of  depth  in  the  water 
column.  Wind  waves  are  simulated  by  first  computing  their  properties  with  the 
Sverdrup-Munk-Bretschneider  method  and  then  by  using  Borgman's  energy 
partitioning  scheme  on  a  two-parameter  Bretschneider  spectrum  to  compute 
component  sine  wave  amplitudes  and  frequencies.  Since  the  component 
Stokesian  waves  are  linear,  the  principle  of  superposition  can  be  used  to  sum 
component  magnitudes  in  order  to  cc  mnute  v  Uor  particle  motions. 

Equations  of  motion  for  the  buoy  are  developed  for  six  degrees  of  freedom- 
three  translational  and  three  rotational.  Hydrostatic  and  hydrodynamic  forces 
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and  moments  acting  on  an  oblate  spheroid  moving  on  the  free  surface  of  an 
infinite  body  of  water  are  investigated  in  detail.  The  set  of  integro-differential 
equations  for  buoy  motions  are  reduced  to  a  set  of  nonlinear,  ordinary  differ¬ 
ential  equations  with  nonconstant  coefficients  by  using  the  Haskind  hypothesis 
to  evaluat^  the  hydrodynamic  force  and  moment  integrals  and  to  represent 
them  as  frequency  dependent  coefficients.  Buoy  motions  are  coupled  through 
hydrostatic,  hydrodynamic,  and  mooring  line  forces. 

Cable  dynamics  are  also  investigated.  A  set  of  coupled,  hyperbolic, 
partial  differential  equations  for  cable  motions  are  developed,  and  character¬ 
istic  equations  are  derived  to  effect  a  method  of  characteristics  solution. 

A  unique  numerical  method  of  characteristics  technique,  based  upon  Hartree's 
method,  is  developed  for  the  solution' of  the  cable  equations  in  the  time-space 
domain.  Buoy  motions,  which  are  dependent  upon  the  cable  tensions,  serve  as 
the  upper  boundary  conditions.  'Lower  boundary  conditions  are  prescribed  at 

the  anchor,  where  there  can  be  no  motion. 

/ 

i 

For  certain  buoy  systems,  where  many  mass  discontinuities  exist  along 

/ 

the  cable,  or  for  shallow  water  moorings,  where  slack  cable  conditions  can 
exist,  a  lumped-mass  method  of  computing  cable  dynamics  is  developed  as 
opposed  to  the  finite-difference  method  just  described.  In  general,  for  cable 
dynamics  the  lumped-mass  numerical  method  is  an  order  of  magnitude  faster  in 
computation  time  than  the  finite  difference  method. 

The  equations  of  motion  developed  for  the  buoy  were  solved  numerically 
in  the  time  domain  using  a  fourth-order,  Runge-Kutta  integration  method.  Cable 
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equations  can  be  solved  either  by  finite-difference  methods  or  by  integrating 

■<4> 

with  the  Eunge-Kutta  algorithm  fcr  the  lumped-mass  model. 

In  order  to  validate  the  numerical  models  developed,  two  buoy  systems 
were  instrumented  and  deployed  in  Block  Island  Sound.  The  motion  data  from 
these  experiments,  along  with  data  published  in  the  literature,  are  compared 
with  simulated  buoy  motion  data.  This  comparison  indicates  that  steady-state 
buoy  system  forces  and  configurations  can  be  predicted  within  approximately 
5  percent  and  that  buoy  system  dynamics  can  be  predicted  within  approximately 
50  percent.  There  are  some  indications  that  the  surge  and  sway  hydrodynamic 
forces  acting  on  the  buoy  are  being  underestimated  by  the  computer  model. 
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I.  INTRODUCTION 


Clean,  green,  windy  billows  notching  out  the  sky, 

Grey  clouds  tattered  into  rags,  sea-winds  blowing  high, 
And  the  ships  under  topsails,  beating,  thrashing  by, 
And  the  mewing  of  the  herring  gulls. 

Dancing,  flashing  green  seas  shaking  white  locks, 
Boiling  in  blind  eddies  over  hidden  rocks, 

And  the  wind  in  the  rigging,  the  creaking  of  the  blocks, 
And  the  straining  of  the  timber  hulls. 

John  Masefield 
"Cardigan  Bay" 


This  study  is  concerned  with  the  analysis  and  simulation  of  the  dynamics 
of  simple  oceanic  buoy  systems.  The  analysis  must  include  the  effect  of  the 
significant  forces  that  act  on  the  buoy  system  and  are  imposed  by  the  ocean 
environment  —  wind,  waves,  and  currents.  Because  of  the  highly  nonlinear 
natui'e  of  the  problem,  numerical  methods  are  favored  in  order  to  provide  a 
realistic  simulation. 

Buoys  have  been  employed  by  mariners  for  centuries  as  aids  to  naviga¬ 
tion  and  to  support  mooring  chains.  In  this  country,  navigational  buoys  were 
in  Service  in  the  Delaware  River  in  1767  and  in  Boston  Harbor  by  1808.  At 
present,  the  United  States  Coast  Guard  maintains  over  24,000  buoys  in  the 
navigable  waters  of  the  United  States  and  its  possessions.  Navigational  buoy 
system  design  is  largely  a  matter  of  employing  "rules  of  thumb"  evolved  over 
decades  of  experience  with  these  buoys.  The  vast  majority  of  navigational 
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buoys  are  moored  in  waters  less  than  100  ft  deep,  and  the  most  common  cause 
of  failure  is  collision  with  a  vessel.  Tims,  these  buoy  systems  are  charac¬ 
terized  by  massive  steel  buoys  moored  with  heavy  chain  to  large  concrete 
clumps.  They  are  serviced  on  a  regular  basis  and  are  recovered  and  over¬ 
hauled  annually. 

In  recent  years,  oceanographers  have  used  buoy  systems  to  support 
current  meters,  thermistor  chains, and  other  oceanographic  instrumentation. 

The  state  of  the  art  in  oceanography  has  advanced  to  the  point  where  oceanog¬ 
raphers  are  no  longer  satisfied  with  data  taken  at  a  single  point  over  a  rather 
short  time  duration.  Multiple  measurements  to  be  made  simultaneously  over 
wide  areas  of  the  ocean  or  long-duration  measurements  are  made  most  econom¬ 
ically  with  a  buoy  system  equipped  with  self-recording  or  telemetering  instru¬ 
mentation.  This  economy  can  be  realized  only  if  the  buoy  system  is  designed 
to  have  a  life  on  station  greater  than  the  desired  measurement  time. 

The  basic  design  philosophy  of  deep-sea  oceanographic  buoy  systems  is 
quite  different  from  that  of  the  navigational  buoys;-  Oceanographic  vessels  are 
usually  small  and  are  not  equipped  for  handling  heavy  objects  over  the  side  at 
sea;  thus,  the  buoys  and  mooring  line  components  must  be  of  relatively  light 
weight.  The  mooring  lines  are  miles  long  and  thus  preclude  the  use  of  heavy 
chains  (except  at  the  bottom)  and  tend  to  be  madt-.  up  of  light,  high-strength 
wire  ropes  or  of  synthetic  fiber  ropes.  Designers  of  oceanographic  buoy 
systems  are  faced  with  the  near-impossible  task  of  designing  a  lightweight, 
highly  compliant  structure  to  survive  for  periods  of  a  year  or  more  in  one  of 


the  harshest  environments  known  to  man. 
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Isaacs'1'  describes  a  mean  time  between  failure  (MTBF)  of  121  days  for 

the  taut-moored  Scripps  Institution  of  Oceanography  ’'Catamaran”  buoys.  The 

observed  system  failures  were  due  to  parting  of  the  nylon  mooring  line  near 

the  surface.  The  writers  hypothesize  that  high  tensile  loads  in  the  nylon  line 

2 

are  caused  by  tensile  waves  propagating  up  and  down  the  cable.  Richardson  of 
theWoods  Hole  Oceanographic  Institution  (WHOI),  anchored  106  buoy  systems 
between  Cape  Cod  and  Bermuda.  The  MTBF  for  these  systems  was  found  to 
be  about  90  days.  The  WHOI  buoy  system  failures  were  attributed  to  mooring 
line  failures,  fish  bite  of  synthetic  mooring  line$  and  theft.  In  1967,  WHOf* 
set  nine  long  term  buoy  moorings  of  which  only  one  was  recovered  on  station 
after  60  days.  Three  of  these  were  found  adrift.  WHOI  was  more  successful 

4 

in  1968,  when  only  3  of  14  long-term  buoy  moorings  failed. 

All  the  oceanographic  buoy  systems  described  above  were  taut-moored 
systems,  the  majority  of  which  utilized  synthetic  rope  in  their  mooring  lines. 
Their  short  life  on  station  and  their  low  recovery  rate  indicate  a  need  for  an 
accurate  engineering  method  of  computing  the  dynamic  response  of  the  buoy 
system  to  the  ocean  environment. 

Slack  moored  buoy  systems  have  a  much  greater  MTBF  and  are  more 
reliable.  The  catenary  of  the  mooring  line  provides  the  necessary  compliance; 
thus,  dynamic  tensions  in  the  mooring  line  are  reduced.  Navigational  buoy 

5 

systems  are  slack  moored  and  are  very  reliable.  Smith  cites  a  long  history 
of  successful  moorings  in  the  Gulf  of  Mexico  for  the  NOMAD  buoy  system. 

Over  a  5-year  period,  a  number  of  the  NOMAD  buoys  were  kept  on  station  for 
periods  of  a  year  or  more.  Smith  also  describes  fifteen,  25-ton  barges  that 


were  slack  moored  for  8  months  in  water  depths  over  12, 000  ft  with  no 
failures. 

Oceanographers  are  also  concerned  with  the  effects  of  buoy  system 

motions  on  their  instrumentation.  For  example,  instrumentation  fastened  to 

the  mooring  line  of  a  slack-moored  system  will  undergo  depth  excursions  that 

6  7 

are  dependent  on  the  current  structure.  Webster  ’  discusses  errors  in  self- 

recording  current  meter  data  due  to  buoy  system  motion.  Webster  shows  a 

current  energy  spectrum  developed  from  a  current  meter  that  is  attached  to 

the  mooring  line  of  a  buoy  system  (figure  1).  The  energy  introduced  into  the 

current  data  by  the  buoy  system  motions  is  far  greater  than  the  energy  of  the 

currents  themselves. 

8 

Marcus  compared  anemometer  data  taken  from  a  NOMAD  buoy  in  the 

Gulf  of  Mexico  with  other  meteorological  observations  in  the  area  over  a 

6-month  period.  The  mean  error  of  the  wind  speed  data  was  0. 2  knot  with  a 

g 

standard  deviation  of  the  error  of  +4. 67  knots.  Huff  shows  a  power  spectrum 
of  anemometer  data  taken  from  another  NOMAD  buoy  system  moored  off 
Bermuda.  This  spectrum  (figure  2)  indicates  that  a  large  amount  of  energy 
was  introduced  into  the  spectrum  by  the  motion  of  the  buoy.  Huff  also  shows 
an  increase  in  the  average  wind  speed  deviation  from  the  mean  with  increasing 
mean  wind  speed  that  levels  out  at  high  wind  speeds.  This  variation  is  char¬ 
acteristic  of  sea  surface  slopest  which  have  an  upper  limit  due  to  gravity.  This 
upper  limit  implies  that  the  wind  speed  error  is  due  to  the  pitch  and  roll  mo¬ 
tions  of  the  buoy.  Day*^  found  that  wind  data  sampled  every  10  min  from  a 
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Figure  2.  An  Anemometer  Power  Spectrum 
(From  Webster,  reference  9.) 
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buoy-mounted  anemometer  had  to  be  time  averaged  over  a  2- hr  period  to 
remove  errors  due  to  buoy  motions. 

The  two-dimensional,  steady-state  configurations  of  buoy  cable  systems 

11  12 

have  been  investigated  by  Wilson.  ’  Wilson  constrained  the  upper  end  of 

the  cable  to  be  at  the  mean  ocean  surface  and  did  not  consider  cable  elasticity. 
13 

Patton,  as  part  of  this  dissertation  research,  developed  a  numerical  method 

to  determine  the  equilibrium  configuration  of  buoy  cable  systems.  The  three- 

dimensional  configuration  for  any  current  structure  (currents  may  vary  in 

strength  and  direction  as  a  function  of  depth),  as  well  as  the  buoy  draft,  is 

computed.  The  elasticity  of  the  cable  is  considered  and  the  stretch  is  also 
14 

computed.  Martin  developed  a  numerical  method  similar  to  Patton’s,  but 

it  is  restricted  to  two  dimensions.  Martin  also  experimentally  investigated 

§ 

the  elastic  properties  of  nylon  rope  and  the  drag  of  buoy  models.  Smith 

presented  a  graphical  method  for  two-dimensional  buoy  system  configurations 

but  neglected  tangential  drag  and  elasticity. 

The  study  of  the  motions  of  bodies  floating  on  the  ocean  surface  and  being 

15 

excited  by  waves  originated  in  1749  with  Euler  in  his  classic  work,  Sclentia 
16 

Navalis.  Froude  was  concerned  with  the  rolling  and  roll  stability  of  ships. 

Froude  recognized  the  nonlinearity  of  the  problem  and  included  viscous  resist- 

17 

ance  in  the  equations  of  motion.  Kriloff  investigated  ship  motions  and  wrote 
coupled  equations  of  motion.  Both  Froud'’  ~nd  Kriloff  assumed  that  the  ship 
did  not  influence  the  waves,  which  allowed  them  to  treat  the  hydrodynamic 
properties  of  the  ship  as  a  body  oscillating  on  a  free  surface.  By  far,  the 


single  greatest  problem  in  ship  dynamics  is  the  description  of  the  hydrodynamic 

18 

forces  acting  on  the  ship.  Lewis  introduced  the  "strip  theory"  in  1929, wh  . .. 

allowed  computation  of  three-dimensional  hydrodynamic  characteristics  from 

19  21 

two-dimensional  theory.  Haskind  ’  assumed  that  the  hydrodynamic  equations 
could  be  linearized  in  such  a  manner  that  velocity  potentials  could  be  super¬ 
imposed.  This  method  allowed  the  use  of  three  velocity  potentials:  (1)  incident 
wave  potential,  i.e. ,  the  velocity  potential  of  the  waves  alone;  (2)  diffracted 
wave  potential,  i.e. ,  the  velocity  potential  of  the  body  fixed  on  a  free  surface 
exposed  to  waves;  and  (3)  forced  heave  potential,  i.e. ,  the  velocity  potential 
of  the  body  oscillating  in  still  water.  Naval  architects  currently  favor  the 
Haskind  hypothesis  as  opposed  to  the  earlier  Froude-Kriloff  hypothesis,  which 

assumes  that  body  dimensions  are  small  compared  with  the  waves. 

22  23 

John  ’  wrote  complete  sets  of  coupled  equations  of  motion  for  floating 

bodies  in  harmonic  waves  and  considered  the  influence  of  the  body  on  the  waves 

(after  Haskind).  John  also  included  an  external  force  term  that  could  be  used 

to  describe  a  mooring  cable.  Heave  and  surge  motions  of  a  sphere  were 

computed  for  various  wave  frequencies. 

24 

St.  Denis  and  Pierson  linearized  the  decoupled  equations  of  motion  for 

a  ship  and  investigated  ship  motions  in  confused  seas  by  summing  the  ship's 

25-27 

responses  to  sine  wave  seas  of  different  frequencies.  Korvin-Kroukovsky 

used  the  "strip"  method  to  compute  the  hydrodynamic  characteristics  of  the 

ship  and  included  cross-coupled  hydrodynamic  forces.  A  complete  discussion 

of  the  state  of  the  art  in  the  prediction  of  ship  motions  is  presented  in  the 

28 

proceedings  of  the  fifth  symposium  on  naval  hydrodynamics.  Current 


8 

research  indicates  that  cross-coupled  hydrodynamic  forces  are  the  same  order 

of  magnitude  as  other  hydrodynamic  forces  and  can  not  be  neglected. 

A  large  oceanographic  buoy  was  built  in  1965  by  General  Dynamics/ 

29 

Convair  Division.  As  part  of  the  design  process,  model  tests  were  conducted 
in  a  towing  tank  for  various  buoy  hull  shapes,  and  buoy  motions  were  simulated 
on  an  analog  computer.  The  analog  computer  simulation  considered  the  dynam¬ 
ics  of  the  planar  motions  of  the  buoy  alone;  the  mooring  line  was  treated  as  an 
elastic  spring.  The  following  quote  is  from  reference  29: 

•  The  simulation  was  not  fully  successful;  some  results  are 
considered  inconclusive.  Due  to  the  difficulty  in  obtaining 
reasonable  agreement  with  the  model  data,  the  analog 
computer  study  was  terminated  short  of  its  goal. 

The  40-ft-diameter  "MONSTER”  buoy  described  in  reference  29  has  proven  to 

be  a  successful  ocean  data  station.  At  sea  motions  of  this  buoy  are  described 
30  31  32 

by  Devereux  and  Uyeda.  Gaul  and  Brown  correlated  buoy  heave  accel¬ 
eration  power  spectra  from  the  "MONSTER"  buoy  and  from  a  small  wave 
sensing  buoy. 

33 

Paquette  developed  a  two-dimensional,  lumped-mass  analog  computer 
model  for  buoy  system  dynamics.  The  buoy  was  assumed  to  follow  an  elliptical 
orbit  (major  axis  vertical  and  equal  to  the  wave  height),  and  its  motions  were 
not  integrated  as  part  of  the  system  dynamics.  The  lumped-mass  cable  elements 
were  acted  upon  by  tensions  on  adjacent  elements,  weight  and  buoyancy  forces, 
and  velocity-squared  drag  forces.  Cable  hydrodynamic  masses  and  linear 


damping  forces  were  neglected,  Paquette  concluded  that  at  least  ten  mass 

elements  are  needed  for  a  deep-sea  mooring  line  to  adequately  describe  the 

system  dynamics  in  the  band  of  ocean  wave  frequencies  that  were  considered 

(0  to  0. 5  Hz).  The  coupling  of  tensile  waves  into  transverse  waves  due  to  the 

steady-state  curvature  of  the  cable  was  also  noted. 

34 

Bivens  and  Swann  also  developed  a  two-dimensional,  lumped-mass  simu¬ 
lation  of  buoy  system  dynamic  but  included  the  buoy  dynamics.  However, 

35 

hydrodynamic  cross-coupled  terms  were  neglected.  Rudnick  measured 

motions  of  the  "FLIP"  spar  buoy  at  sea  and  compared  motion  power  spectra 

with  the  power  spectra  predicted  from  a  linear,  decoupled  buoy  motion  model. 

36 

Blumberg  and  Osborn  developed  a  digital  computer  simulation  for  submerged 

buoy  motions.  The  mooring  line  was  considered  to  be  a  rigid,  massless,  and 

dragless  link.  Hydrodynamic  forces  acting  on  the  buoy  included  no  cross- 

37 

coupled  terms, and  the  equations  of  motion  were  linearized.  Millard  describes 
tension  measurements  made  at  sea  as  part  of  the  Woods  Hole  Oceanographic 
Institution  buoy  reliability  program.  Tension  amplitudes  were  correlated  with 
recorded  currents  and  wind  speeds.  Millard's  data  indicate  that  the  tension 
amplitudes  are  attenuated  with  length  down  the  mooring  cable.  A  very  compre¬ 
hensive  study  of  buoy  system  dynamics  has  been  conducted  by  Prof.  Nath,  of 

38 

Oregon  State  University.  The  two-dimensional  motions  of  a  buoy  and  its 
mooring  cable  were  considered.  Buoy  motions  were  solved  by  use  of  recur¬ 
rence  formulas  and  served  as  boundary  conditions  for  the  cable.  Cable  dynam¬ 
ics  were  solved  by  using  a  numerical  method  of  characteristics  solution. 


Hydrodynamic  forces  acting  on  the  buoy  and  cable  were  included, and  nonlinear 
stress-strain  properties  of  the  cable  were  used.  Transfer  functions  between 
wave  spectra  and  line  tension  spectra  along  the  cable  were  developed  and 

compared  with  MONSTER  buoy  data. 

39 

Hsu  and  Blenkarn  utilized  momentum  flux  equations  to  compute  the 

hydrodynamic  forces  acting  on  a  moored  ship.  Equations  of  motion  were  solved 

numerically,  and  the  mooring  lines  were  considered  as  elastic  springs.  Each 

wave  was  assumed  to  impart  an  impulse;  thus,  the  forcing  function  was  composed 
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of  a  series  of  impulses  acting  on  the  ship.  Burke  assumed  that  a  set  of  linear 
response  functions  for  the  vessel  were  known  and  developed  sets  of  statistical 
relations  for  vessel  motions  in  a  random  sea.  This  technique  was  applied  to 
predict  drilling  vessel  motions,  and  the  results  were  compared  with  drilling 
vessel  motions  recorded  at  sea. 

The  present  investigation  has  produced  a  three-dimensional,  numerical 
model  for  buoy  system  dynamics.  The  model  includes  cross-coupled  hydro- 
dynamic  forces  and  can  be  excited  by  wind,  current,  and  wave  forces.  Cable 
dynamics  are  investigated  with  both  finite- element  (lumped  masses)  and  finite- 
difference  (distributed  mass)  methods.  In  this  study,  finite-element  methods 
were  found  to  be  attractive  because  of  their  relative  economy  with  regard  to 
numerical  computational  time.  Finite-difference  methods,  while  more  rigorous, 
are  more  involved  numerically  and  require  relatively  large  amounts  of  computer 
time. 

The  deterministic  model  is  excited  by  a  numerical  wave  model  having  the 
same  spectral  characteristics  as  the  ocean  waves.  The  computed  buoy  system 
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response  is  then  sampled  to  provide  motion  spectra  by  using  Fast  Fourier 
Transform  (FFT)  techniques. 

To  validate  the  model,  two  oceanographic  buoys  were  equipped  with 
motion  sensing  instrumentation  and  installed  in  Block  Island  Sound.  Buoy 
motions  were  monitored  and  recorded  for  various  wind,  current,  and  wave 
conditions.  These  data  are  compared  with  buoy  motions  computed  with  the 
numerical  model  for  the  same  environmental  conditions. 
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II.  PROCEDURE 


In  order  to  predi^*-  the  response  of  the  buoy  system  (figure  3)  to  the  ocean 
environment,  a  deterministic  model  of  the  system  dynamics  must  be  constructed 
and  excited  by  a  random  model  of  the  oceanic  conditions.  The  two  major  struc¬ 
tural  components  of  the  system  (the  buoy  and  the  mooring  line)  are  treated 
separately  and  then  are  combined  to  form  the  deterministic  buoy  system  model. 

2. 1  System  Dynamics 

Buoy  motions  can  be  described  by  the  equations  of  motion  for  a  body  with 
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six  degrees  of  freedom  floating  on  the  free  surface  of  a  fluid.  The  major 
problem  encountered  in  the  solution  of  the  set  of  six,  coupled,  elliptical,  differ-, 
ential  equations  of  motion  is  the  description  of  the  hydrodynamic  forces  acting 
on  the  buoy.  In  their  most  rigorous  form,  the  buoy  equations  of  motion  would 
be  integro- differential  equations  since  the  dynamic  pressures  must  be  integrated 
over  the  immersed  surface  of  the  buoy.  Analytical  solution  of  these  equations 
of  motion  for  an  arbitrary  body  in  a  random  sea  state  has  not  been  accomplished 
up  to  this  time.  If  the  hydrodynamic  forces  can  be  expressed  as  variable  coef¬ 
ficients  in  the  equations  of  motion,  the  equations  can  be  written  as  a  set  of  six, 
ordinary  differential  equations  that  can  be  solved  by  using  the  approximate 
methods  of  numerical  techniques.  The  problem  now  is  to  define  the  variable 
hydrodynamic  coefficients.  If  the  fluid  is  assumed  to  be  incompressible  and 
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* i 

*  *  irrotational,  the  hydrodynamic  properties  of  certain  simple  two-  and  three- 

dimensional  bodies  can  be  computed.  Also,  if  the  waves  on  the  free  surface 
of  the  fluid  are  deterministic  and  linear,  the  velocity  potentials  for  various 
body  motion  modes  can  be  superimposed  to  construct  the  case  of  a  body  floating 
on  a  free  surface  and  responding  to  waves  propagating  on  that  free  surface. 

This  study  employs  the  technique  described  above.  Equations  of  motion 
are  written  as  a  set  of  six,  coupled,  ordinary  differential  equations  with  variable 
coefficients;  hydrodynamic  coefficients  are  computed  by  assuming  that  the  fluid 
is  incompressible  and  irrotational  and  that  velocity  potentials  can  be  super¬ 
imposed;  aerodynamic  and  hydrodynamic  viscous  forces  are  assumed  to  follow 
a  velocity-squared  law;  and  the  body  is  assumed  to  be  axisymmetric  about  a 
vertical  axis  (as  are  most  oceanographic  and  navigational  buoys),  which  simpli¬ 
fies  the  computation  of  the  hydrodynamic  coefficients. 

Dynamics  of  cables  are  investigated  and  simulated.  The  most  direct 
approach,  i.e. ,  solution  of  the  cable  equations  of  motion  by  a  finite-difference 
method,  is  developed  first.  Since  the  cable  equations  are  a  set  of  nonlinear, 
hyperbolic,  partial  differential  equations,  analytical  solutions  are  intractable, 
and  a  numerical  method  of  solution  is  devised.  Although  more  accurate,  the 
finite-difference  method  can  be  very  expensive  with  regard  to  digital  computer 
time.  A  lumped-mass  simulation  of  cable  dynamics  is  also  investigated  and 
developed.  Lumped-mass  methods  offer  significant  savings  in  computational  time 
at  the  expense  of  truncation  of  the  higher  frequency  cable  dynamics. 

The  buoy  equations  of  motion  and  the  two  sets  of  cable  equations  (finite- 
difference  and  lumped- mass)  are  then  coupled  and  solved  numerically  on  a 
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UNIVAC  1108  digital  computer.  A  numerical  model  of  ocean  waves  is  used  to 
excite  the  buoy  system  dynamics  model.  Steady-state  buoy  system  configura¬ 
tions  are  solved  as  the  zeroth-order  case  of  buoy  system  dynamics. 

2. 2  Experimental  Validation 

In  order  to  validate  the  computer  model,  two  oceanographic  buoys  were 
equipped  with  motion  sensing  instrumentation  and  were  monitored.  The  recorded 
buoy  motion  data  have  been  reduced  in  statistical  form  and  will  be  correlated 
with  buoy  motions  predicted  from  the  computer  models.  Also,  buoy  motion  data 
reported  in  the  literature  have  been  used  to  validate  the  computer  models. 


III.  ANALYTICAL  DEVELOPMENT  AND  DISCUSSION 


3.1.  Buoy  Dynamics 

Consider  an  axisymmetric  buoy  having  six  degrees  of  freedom,  floating 
on  the  free  surface  of  a  fluid,  constrained  by  a  mooring  line,  and  exposed  to 
wind,  waves,  and  currents  (figure  4).  The  buoy  is  being  acted  on  by  the 
following: 

Inertial  forces  and  moments 
Hydrostatic  forces  and  moments 
Hydrodynamic  forces  and  moments 
Wind  forces 
Mooring  line  tensions. 

The  inertial  forces  can  he  separated  into  those  due  to  gravity  (weight) 
and  those  due  to  the  motion  of  the  buoy.  Hydrostatic  forces  can  be  obtained  by 
integration  of  the  hydrostatic  pressure  acting  on  the  submerged  surface  of  the 
buoy.  Likewise,  hydrodynamic  forces  can  be  obtained  by  integration  of  the 
hydrodynamic  pressures  acting  on  the  submerged  surface  of  the  buoy.  Hydro- 
dynamic  forces  are  classified  as  inertial  or  dissipative.  Energy  is  being 
dissipated  through  viscous  effects,  radiation  of  pressure  waves,  and  radiation 
of  surface  waves  generated  by  the  motion  of  the  buoy.  Energy  is  introduced  to 
the  system  through  hydrostatic  and  hydrodynamic  forces  due  to  currents  and 

f 

<  surface  waves. 
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CENTERS 

1.  CENTER  OF  GRAVITY 

2.  WIND  CENTER  OF  PRESSURE 

3.  CENTER  OF  BUOYANCY 

4.  WATER  C  NTER  OF  PRESSURE 

5.  MOORING  LINE  CONNECTION  CENTER 

A  FREE  BODY  OF  THE  BUOY 


COORDINATE  SYSTEMS  *° 


Figure  4.  Buoy  Coordinate  Systems 
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Wind  forces  acting  on  the  exposed  surface  of  the  buoy  also  introduce 

> 

energy  into  the  system.  Mooring  line  tensions  indicate  a  path  of  energy  removal 
from  the  system.  If  the  buoy  were  treated  as  a  "black  box"  that  transforms 
energy  from  one  form  to  another,  we  can  draw  a  schematic  as  shown  in  figure  5. 

From  figure  4,  it  is  seen  that  four  sets  of  coordinates  must  be  considered. 
It  is  desired  to  solve  for  the  coordinates  of  the  center  of  gravity  of  the  buoy  in 
inertial  coordinates,  (x^,  yog,  z^),  but  the  hydrostatic  and  hydrodynamic 
forces  and  moments  are  due  to  fluid  motions  relative  to  the  buoy.  Any  point  in 
space  has  coordinates  x.«  z^  relative  to  (i  =  0,  1,  2,  3);  therefore, 


Xo“Xog  +  xi 

y„ = yog  +  y, 

Z0°Z0g+Zi  ' 


The  only  difference  between  R  and  R  is  a  space  rotation  about  the 
axes  of  the  buoy.  Thus,  we  have 


*. 

y* 

=  n- 

yt 

.2, 

(1) 


where  SL  is  a  3-by~3  orthogonal  rotation  matrix, 

'cos  y  cos  a 


A  - 


-SIN  Y  c  os  A 


-SIN  A 


-  COSrSiNJS  SIN<K 
+  SIN  Y  cosof 


COS  If  SINJB  COScK 
+  SINYSIN<X 


SINYSINA  SWO f  - SIN  YSINAcOSK 

+  COSY  COS  +CO$YSINo< 


-  cos  A  s,A/<* 


COS^S  COS  o< 


(2) 
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ENERGY  TRANSMITTED 
(CABLE  STRAIN  ENERGY) 

*  ENERGY  DISSIPATED 
T  1 .)  PRESSURE  WAVE  GENERATION 

2. )  SURFACE  WAVE  GENERATION 

3. )  VISCOUS  DAMPING 


MATRIX  EQUATIONS  OF  MOTION 
MQ  =»  MG  -B-H+W-T 
WHERE: 

M  -  MASS  MATRIX 
Q  -  ACCELERATION  VECTOR 
G  -  GRAVITY  VECTOR 

B  -  HYDROSTATIC  FORCE  &  MOMENT  VECTOR 
H  -  HYDRODYNAMIC  FORCE  &  MOMENT  VECTOR 
W  -  WIND  FORCE  VECTOR 
T  -CABLE  TENSION  VECTOR 


ENERGY  IN 


WIND 

WAVES 

CURRENTS 


Figure  5.  Energy  Flow  in  Buoy  System 
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<  * 


Since  fl  is  orthogonal,  we  see  that 


;  thus, 


COSTS  COS  £  -SIN  Y  COS J  SIN; 8 


-COSVSINJZSIN# 
+$IN  y  COS  <X 


SINS  SIN  £  SIN  Of  -COSJ&  SINK 

+  COS  y  COS  X 


COS  tSINjCOSX  - SIN  y  SINA  COS *  C<?5^  COS  * 

iSlNTSlNX  +  COO'S  SINK 


(3) 


The  spatial  coordinates  become 


*0  -  XoS  +  x*  cos  y  cos^a  -YtS! w  y  cos^  -  H*  sw>S  ;  (4A) 

>0  5=1  >05  +  X2  (-COS  YSINJ&  SIN  tx  -f-  S»M  ^  COS°c) 

+  YZ(SIN  K$lN/6COSoc  +  C0SV'SlHoc)  +  Z^-COSfc  SlN<x)  )  (4B) 

rW 

-2(7  *  -2(3^  +  Xj.  (cos  y  COS X  +  SIN  y  SIN  ex) 

+  (-SINTfSlNjl  cosxi-  cosy  SiNck)  +  Zz(cosJl  COSx).  (4C) 

The  center  of  the  coordinate  system,  which  is  aligned  with  the  waterplane 
area,  is  located  directly  above  the  center  of  gravity  of  the  buoy.  Motions  of 
fluid  particles  due  to  waves  are  described  relative  to  this  coordinate  system. 

The  dynamical  equations  of  motion  will  be  written  in  buoy  coordinates,  but 
displacements  will  be  transformed  to  the  Rq  coordinate  system  in  order  to 
solve  for  cable  tensions. 

Using  the  free  body  of  the  buoy  (figure  4)  and  applying  Newton's  Second 
Law,  we  can  develop  the  equations  of'  motion  for  the  buoy.  In  matrix  form,  the 
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equations  of  motion  for  the  buoy  are 


MQ  =  MG  -B  -H  -W  -T  ,  » 

where 

M  is  the  structural  mass  matrix 

#• 

Q  is  the  acceleration  vector 
G  is  the  gravitational  acceleration  vector 
J3  is  the  hydrostatic  force  vector 
'}{  is  the  hydrodynamic  force  vector 
w  is  the  wind  force  vector 
T  is  the  mooring  line  tension  vector. 

Each  of  these  forces  will  be  considered  in  turn. 

Using  a  coordinate  system  with  the  origin  located  at  the  center  of  gravity 
of  the  buoy  and  including  moments  and  products  of  inertial,  we  can  write  the 
structural  mass  matrix  as 


M  = 


0 

0 

0 

0 

0 

0 

0 

0 

0 

•U* 

u 

Irt 

(6) 


(10) 


L  ~iJ>n  ClAzmv\ 

where 

D  is  the  drag  force 

L  is  the  lift  force 

J*  is  the  air  density 

C  is  the  drag  coefficient  (subcritical,  5 x IOS  ) 

*  S** 

dL  is  the  lift  coefficient  (subcritical,  ) 

A,  is  the  vertical  projected  area 
A2  is  the  horizontal  projected  area 
w  is  the  wind  velocity 
A*  is  the  absolute  viscosity  of  the  air 
fa  is  the  buoy  width  . 

Because  of  the  axial  symmetry,  the  drag  and  lift  coefficients  of  the  buoy  are  the 
same  in  the  y^  and  directions.  Given  wind  velocity  components  WV  and 
WV  in  the  Rq  coordinate  system  and  assuming  that  the  wind  velocities  are 
an  order  of  magnitude  greater  than  the  displacement  velocities  of  the  buoy,  we 
can  transform  to  the  R  coordinate  system  and  can  compute  the  wind  forces 

Li 

and  moments.  For  forces  acting  on  a  point  that  lies  on  the  axis  of  symmetry  of 
the  buoy,  the  transform  from  inertial  coordinates  to  buoy  coordinates  is  inde¬ 


pendent  of  a  rotations. 
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The  wind  components  in  the  R  coordinate  system  are 


W* 

WM 


.  w„ 
wv; 


Neglecting  the  small  wind  velocity  component  acting  along  the  axis  of  symmetry, 
the  magnitude  of  the  wind  velocity  in  the  R  coordinate  system  is 


and  the  wind  forces  become 


WI  -  iJ>  c.  A,  wv%  WVr, 


WL,  -1^  CL  A,  W*8 1 WV-J  # 


(Note  that  velocities  squared  is  written  as  the  product  of  the  velocity  and  its 
absolute  value  in  order  to  maintain  the  sign  convention. )  Resolving  the  wind 
drag  force  into  y  and  zf  components,  we  find  that 


w  =  .  ur 


.  M 


i 


4  4 
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and 


i»r  *WVia  -r 

w'*  W 


(16) 


The  wind  moments  are  computed  by  usi'.g  the  wind  forces  and  the  moment  arm: 


W^  =  o 


(17A) 


and 


Vfjtt  =  -  V4,  •  Hw 


(17B) 


W*  =  -Wy2  •  Hw 


(17  C) 


where  is  the  height  from  C.  of  Of-  to  the  wind  force  center  of  pressure. 
In  the  R  coordinate  system,  the  wind  forces  and  moments  are 


= 


w 


% 


TV, 


w, 


w. 


(18) 


3. 1. 2  Mooring  Line  Forces 

Mooring  line  tensions  are  acting  on  the  buoy  at  the  mooring  line  termi¬ 
nation  point.  This  point  is  taken  to  be  below  the  center  of  gravity  and  along  the 
axis  of  symmetry  a  distance  from  the  center  of  gravity.  If  the  space 

orientation  of  the  cable  is  described  by  angles  0  and  <f>  relative  to  the  Rq 
coordinate  system  (figure  4),  we  can  develop  a  3-by-3  orthogonal  rotation 
matrix  to  transform  from  inertial  to  cable  coordinates.  The  inverse  of  the 
matrix  can  be  used  to  compute  the  tension  components  acting  on  the  buoy.  This 
rotation  matrix ,  (A),  is  developed  in  the  next  section  on  cable  dynamics  but  is 
used  here. 

The  tension  components  at  the  buoy  end  of  the  cable  are 


_ I 

V 

* 

.-3 

O 

ll 

0 

XT 

0 

In  the  R  coordinate  system,  they  become 

Ct 


si'X  - 


Be 


Tx* 


The  moments  due  to  the  mooring  line  tension  are 


(20) 


ft 
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TA  =  Ti2  •  H* 

^Vi  =  T'yz  '  ^Mt  , 

The  forces  and  moments  due  to  cable  tensions  are  then 


(21B) 


(21C) 


»»  — 

V 

T* 

T* 

0 

Th 

.  \ 


(22) 


3.1.3  Hydrostatic  Forces 

The  hydrostatic  forces  are  considered  next.  If  the  displacements  of  an 
elemental  volume  of  fluid  just  below  the  free  surface  and  in  the  immediate 
vicinity  of  the  buoy  are  given  relative  to  the  inertial  coordinates  (Tj),  and  if 
the  slope  of  the  free  surface  above  this  particle  is  also  given  (£),  the  buoyant 
forces  and  moments  can  be  computed.  The  assumptions  that  the  buoy  diameter 
is  small  relative  to  the  wavelength  and  that  the  presence  of  the  buov  does  not 
influence  the  shape  of  the  free  surface  (Froude-Kriloff  hypothesis)  allow 
representation  of  the  sea  surface  as  a  plane  intersecting  with  the  body  volume. 
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Most  oceanographic  buoys  have  axial  symmetry,  and  their  shape  can 
be  approximated  by  an  oblate  spheroid.  Consider  an  ellipsoid  of  revolution 
(figure  6)  with  a  major  diameter  of  2b  and  a  minor  diameter  of  2a.  The 
equation  of  the  surface  is 


Let  a  plane  intersect  the  oblate  spheroid  at  a  height  from  the  geometric 
center  at  an  angled  .  The  intercepts  of  the  plane  are 

x  =  -  H„ 

y  =  oo 

z.  =  -  Hd  /  tm  A* , 

The  equation  of  the  plane  becomes 


__  4 

-Hd 


Xy  +  Zy  +  Hp  —  0 


(24  A) 


The  intersection  of  the  body  and  the  plane  is  therefore  given  by 


=  Xy  +  +  Hp 


(25) 


i H 


2  DIM'L  SECTION 
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^  +  'i  +  T  -Xi-KitiM-tf-Vl  =0,  «W) 


which  is  an  ellipse.  In  order  to  compute  the  volume  of  the  "cut"  oblate  spheroid, 
the  area  of  any  section  parallel  to  the  cutting  plane  must  be  defined.  In  the 


-  z^  plane,  we  find  that  on  the  ellipse 


X*  =  n* 


and  on  the  line 


5=5  zt  tmzJ3'+  aKp  zb  tnr)J3'+  Hf  # 


/  t  r  2 


Equating  (26)  and  (27),  we  can  compute  the  z  coordinates  of  intersecting  points 
P  and  P  by 

X  £j 

(twn2., 8  +  +  Z .  Z(,HD  tmj}'  +  (Hp-  ftz)-0 


Solving  for  ,  we  find 


-H et*nJB’±  %  /lFfaVF+  Hf 


(tHnlj5'  4  nVt4) 


For  P^  ,  the  coordinates  are 


~Hd  tmJ3'  +  ^4 fl2- H» 
(tnn&jB'-h  nz/hz) 


(31) 


*b,  588  ~  zi>)  ~  H*  j 


for  ,  the  coordinates  are 


H,  = 

i>z 


W2-  H/' 

(tnn^S'r  w^e) 


(32) 


and 


V  =  tnn^S'-  H„  J 


(33) 


thus, 

z  _2  2  *4  Jb*tM*js'+n*-  H? 

^  bt  (tnnsfa+  «Vb*)  #  <34) 

The  b1  axis  of  the  intersecting  ellipse  is 

L'  _  JL  fa,  -  gjj 
2.  cos  Ji' 

From  equation  (34),  we  substitute  and  find 

,/  %  J  b*  t^jf  +  nz-H? 

cos  3 '  ( tnn  *J3 '  +  ft* /y1 )  (35) 


The  z  coordinate  of  the  center  of  the  elliptical  section  is 


-  Hp  tftnM' 
(tmzJ&'+  *z/y*)  * 


2C  “  Z‘t>(  -  k'  cosJS'  =5. 


(36) 
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S.  The  x  coordinate  of  the  center  of  the  elliptical  section  is 

y  —  Hp  trir)zJ3  it 

°  «Vi‘)  "  '  <37> 

The  dimension  of  the  a’  axis  of  the  intersecting  ellipse  is  one- half  the  y 
dimension  at  the  center.  The  y  dimension  at  the  center  is 


v  ~ 

^  /  (tmxJ 84 


/ Hd  tQn7J$' 


(38) 


) 


and  the  axis  length  is 


H7  ~ 


b  - 


h; 


(39) 


The  area  of  the  intersecting  ellipse  is 


a  _  yr  rfb' s  tt  /b* — ^  ^ ~ jUL 


or 


A  -  X."b 

cos 


h; 


-s'  l  (tnn2>S/+ 


(40) 


The  volume  of  the  "cut"  oblate  spheroid  is  given  by 
/>»>' 

V  *  /  A  JH' 


-H 


to 


(41) 


T 


,• : 


I  COSTS'  a  1 '  > 

'-COSfi'JpttfJi'+tf 
or,  when  integrating,  it  becomes 


V=rr*b>[f  +j- 


Ho _ Hp9 


The  location  of  the  centroid  is  given  by 


,/  _  fH'Jv 

CB  7  . 


Expanding,  we  find 


m' _ I  / TThI»  (  1  H'^/cos*#  \u'JU‘ 

ncs  ~  V  /  co*-8'  U*w.V+  "4>)4_  »3j.)* ) 

-  cosp'jg  t  QrfJB'  +  w  3 

Thus,  carrying  out  the  integration,  we  see  that  the  centroid  location  becomes 

cosji  I  H„  2  I 


6  tnnVS'+H*  +  WA - - a* - 

n,>  3ft  WA  n 


34 


The  H  dimension  of  the  centroidal  plane  is 


.  f ,  H> _ u2  .  +  ft*)  7 

ZlZ(bxt*r?Jl>'+tf)  n> _ 2 _ J_ 

[  §  »*'  +  H,  -  -j ^  t  “njy  +  «*)  J 


(45) 


The  x  coordinate  of  the  centroid  is 


■  =  /  tun^'  _  ,  \ 

CS  [(tnn‘ye't  «>t.)  J/' 


The  z  coordinate  of  the  centroid  is 


The  moment  arm  for  the  buoyant  force  is 


d  - 

d,  + 

4 , 

where  the  distances  d 

'l  and  d2 

are  given  by 

4  =7 

/"Tra¬ 

ces 

...  2  * 

♦  cos(  t*'+  JS') 

or 

cos  JS'  -f  XC3  SIN  , 

and 

<4  = 

Hc, 

SinJZ' 

• 

The 

angle  cx 

*  is  defined  as 

/ 

= 

:  tnr)"1 

2-cb 

f 

(46) 


(47) 
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The  moment  arm  is 

<A  *  Z„  COtfi'  (Hca  -  XCJ  SW*',  (4#( 

For  the  special  case  of  a  sphere  of  radius  A  ,  the  immersed  volume  is 

V  =  n  [§■  +  -5s  CM*'  -  COS  **  'J  ,  (49) 

The  height  of  the  plane  intersecting  the  center  of  buoyancy  is 


3  [4-  cos*J'  -  2  Hi  +  h2  sec*J?'] 
C°  V  [2  H-sec*'  +  3Hf-  ^  cos **'J 

The  coordinates  of  the  center  of  buoyancy  are 

xC3  ■  -  cos**'-  Hcs 


•2C»  =  -  Sl«*'  COS*'-  Heg  f  (E2) 

The  buoyancy  moment  arm  is  then 

d  ~  zc8  c'oSJ3'+’  (  Hc$  -  XCfi  )  SIN  JZ  #  (53) 
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If  the  buoy  is  pitched  and  rolled  with  angles  jQ  and  and  the  sea 
surface  slope  at  tbe  buoy  is  A  and  Tf#  ,  the  slope  of  the  sea  surface 
relative  to  buoy  coordinates  is  and  ~  #w  •  The 

angle  between  the  sea  surface  plane  and  the  buoy  vertical  axis  is 
Jft'  —  CO  s'  ( cos  2fs  COS  JBS  With  the  buoy  tilt  angle  defined,  along 

with  the  location  of  the  buoy  center  of  gravity  relative  to  the  sea  suriace,  the 
buoyant  force  and  moment  can  be  computed  for  an  oblate  spheroid.  The  buoyant 
force  B  acts  normal  to  the  plane  of  the  sea  surface  through  the  center  of 
buoyancy.  In  buoy  coordinates,  the  buoyant  force  vector  is 


~B  cos  y5  COS  A 

-B  sin  rs 
~8  cos  Vs  SinJZs 


The  buoyant  moment  becomes 


(54) 


Jvts  “  (B  CDS  Ys  COSA)(zot  «) ~(fW-*c,)(6  »»  ,(55) 

and  the  buoyant  moment  becomes 

JVly  =  (B  Ys  COSjis)(ZCB  SIN£*)~  (Hc©""*e*X  ^  Slli  n  (56) 


The  same  transformations  as  used  above  can  be  applied  to  any  axisymmetric 
buoy  if  the  buoyancy  can  be  defined  as  a  function  of  the  draft  and  tilt  angle  of 
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the  buoy.  The  angular  stability  of  axisymmetric  buoy  hulls,  as  defined  by  the 
locations  of  the  centers  of  buoyancy  and  gravity,  is  discussed  in  appendix  A. 

3.1.4  Hydrodynamic  Forces 

The  hydrodynamic  forces  acting  on  the  buoy  that  are  due  to  the  waves 
incident  on  the  buoy  and  the  motion  of  the  buoy  in  the  fluid  must  be  included  in 
the  equations  of  motion.  Ideally,  the  computation  of  these  forces  should  be  made 
for  a  buoy  moving  in  a  viscous  fluid  exposed  to  a  random  sea  state.  However, 
the  solution  to  this  general  problem  is  not  tractable  and  the  forces  acting  on  a 
buoy  moving  sinusoidally  in  an  ideal  fluid  is  considered  in  this  study.  These 
forces  are  considered  as  being  composed  of  two  components  —  inertial  and 
dissipative.  Energy  is  dissipated  from  the  buoy,  which  is  moving  in  an  ideal 
fluid,  by  the  generation  of  surface  waves  that  radiate  cylindrically  to  infinity. 
Dissipative  forces  due  to  viscosity  will  be  included  in  the  equations  of  motion  as 
separate  force  components. 

The  separation  of  dissipative  forces  into  those  due  to  surface  wave  gener- 

41 

ation  and  those  due  to  viscous  drag  is  supported  by  Havelock,  who  concluded 
through  dimensional  analysis  of  the  decay  of  oscillations  of  a  prism  on  a  free 
surface  that  the  viscous  damping  is  an  order  of  magnitude  less  than  the  damping 

no 

due  to  surface  wave  generation.  Ogilvie  cites  other  model  experiments  which 
support  Havelock's  conclusions.  Ogilvie  also  lists  generalized  equations  of 
motion  for  ships  in  a  seaway  in  which  the  viscous  damping  is  either  neglected 
or  included  as  a  separate  force  term,  which  is  the  current  practice  among 


naval  architects. 
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22  23 

The  analysis  of  the  motion  of  floating  bodies  conducted  bv  John  ’ 

illustrates  the  computation  of  the  hydrodynamic  forces.  Consider  a  mechanical 

system  consisting  of  a  liquid  and  a  partly  immersed  body  B .  The  liquid  is 

assumed  to  be  incompressible  and  to  have'irrotational  motion.  The  free  surface 

extends  to  infinity  in  all  directions  (figure  7).  The  body  B  is  assumed  to  be 

rigid  and  to  describe  a  forced  motion  under  the  influence  of  external  forces. 

The  state  of  the  liquid  is  described  by  the  velocity  potential  4  (x,  y,  z;  t), 

which  satisfies  LaPlace's  equation.  The  boundary  condition  that  the  normal 

velocity  of  the  particles  along  with  the  pressure  is  continuous  across  the  surface 
% 

must  be  satisfied.  In  addition,  the  pressure  on  the  free  surface  is  assumed 
constant  and  equal  to  the  atmospheric  pressure.  Under  these  conditions,  energy 
is  gained  or  lost  by  the  system  only  through  waves  arriving  or  departing  at 
infinity  or  through  the  external  forces. 

The  difficulties  arising  from  the  fact  that  the  velocity  potential,  (j)  ,  is 
a  solution  of  the  potential  equation  determined  by  nonlinear  boundary  conditions 
on  a  variable  boundary  force  linearization  of  the  problem  in  order  to  make  it 
tractable.  Restricting  the  analjJis  to  infinitesimal  motions,  note  that  the  bound¬ 
ary  conditions  become  linear  conditions  for  the  potential  function  <j>  on  fixed 
surfaces  corresponding  to  tiie  rest  or  equilibrium  position.  The  average  free 
surface  lies  in  a  horizontal  plane,  and  the  average  immersed  surface  S°for 
the  body  B  is  for  a  position  of  equilibrium  for  B .  On  the  average  free 
surface,  y  =  0  ,  the  wave  equation  is 


tit  +  5<Py  «0 


(57) 
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<  * 


On  the  average  immersed  surface,  we  find 


a4 

9n 

where 


x,y,z  are  the  coordinates  of  the  center  of  gravity  jf  B 
are  the  angular  displacements  of  B 
Y\  is  the  unit  normal  of  S* 


(58) 


p  are  the  components  of  n 

K 

are  the  components  of  the  moment  of  n  about  the  center  of 
gravity. 


Sue  differential  equatiors  of  motion  for  B  must  also  be  written.  For  small 
perturbations,  they  are  linear,  second-order  differential  equations  in  x  ,  y  ,  z, 

-V  //  m  i 

“  y  & )  &  with  constant  coefficients  and  integrals  of  (p  in  the  inhomogeneous 
part.  They  are  of  the  form 

fit  »  -fff.  1 &  -  S  [rA(r-y) + m 

where 

M  is  the  mass  of  B 

is  the  fluid  density 

rA 

1  is  the  area  of  A  ,  the  intersection  of  the  body  and  free  surface 

r  A 

lx ,  J*  are  the  moments  of  A  about  v  irtical  planes 


Gf  is  the  acceleration  of  gravity. 


«f  - 
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The  integral  expresses  the  hydrodynamic  forces,  whereas  the  right-hand  term 
represents  the  hydrostatic  forces. 

For  a  buoy  undergoing  simple  harmonic  oscillations,  <p  can  be  defined  as 

<p(x,y,Eit)  -  Re[V(*>7>z)  6  j  ^  (60) 

where 

<7*  is  the  angular  frequency  of  the  oscillation. 

For  an  incompressible  fluid,  V*  is  a  solution  to 

vzV~(x>y>z)  =  o  in  y<o  ? 

where  V  is  complex  valued.  For  small  oscillations  of  the  buoy,  the  amplitude 
of  induced  wave  motion  will  be  small  compared  with  the  wavelength.  Thus,  the 
linearized  dynamic  condition  for  <p  on  the  free  surface  is 

^  (x>  2-  S  fc)  +  fyC  @  "  0  5  (611 

where 

%  is  the  f  x  Cv  surface  elevation. 

Equation  (61)  with  the  linearized  kinematic  condition, 


\t  (X,  Bit)  ,  yields 

-  k~V(x>  °>2-)  —  o  on  y  *=  0  }  (62) 


where 


k  is  the  wave  number,  k  =  <ry$  -  ) 

7\  is  the  free  wavelength. 

The  normal  velocity  across  the  immersed  surface  of  the  buoy  is  continuous; 
thus,  we  find 


,  (63) 

which  is  equation  (58)  restated  in  vector  form  where 

Q  is  the  position  vector  of  the  body  C.  of  Gr.  ;  **■  'Y /j 

@  is  the  rotation  vector  of  body  B  ;  0  —  @  ^  @  -f  (?/;/  ^  ^ 

^  is  the  position  vector  of  some  point  on  the  immersed  surface  of  the 
body. 


The  kinematic  condition  is  to  be  satisfied  on  the  immersed  surface  in  the 
undisturbed  position,  i.e., 

6 

ihVfijrf-  -ler[Q'.n+tffa <n)](64: 


Applying  a  Sommerfeld  radiation  condition  at  infinity,  we  find  that  a  disturbance 
in  the  finite  region  should  only  produce  an  outgoing  wave  at  a  lr  ge  distance: 


V(lY,? )  -AW^e^o  }  (65, 

y'  “  £wnw  (z/*)  . 


where 
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To  put  the  equations  in  dimensionless  form,  let 

ft «  nk  the  frequency 

X  -  X/rt 

y/n  space  variables 

£■  =  a/n 

Q  —  Q  /ft  the  space  parameters, 

where  a  is  a  typical  buoy  dimension.  Introduce  the  pressure  function  Uj 
by 

i crVj(x,y,*)/$n  Qj  =  RUjfay,*)  j=l,£,3  (eeAj 

and 

i<rVj(x>y>z)/<rndj  =RUj(x>>',2)  (66B) 

The  boundary  value  problem  is  to  find  a  potential  Uj  (x,  yy 2)^  j  ••••,,  6 

that  is  continuous  in  the  fluid  space  in  such  a  manner  that 

V20.j  (x,y,z)  =  0  in  y  <0  > 

2-  Uj  (x,0,z)-  AUj  (A)O,z)  =  0  owtsiae  5% 

f„  Uo  (*>  y>z)  =  Hj  (X,  y,  z)  on  5°  , 

uj(*i,t>r)-  AjCrU'^e*^-^ o  w  } 

where  h ;  represents  the  prescribed  function  that  depends  on  the  mode  of 
oscillation.  The  hj  are  given  as 

h,(x,y,2)  =  h*  hvU,y,2t)=:  /hj* 

h2(x,y,z)  =  hy  (-s.(x,y,2)  =  zn< -xnE 

h,(x,y,H)  =  Ha  hs(x,y,z)=xny-ynx. 
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The  source  potentials  G  of  unit  strength  in  the  lower  half-space  which  satisfy 
the  sets  of  boundary  conditions  are 

=  Tf1*  R ' -tt « e(y*n)[S/*&)+XM-i2 
"y+% 

where 

co=j(x-%f+(y- \f]*  cu'x^(x-$f+(y+nfj  2 

S0fa  w)  is  the  Struve  function  of  order  0 

)Xofa&)  is  the  Bessel  functions  of  first  and  second  kind  of  order  0  . 

The  solution  to  the  boundary  value  problem  is  now  in  the  form 


(68) 


where  j*  is  the  strength  of  the  distributed  sources  over  the  immersed  surface 
and  is  a  continuous  complex  function.  If  the  forces  and  moments  are  written  as 
components  in  phase  with  the  acceleration  and  velocity,  we  find  that 


5  =  -/«3MQ.-^<rfi5j\fQ  j=  i,2,3  (69) 

w  J  • 


lor  simple  harmonic  motions,  the  body  vectors  become 
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•*  * 
%JT 


QjCt)  =  Ref  QjJrt)e  ]  j=i>2>3 


and 


©.  (t)  -  Re  [  0i  (t)  e  J  j  w  V,5,  £  (70) 


The  moments  are 

Then,  M  and  i'T  are  the  hydrodynamic  mass  and  the  linear  damping  coef¬ 
ficients,  and  I  and  H  are  the  hydrodynamic  mass  moment  of  inertia 
and  rotational  damping  coefficients.  The  dimensionless  coefficients  are 


RrU 


JSf  =j~rg3  "  Im  j=i)z>3  > 

I  =  -jXt-  =  Rejjy,Uj(x*y,a)'(rxn)ds]  j*v,s;6) 


(72) 


(73) 


j ■=lijS>6  (74) 


42 


Kim  has  evaluated  these  integrals  numerically  for  spheroids  of  various 

aspect  ratios  and  has  plotted  the  dimensionless  hydrodynamic  coefficients  versus 

43 

the  frequency  parameter.  Barakat  evaluated  the  Fredholm  integrals  with  an 
approximate  analytical  solution  for  the  case  of  a  sphere  on  the  free  surface. 
Kim's  data  are  shown  in  figures  8,  9,  and  10.  A  curve-fitting  program 


(CUUFIT)  built  into  the  Government  Services  Administration  remote  terminal 


FREQUENCY  PARAMETER  -  A 


FREQUENCY  PARAMETER  -  A 


Hydrodynamic  Mass  and  Wave  Damping  for 
Heaving  Oblate  Spheroids 

(From  Kim,  reference  42.) 


DIMENSIONLESS  PITCH-WAVE  DAMPING 


FREQUENCY  PARAMETER  -  A 
SPHEROID 
H  “  HALF  BEAM 
DRAFT 

A=  ACT2/g 
H^  =  H/  pff  A’4 

0.12  - 


H  “  00 


FREQUENCY  PARAMETER  -  A 


Figure  10.  Hydrodynamic  Mass  Moment  of  Inertia 
For  Pitching  or  Rolling  Oblate  Spheroids 

(From  Kim,  reference  42.) 


computer  system  was  used  to  develop  approximate  equations  for  Kim's  data  for 
a  sphere.  Program  CURFIT  fits  six  curves  to  the  data  by  a  least-squares  fit  of 
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the  candidate  curve's  linear  transform.  The  six  curves  are  of  the  following 
types: 

Y  =  A  -  B-X 

r  -  Atsx 

Y  =  AXB 

Y  =  A  +  B/X 

Y  =  1/(A+B‘X) 

Y  =  X/(A  +  BiX) 

Coefficients  for  each  curve  and  an  index  of  determination  (best  fit)  are  computed. 
Kim's  data  for  the  sphere  are  approximated  by  the  following  functions: 

Sway- Surge 

HYDRODYNAMIC  MASS 


=  1,0*1  +  O.szi  n' 

0<fl'<  0.71 

(75A) 

-  i/(- 0.031 8+0.%1  ft') 

0,71  <  *\'<  3.1 

(75  B) 

DAMPING 

X* 

O 

|{ 

s? 

ll 

0<  rt'<  0.1 

(76  A) 

=  -0.061  +  0,71  n' 

O.l  <  K)'<  1.37 

(76  B) 

4  -O.HlS-R' 

=  1,515  Q 

1,37  <  »'<  3.1 

(76C) 

Heave 

HYDRODYNAMIC  MASS 

My 

-1,  is 

,-o.a?6 

0<  *'<  0,1 

0,1  <  H'<5,1 

(77A) 

(77B) 

*"•*38 


I 


I 
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DAMPING 

=  0,116  +  1.7 

0<  K)'<  OH 

(78A) 

-use-*’3"' 

0M< 

• 

(78B) 

Dimensionless  hydrodynamic  mass  and  wave  damping  coefficients  for  a  sphere 
based  upon  Kim's  study  along  with  the  above  approximate  curves  are  shown  in 
figure  11.  The  set  of  approximate  functions  will  be  used  in  the  computer  simu¬ 
lation  of  buoy  dynamics  for  a  spherical  buoy. 
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Birkoff  has  investigated  the  influence  of  body  symmetry  on  the  hydro- 

dynamic'  mass  dyadic.  An  oblate  spheroid  possessing  an  axis  of  symmetry  has 

five  hydrodynamic  mass  values  along  the  main  diagonal:  heave,  surge,  sway 

(same  as  surge),  pitch,  and  roll  (same  as  pitch).  For  an  ideal  fluid,  the  yaw 

45 

hydrodynamic  mass  is  zero.  However,  Lamb  studied  the  rotational  motion  of 
a  sphere  in  a  viscous  fluid  and  identifies  a  force  proportional  to  angular  accel¬ 
eration  that  can  be  considered  as  a  yaw  hydrodynamic  mass.  For  a  fully 
immersed  oblate  spheroid  in  an  ideal  fluid  with  the  centers  of  gravity  and 

pressure  coincident,  all  the  off-diagonal  terms  would  be  zero.  However,  for 

\ 

a  half-immersed  oblate  spheroid,  the  following  hydrodynamic  forces  are  coupled: 
Surge- Pitch 
Sway-Roll 
Pitch- Surge 
Roll- Sway. 

Thus  far,  hydrodynamic  mass  and  wave  damping  have  been  computed  for 
heave,  sway,  surge,  pitch,  and  roll.  Using  Lamb's  analysis  for  a  rotating 
sphere,  we  can  approximate  the  yaw  hydrodynamic  mass  for  an  oblate  spheroid 


3 
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w 


in  a  viscous  fluid: 


M*-  »  ^7 Tj>ftS> - LlA*. 

<X  3  •'  V  o  rt  —  .  <* 


3  "y  ”  /  + 
The  viscous  damping  in  yaw  is 


(79) 


1  ^  . (80) 

3  ^  1  +  Z£n  +  ZJZft 


K*  - 


where 


H  is  the  major  diameter 
J>  is  the  fluid  density 
M  is  the  viscosity 
* l /  is  the  kinematic  viscosity 
jS  is  defined  as  J 

Q~  is  the  angular  frequency  of  oscillatory  motion. 

Since  the  center  of  pressure  and  center  of  gravity  of  the  buoy  do  not 
necessarily  coincide,  force  components  due  to  hydrodynamic  mass  or  damping 
will  induce  moments  about  the  center  of  gravity.  The  projected  area  of  an 
oblate  spheroid  in  the  y-x  or  z-x  body  planes  is  a  semiellipse.  For  HD<.  0, 
i.e. ,  less  than  half  immersion,  the  center  of  pressure  is  located  a  distance 


HCP  3  jq  /(tf-Hp)3' 


(81) 


below  the  centroid  of  the  oblate  spheroid.  For  >  0,  i.e. ,  more  than  half 
immersion,  the  center  of  pressure  is  located  a  distance 


A 

* 


below  the  centroid. 


In  the  coordinate  system  shown,  a  positive  y  acceleration  will  induce  a 
negative  If  moment  and,  conversely,  a  positive  If  acceleration  will  induce  a 
negative  y  force  if  the  center  of  gravity  is  below  the  center  of  pressure.  Thus, 
the  coupled  roll- sway  hydrodynamic  mass  moment  of  inertia  is 

~  -^byr  ~  ~  M-hyy  (  HC6.~  Hcp)  ,  (83) 

Also,  a  positive  z  acceleration  will  induce  a  positive^  moment  and 
vice  versa.  Thus,  the  coupled  pitch-surge  moment  is 

~  -Mhy,  =  Mha  (H£S-  HcJ  .  (84) 


Summarizing  the  inertial  hydrodynamic  force  coefficients,  i.e. ,  the 
elements  of  the  hydrodynar  '.c  mass  dyadic  for  a  ^alf-immersed  oblate  spheroid, 
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I  Again,  for  the  special  case  of  a  sphere,  the  elements  of  the  dyadic  arc 


n/T,v  ~  (0.tZ6  +  1.7'*')  jxrft3  0<  n'<o.*f 

=>  (1.1%  £  013 *  )  jzcrft3  O.H<\ft'<3.1 

Y\yy  =  ~  0/  0<  fi'<  o,i 

-(-0.061 +o,7i  n')  j>cr  ft3  o.i<n'<i.37 

=  (1. 515  eat/ls‘n)  J)(T  ft3  1.37 <  ft '< 3. 1 

-  <? 


Klyy  —  flyy - Wyy  (  ~~  ^cp') 

~  ^i*Z  *=  ^Z2  (  “*  Hep)  » 

Dissipative  forces  due  to  viscosity  are  assumed  to  follow  a  velocity- 

squared  drag  law.  There  is  some  question  as  to  the  validity  of  representing 

an  unsteady  force  with  a  coefficient  based  upon  steady  flow  experimental 
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measurements.  Martin  has  found  that  the  mean  drag  coefficient  for  a  plate 

started  impulsively  from  rest  is  an  order  of  magnitude  greater  than  for  steady 

flow.  However,  there  is  no  general,  analytical  method  available  to  compute  the 

47 

viscous  forces  for  a  fully  turbulent,  oscillatory  flow.  Schlicting  cites  use  of 
a  method  of  successive  approximations  for  unsteady  laminar  flow.  Since  viscous 
forces,  f(a2),  are  an  order  of  magnitude  less  than  hydrodynamic  inertia  and 
wave  damping  forces,  f(a3),  and  because  steady  currents  are  acting  on  the 
buoy,  the  viscous  forces  are  assumed  to  follow  a  velocity- squared  drag  law  for 
subcritical  Reynold's  numbers. 

Analagous  to  the  hydrodynamic  mass  dyadic,  the  viscous  force  matrix  will 


contain  ten  elements.  The  viscous  force  coefficient  matrix  is 
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where 


0  o 


0  0 


0  0  0 


D  = 


0  0  0 


0  0 


0  0 

dkti  0 


0  0  0 


dxx  —  2.  ^dh  Ah 

Jyy  -  olZi ,  -  ^  Cps  As 


J  _  /  _  /  _  «  TV  rn3  3+ +  z .g V 

3"/**  t+£JS*+a,B*A* 

dy-ff  ~  dyy  —  “  dyy  (  flgp) 

ol^jg  -  djz  --  {  HC6.  ~  Hcp) 

The  projected  area  in  heave  is 


A„  S7Ti‘ 


and  in  surge  or  sway  is 


A5  =  f«i>  +  ^/WhT  +  «  b  sm'  A  , 

d.Q 

For  half-i:nmersed  spheroids,  Hoerner  “  shows  plots  of  drag  coefficient 
(surge  or  sway)  versus  the  ratio  of  vertical  to  horizontal  axis  dimensions. 
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Hoemer's  plot  for  a  subcritical  drag  coefficient  is  fitted  by  the  function 

CDS  ~  0.35 Y  ft/b  fa  <  10  j  subcritical 

for  various  spheroids.  The  drag  coefficient  in  heave  is  approximated  with  one- 

half  the  value  for  a  sphere.  C  "0.3  . 

Dn 

3.1.5  Ocean  Waves 

In  order  to  compute  the  magnitude  of  the  hydrodynamic  forces  with  the 
coefficients  just  derived,  the  relative  motion  of  the  fluid  surrounding  the  buoy, 
relative  to  the  buoy,  must  be  computed.  A  mathematical  model  of  the  sea  state 
must  be  developed. 

24 

The  following  hypothesis  offered  by  St.  Denis  and  Pierson  in  1953  has 
been  verified  by  Dalzell,49,  °°  Gerritsma,51  and  others  for  the  motions  of  ships 
in  a  random  seaway : 

1.  Assume  that  the  sea  can  be  represented  as  the  linear  sum  of 
elementary  waves,  each  traveling  in  the  manner  described  by  the 
classical  Airy  formulas  of  linearized  water  wave  theory.  Each 
component  wave  train  will  have  random  phase. 

2.  Weight  the  compon^  *  waves  to  have  the  same  spectral 
characteristics  as  the  observed  sea  state. 

3.  Assume  that  the  body  response  to  a  random  sea  is  the  sum 
of  its  responses  to  the  various  frequency  components. 

Within  the  constraints  imposed  by  the  assumptions  made  in  the  derivation  of 
the  hydrodynamic  coefficients,  i.e. ,  body  dimensions  are  small  compared  with 


a  wavelength,  the  St.  Denis- Pierson  hypothesis  should  be  better  suited  to  the 
case  of  a  buoy,  with  dimensions  on  the  order  of  10  ft,  than  to  ships,  with 
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dimensions  on  the  order  of  100  ft.  (For  example,  as  a  worst  case,  a  10— ft 
buoy  in  waves  with  100— ft  wavelengths  would  only  cause  a  peak  error  of 
3. 3  percent  in  the  computed  elevation  of  the  mean  waterplane.  This  error  will 
decrease  as  the  wavelengths  become  longer.) 

From  the  Airy  formulas,  water  particle  motions  in  deep  water  for  waves 
traveling  along  the  z  axis  are  described  as  follows: 

Wave  Height 

Xw  -  SIN  (kz  -  cr  t)  (90A) 

Vertical  Velocity  Component 

Xw=~'B!^i  C0S(kH-<rt)  (90B> 

Horizontal  Velocity  Component 

Aw  —  ~Crt)  (90C) 

Vertical  Acceleration  Component 

Xw  =•  SiN(kz  -<rt)  (90 d) 

Horizontal  Acceleration  Component 

~  C0$(kz~<rt)  (90S) 
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Wave  Slope 

JBW  -  -TAN](nwk  cosCka-crt^ 
=  -  cos(  k&  -crt) 

Angular  Velocity  of  Free  Surface 


A 


j3w  £  -nwkcr  5/m(Kz- cr t) 

Angular  Acceleration  of  Free  Surface 


(90  F) 


(90G) 


(90  H) 


where 


JSW  k  (T* cos(kTL-cr1>) 

K  is  the  wave  number  (  k  ~  Z^/l-o) 

La  is  the  wavelength 
CT  is  the  angular  frequency 
£  is  the  time. 

The  assumption  that  body  dimensions  are  small. compared  with  wavelengths 
implies  that  body  displacements  are  small  compared  with  wavelengths.  Thus, 
we  find  that 


3  SIN  Crt  > 


K  C0S  > 


(91A) 


(91 B) 


Z.w  &  -^wcr  SIN  crt 


(91C) 
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Xw  s  crz  sin  crt  ) 


(DID) 


£  ~nw(Tz  COS  crt  , 


(DIE) 


>SW  =  >V  K  COS  crt  , 


(DIF) 


^  =  -'Wwkcr  SIM  <rt  J 


(DIG) 


,3W  «  -n„k<rzco$crt  t 


(91 H) 


Transforming  the  water  mass  velocities  and  accelerations  to  body  coordinates, 
we  see  that  the  velocities  are 


-  n 


and  that  the  accelerations  are 


(92A) 


XWz 


(92  B) 
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The  location  of  the  intersection  of  the  free  surface  and  the  "vertical"  axis 


of  the  buoy  is 


Hfi  =  CoSJB  COS  r  -f  -  Xw  C0SJ2  COS  Y 


In  order  to  develop  a  random  wave  model,  the  statistical  properties  of 
the  sea  state  must  be  described.  The  development  of  wind  waves  on  a  body  of 
water  is  either  limited  by  the  distance  to  land  in  the  direction  from  which  the 
wind  blows  (fetch  limited)  or  by  the  length  of  time  during  which  the  wind  acts  on 

52 

the  water  surface  (duration  limited).  For  the  fetch  limited  case,  Bretschneider 

53 

normalized  the  original  wind  wave  forecasting  relations  of  Sverdrup  and  Munk 

and  included  much  additional  data.  Bretschneider's  dimensionless  curves  have 
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been  approximated  by  piecewise  linear  functions  by  Pattoa  The  approximate 
functions  are  as  follows: 

Significant  Period 


=  2, 17  xio 


F<  2.0*10  (94A) 


LOG-  n  -  -L136  +  0.ZZ3  Log  F  £*ios<  F<  1.2  *!05  (94B) 

'3 


r  =  2.o 

<3 


F>!,Z*\0*  (94C) 


Significant  Height 


Hi,  -S.7HXI0' 


F<t.S  K  iO  (95A) 
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LOG-  H,  =■  -Z.5  t-O.ilS  LOG-  F  L6xio'z<F<  5.0x10*  (95B) 

H.  =  £.?£*  /(?"'  F  >5,0x10*  (95C) 
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Weigel  extends  Bretschneider’s  dimensionless  curves  to  the  case  of 
duration  limited  wind  waves.  These  curves  are  approximated  by 
Significant  Period 


-  -  as  +  asav  log 


> 


(96) 


Significant  Height 


LOG  Hys  = -2.272+  C7, 32 ?2  iO&D  ; 


(97) 


where 


Mj  ZttW> 


u  -  k 
WVZ ) 

f  ~  JlEL 
r  W*> 


sD 

w  > 


the  dimensionless  significant  period 


the  dimensionless  significant  height 


the  fetch  parameter 


the  duration  parameter 
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77  is  the  significant  period 

4 


is  the  significant  wave  height 


f-  is  the  fetch 

J)  is  the  duration 

Cj  is  the  gravitational  constant 

w  is  the  wind  velocity 

In  order  to  determine  if  the  waves  are  fetch  or  duration  limited,  the 
minimum  duration  for  a  given  fetch  must  be  determined.  From  Bretschneider’s 
curves,  the  minimum  duration  is  given  by 


LOG 


=  1,HTT-  0.255  LOG  F 


(98) 


If  the  duration  (D)  is  less  than  the  minimum  duration  (D^^)  the  waves  are 

duration  limited,:  if  greater,  the  waves  are  fetch  limited. 

For  a  given  wind  speed,  fetch,  and  duration,  the  significant  wave  heighfc 

and  period  canbe  computed  byusingthe  above  equations.  Longuet-Higgins56  and 
52 

Bretschneider  have  shown  that  the  distributions  of  wave  heights  and  squared 
periods  can  be  represented  by  Rayleigh  distributions.  Thus,  the  wave  height 


64 


distribution  is 


and  the  wave  period  distribution  is 


rp3  -0.67s($r)'t 

p(t)  =  2.  7-fae 


The  mean  wave  height  and  periods  are 


H  =  o.  6Z5  Hys 

and 

T  =  o.7SSTh 

Ocean  wind  wave  amplitude  spectra  are  of  the  form 


(99) 


(100) 


(101) 


(102) 


(u )  =  A  eSn  (W3) 

# 
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Pierson  and  Moskowitz  analyzed  54  data  sets  and  evaluated  A,  B,  m,  and  n  . 

The  exponent  m  was  set  equal  to  -5  and  n  was  found  to  vary  between  2  and  4 
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depending  on  the  wind  speed.  Kcttler  optimized  the  parameters  by  using  a 
least-squares  fit  of  each  data  set  used  by  Pierson-Moskowitz  and  proposed  the 


form 
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where  £  (ft>)  is  a  random  variable  whose  values  are  equally  probable  for  any 
value  between  0  and  2  rr .  The  integral  is  not  an  integral  in  the  Riemann 
sense  since  the  function  is  discontinuous  because  of  the  random  variable.  The 
expression  indicates  that  the  random  wave  heights  can  be  represented  as  a 
finite  number  of  cosine  waves  of  different  frequency,  each  having  random  phase. 
Let  the  spectral  density  (co)  be  partitioned  into  K  frequency  bands 

(figure  12)  in  such  a  manner  that 

0  <  cv0  <  OJ,  <  Cd^  •  •  .  <  =  Wy 

where  the  spectral  density  is  essentially  zero  if  W  is  greater  than  DO* . 

The  width  of  the  nth  band  is 

A  6Uh  ~  C0H_,  (107) 

and  the  mean  angular  frequency  of  the  nth  band  is 

"n  =  2  (H  +  W„_, )  ,  (108) 

The  random  wave  height  model  is  now  defined  as 

£  _ _ 

Xw  ( =  aW  Cos(wht  +  €n)  } 

h*i 

where  ^  K)  ==  £,  3,  .  •  •  ^  jN"  are  independent  random 
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variables  distributed  uniformly  over  the  interval  0  to  2  TT  .  Borgman 
indicates  that  an  equally  spaced  subdivision  A  OJn  —  foX/N*  will 
result  in  Xw^t)  repeating  itself  with  period  &rr/oj^  .  Borgman 


Figure  12.  Borgman's  Energy  Partitioning  in  a  Bretschneider  Spectrum 


partitions  the  spectrum  on  an  energy  basis  using  the  cumulative  spectrum  given 


by 


(A) 


r 

$Hz  ((*>)  ~  t  SHz  M  da) 

*  J 


(109 


thus, 


- —  £  (fSHz  ((M>)  y 


(110 


hence, 


a 


*w(t)  =  jf  2Js^)-E„mS  cos  (w„t  +eh)  f 


(111 


n=i 


The  periodicity  is  avoided  if  the  set  of  l*)n  values  are  chosen  to  make 
SHz(M»)  -  Shz(m  constant  for  all  n  .  Let 


SHz(°^)  ~  SH*(0On„,)  =  h  $ 


(Hi 


then,  the  instantaneous  wave  height  is 


Xw(fc)  =  ^^C0S(wnfr+O, 

«■»! 

The  frequencies,  COn  ,  are  given  by 


(ii; 
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which  corresponds  to  an  equal  subdivision  of  the  energy  coordinate  axis  for 

(w)  ' 

Using  the  Bretschneider- Pierson  spectral  density  of  the  form 


C  ,  a  AB 


which  is  directly  integrable,  we  find  that 


.w 


-e^Js 

"0 

Integrating  out  to  infinity,  we  find  that 


A  pB/“H 

a  r 


A 

a 


and 


B 


w"  ~  ' wST(w/„) 


where  the  coefficient  is  given  by 


(115) 


(116) 


(117) 


w/jy , 

Next,  we  define  the  various  elements  of  the  random  sea  state  model: 
Wave  Height 


(118) 


IT 

xjt sm ( a„ t  ■+  e„) 


(119) 
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Vertical  Velocity  Component 


Aw(t)  =  g  -MnfgHt(U„)AWn'  COS (w„t  +  £„)  (120) 

fl-l 


Horizontal  Velocity  Component 


KwCt)-  2.  SIN  (faJn  £  +  €M) 

0=1 


Vertical  Acceleration  Component 


(121) 


xw(t)  -  i^Js^aMul  siN(u„t+en) 


h*»  i 


Horizontal  Acceleration  Component 


.W 


Hw(t)  =  I  ^-w„y^8(w„)iWn  C05(w„t+e„) 

h^i 


Wave  Slope 


if 

A(t>=  cos(unt+e„) 

0SJ 


Angular  Velocity  of  Free  Surface 


Angular  Acceleration  of  Free  Surface 


/£  _ | 

A/t)=  w„  cos(u»t+en) 


(122) 


(123) 


(124) 


Aft)  =  §;  4  -k„ w„/5Hl!(®4«r  S,N Cw»4 +£»)  <125) 

h=i  V 


(126) 
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•V'  Again,  the  wave  height,  velocity  components,  and  acceleration  components 

must  be  transformed  to  body  coordinates  according  to  equations  (92A),  (92B), 
and  (93). 

A  ten-component,  random  sea  state  model  was  programmed  in  FORTRAN 
for  use  with  the  UNIVAC  1108  digital  computer  at  the  Naval  Underwater  Systems 
Center,  New  London  Laboratory.  This  program  is  shown  in  appendix  B  as 
subroutine  "RWAVE. "  The  ten-component  model  was  found  to  truncate  the 
low  and  high  ends  of  the  spectrum  (the  first  and  ninth  frequency  components). 
Since  the  low-frequency  end  of  the  spectrum  is  important  in  the  computation 
of  the  response  of  mechanical  systems,  Borgman's  frequency  partition  method 
was  modified  to  include  the  low-frequency  energy  by  the  following  numerical 
scheme; 

1.  Use  a  trial  and  error  method  to  compute  the  frequency  at 
which  the  value  of  the  spectral  density  rises  above  some  threshold  value, 
for  example,  Sy*  (b>)  >0iOi  ft?  .  Let  this  frequency  be  denoted 

C^o 

2.  Thus, 

00 1  ~  Mo) 

and 

A(0,  =  00,  -U)0  , 

3.  Use  the  energy  computed  for  the  second  partition  Syzity  (tOi) 
to  compute  the  component  amplitude. 

This  method  accounts  for  the  energy  in  the  low-frequency  end  of  the 


spectrum. 


~XW  COSjS  cosy  —  A0&  COS^S  C(?sy  J  (127) 

where  X0&  is  the  vertical  height  of  the  center  of  gravity  of  the  buoy  below  the 
mean  waterplane  with  no  buoy  pitch  or  roll.  The  transformed  wave  height  must 
be  subtracted  from  the  heave  motion  of  the  buoy  in  order  to  apply  equation  (93), 
the  location  of  the  intersection  of  the  waterplane  and  the  "vertical"  axis  of  the 
buoy.  Restricting  the  study  to  waves  traveling  along  the  x  axis  of  the  coordi¬ 
nate  system,  the  wave  slope  y8w  in  the  inertial  coordinate  system  is  given  by 
equation  (91F),  for  sinusoidal  waves  and  by  equation  (124)  for  the  random  wave 
model.  The  attitude  of  the  waterplane,  relative  to  the  buoy  coordinates,  is 
described  by  the  angles  jQs  and  ^  .  The  slope  of  the 

sea  surface  is 


X  ^  cos' {cos  Xs  COS  a) 


(128) 
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<g> 
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| 

I 

K’ 
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g.  f 

Ift 


With  the  location  (equation  (93))  and  slope  (equation  (128)),  the  buoyant  forces 
and  moments  for  a  sphere  are  given  by  equations  (49),  (55),  and  (56). 

In  a  similar  manner,  the  hydrodynamic  forces  are  computed  by  consid¬ 
ering  water  mass  movements  relative  to  the  buoy.  With  the  assumption  that  the 
buoy  dimensions  are  small  relative  to  the  wavelength  and  applying  the  Froude- 
Krvloff  hypothesis*  for  the  viscous  forces,  we  find  that  water  mass  velocities 
and  accelerations  are  given  by  equations  (91B),  (91C),  (91D),  (91E),  (91 G),  and 
(91H)  for  a  sinusoidal  wave  model  and  by  equations  (120),  (121),  (122),  (123), 
(125),  and  (126)  for  a  random  wave  model.  Water  mass  accelerations  are 
simply  transformed  to  buoy  coordinates  (equation  (100 B))  and  subtracted  from 
buoy  accelerations  to  compute  the  motion  of  the  buoy  relative  to  the  water  mass. 
The  relative  acceleration  is 


Q 


/ 


Q-fl-Q 


W 


(129) 


the  acceleration  vector  relative  to  water  mass 
the  buoy  acceleration 
the  water  mass  acceleration. 

Water  mass  velocities  are  more  complex  because  ocean  currents  exist  and  must 
be  included.  There  is  a  difference  in  reference  frame  that  must  be  resolved 
since  buoy  velocities  cause  a  reactive  force  (i.e. ,  a  positive  buoy  velocity 


where 


Q  is 

Q  is 

4- 


*The  presence  of  the  buoy  does  not  appreciably  change  water  mass  move¬ 
ments  in  the  vicinity  of  the  buoy. 
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causes  a  negative  force)  while  waves  and  currents  cause  active  forces  (positive 
velocities  cause  positive  forces).  The  velocity  vector  of  the  buoy,  relative  to 
the  water  mass,  is  given  by 


0  —  Q--n*Qw-  fl'Qc  > 


(130) 


where 

# 

Q  is  the  buoy  velocity  vector 

# 

Qw  is  the  particle  velocity  vect  r  due  to  waves 
Qc  is  the  current  velocity  veer  ?  c. 

For  current  velocity  components  C  and  CW  acting  in  the  positive  y^  and 
directions,  respectively,  the  co  gent  velocity  vector  is 


o 

cv 

cw 


(131) 


For  viscous  dissipative  forces,  which  are  functions  of  the  velocity  squared,  the 
velocities  are 


The  absolute  value  is  used  in  order  to  preserve  the  sign  convention. 

If  the  buoy  were  free  floating,  not  moored,  the  equations  of  motion  for  the 
buoy  could  now  be  integrated  for  a  given  set  of  environmental  conditions  to  solve 


for  the  buoy  response. 


'Y^y^S-SSSWBBW, 
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The  equations  of  motion  for  a  spherical  buoy  are  summarized.  Equation  (5)  is 


M$  -MG  -B  -  H  -V.-'  -T. 

In  the  R  coordinate  system,  the  heave  acceleration  is  given  by 

x  =  ft;  [  m  g  cos  yens  JR  -  B  cos  r  cos,£ 

+  7^  cosrcosjB  +Trsihnr  +  T*  cos  f  sinjS 


The  sway  acceleration  is  given  by 

y  =  rnfm5  SIN  y  CoS'g  -  Bsmrcosja  -  mhyy-  y' 
-mhy/  y'  -  rv  y'  -  nYr-  f'  -  <An- y"(y'| 

-  •  ?*l  t"l  -wys,-%  sm  v  COS Ji 

+T J,  cos  y  -  Te  s/m  r  sm jb  1 


The  surge  acceleration  is  given  by 


=■  ^  [-mg  smfi  +  £  sinJI  -  mhw- z'  ~ 

~  nzj  ’J'  -  £"l 


~  Tr  +  %  COSJ 


1. 
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<  * 


The  yaw  acceleration  is  given  by 

*  -  x.  f~  ^ '  <*  -  <4«-  <*  /<*/  ]. 

The  pitch  acceleration  is  given  by 

J  -  -f;  f-^  -  »V  * 

-W.£ 

-Hm  ("71  -sw,8  +  7;  cos^e)] ' 


The  roll  acceleration  is  given  by 

*  =  Tm  [~Mv  ~  Kr'*  ~  "V*  y'~  rV*' 

df/-X"ly1  -W* 


(~T!siN}rcoS'fe+TyCos7r-TxsiNXsiNj>) 

« 

Hydrodynamic  inertias  and  wave  damping  are  computed  for  motions  of  the  buoy 


relative  to  the  water  mass. 

Thus,  for  waves  traveling  along  the  z  axis,  ’ 

see  that 

•»  t  •#  ,, 

J 

X  =  X  -  x„t 

K'  =  X-KVk 

9'-9 

Y'~y 

z'  =  z  - 

Z'  =  Z-Zvt 

«  /  «. 

(X  S(X 

#  ^  # 

CX  -  (PC 

oc  =■  cx 

*1 

A*J~A 

J3'~JS- X, 

f  = 

r-t 

■V 

tt 

°< 

V 

\ 


77 


where  the  subscript  W&  indicates  that  the  computed  wave  induced  water 
particle  motions  are  transformed  to  the  buoy  coordinates.  The  viscous  drag 
forces  and  moments  must  include  a  contribution  due  to  the  steady-state  currents; 
thus, 


k"  =  A ' 

y y'-  cv 


«  //  • 

CK  =  CK 


/ 


The  draft  of  the  buoy  is  H0  =*  Hc<-  )('  where  f-f  is  the 

height  of  the  center  of  gravity  of  the  buoy  from  the  bottom  of  the  buoy  hull.  The 
buoyant  force  is 

Bv>§tt  +  ife  cosj cos3#' ]' 

The  pitch  and  roll  buoyant  moments  are  given  by 

=  £■  (zce  costs'  +  ()is  -  xj  srn# ) 

and 

=  cosy'  +(HC5.- xce)siNt") , 

The  hydrodynamic  inertia  coefficients  are  summarized  in  equation  (85).  Wave 
damping  coefficients  are  summarized  in  equation  (86).  Viscous  drag  force 
coefficients  are  summarized  in  equation  (87).  Finally,  wind  forces  are  shown 


in  equation  (18). 
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With  all  the  coefficients  in  the  equations  of  motion  defined  and  with  the  sea 

» 

state  model  available  as  a  forcing  function,  the  dynamics  of  a  free-floating 
buoy*  (no  mooring  line)  can  be  investigated  by  numerical  solution  of  equation  (5). 
The  matrix  equations  of  motion  can  be  rewritten  to  yield  six,  coupled,  second- 
order,  ordinary  differential  equations  with  nonconstant  coefficients  —  three 
force  equations  and  three  moment  equations.  They  can  be  rewritten  as  twelve 
first-order  equations. 

These  equations  can  be  solved  numerically  in  the  time  domain  for  the  six 
buoy  motions  using  a  standard  higher  order  numerical  integration  algorithm 
(usually  Milne  Predictor- Corrector  methods  or  Runge-Kutta  methods).  In  this 
study,  a  fourth-order,  Runge-Kutta  methoi.  was  used  to  integrate  the  buoy  equa¬ 
tions  of  motion  on  a  UNIVAC  1108  digital  computer.  Although  a  bit  slower  than 
Predictor- Corrector  methods,  the  Gill's  Meuiod  Runge-Kutta  subroutine  is 
available  on  the  UNIVAC  1108  and  was  found  to  be  easier  to  implement  for  large 
numbers  of  coupled,  second- order  differential  equations. 

g 

Errors  in  the  fourth-order,  Runge-Kutta  method  are  on  the  order  of  h  , 

where  h  is  the  integration  step  size.  For  the  smaller  oceanographic  buoys,  it 

_2 

was  found  that  a  step  size  on  the  order  of  10  sec  was  required  for  numerical 

-3 

stability.  If  a  step  size  of  5  x  10  were  used,  the  error  would  be  on  the  order 
of  3. 125  x  10  12. 

The  computer  programs  used  are  described  in  detail  in  appendix  B, 

*The  mooring  line  tension  force  will  be  developed  in  the  next  section  on 
cable  dynamics  and  the  equations  of  motion  for  the  buoy  will  be  coupled  with 
those  for  the  cable. 
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O  3 . 2  Cable  Dynamics 

A  study  of  the  dynamics  of  cables  in  the  ocean  environment  must  include 
hydrodynamic  forces  acting  on  the  cable  and  be  capable  of  including  nonlinear 
stress- strain  properties  for  the  cable.  Since  oceanic  currents  may  flow  in 
various  directions  at  different  depths  and  the  excitation  at  the  ends  of  the  cable 
may  be  described  by  a  three-dimensional  force  vector,  the  analysis  must  be 
conducted  in  three  dimensions. 

3.2.1  The  Cable  Equations 

Consider  a  free  body  of  a  cable  segment  of  length  dS^  (figure  13)  being 
acted  upon  by  an  external  force  per  unit  length,  Q  ,  and  assume  perfect  flexi¬ 
bility.  The  geometric  center  of  the  cable  segment  lies  at  a  point  (xq  ,  yQ ,  z q) 
in  a  cartesian  coordinate  system.  The  external  force  Q  dS^  acting  on  the  cable 
segment  can  be  resolved  into  components;  they  are 

Q  dS0  ~  X  dSo  *  I  +  Y"  JS$  '  $  +  %  clS0  •  k  ,  (133> 

Writing  Newton's  Second  Law,  we  see  that  in  the  x  direction 

yU (S,) AS,  f§*=X  AS'  +  [f(&+AS.)- ?©,)] -I  .  <134A> 

in  the  y  direction 

yU(S0)  d$0  |p  —  Y +  [T  (So  +  dSj-TfSlJ'J  .  (134B) 


and  in  the  z  direction, 
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The  unit  vector  tangent  to  the  cable  is 


_ _ 'IS, _ 

/fiRslW 


At  point  ,  the  tension  T  is  tangential  to  the  cable  if  the  cable  is  assumed 
to  be  perfectly  flexible  (i.e. ,  the  cable  cannot  support  bending  moments).  The 
tension  is 


T  «  T 


where  T  is  the  magnitude  of  the  tension.  The  x  component  of  the  change  of 
tension  along  the  cable  is 


a  /lv/\  2 


II.  ? 


3  (JL.M. 

Sllgl  35‘ 


Substituting  equation  (136)  into  equation  (135A),  we  find  that 


(137A) 


and  in  a  similar  manner  for  the  y  and  z  forces,  we  find  that 


/*($,) 


dey  __ 


-  y  +  f  X- *  iZ 
~  \  |iS]  ^ 


(137B) 


and 
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(137C) 


The  strain,  (£ )  ,  is  defined  as 


c6S ~~  <^5g  _  |3r  I  j 

€  ~  us.  -ter1  ) 


<W8) 


where  S  is  the  strained  length  of  the  cable.  Substituting  equation  <138)  into 
equations  (13 7A).  (137B).  and  (137C),  we  can  write  the  equations  of  motion  as 
follows: 


Ip  “  h,(£e X  m  0  >  (18»A> 


(139B) 


JftG+e  *  l&)~  *  ”  0  . 


In  addition,  the  auxiliary  equation  (138)  becomes 

£  =  —  -  1 
C  JSe  1  . 

A  constitutive  relation, 


T  *  T(e)  } 


(139C) 


(140) 
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must  also  be  defined. 

This  set  of  equations  is  standard  for  cable  systems  and  has  been  devel¬ 
oped  by  Cristecu, ****  Schramf1  Whicker,**2  Lindsay,**2  and  others. 

It  is  desirable  to  transform  the  equations  of  motion  to  a  "natural"  coor¬ 
dinate  system  (i.  e. .  a  coordinate  system  aligned  normal  and  tangential  to  the 
cable),  because  hydrodynamic  forces  are  usually  defined  as  normal  and  tangen¬ 
tial  to  the  cable.  The  coordinate  systems  are  shown  in  figure  14.  Rotating 
about  the  y  axis,  the  transform  is 
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Rotating  about  the 

z  axis,  the  transform  is 

I 

I 

[' 

7 

1 

~C0S(p 

SlNp 

0 

i 

f 

f. 

Aa  - 

-S\N‘(p 

cosp 

0 

tt 

1 

0 

0 

1 

m 

for  both  rotations,  w'  'hat 

A  =  A2- A1  5 


thus, 


A  = 


cox  p  cose 
-sin $  cose 


SIN  Ip  COSpSlNd' 
c OS#  -SIN  p  SIN e 

0  cose 


(141) 
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A  is  the  transform  matrix  from  the  x ,  y ,  z  coordinate  system  to  the  x" ,  y" , 
z"  system.  Since  A  is  orthogonal,  we  see  that  A  1==At;  thus, 


A’1  = 


costpcose  -sin  <p  cose  -sino 

$<Ntp 

cos  <p  sin  e  - $iN<t>  sin  e  cos  e 


(142) 


In  the  double-primed  system,  let  x”  be  along  the  cable  and  z"  and  y"  normal 
to  each  other  and  the  cable.  Since  ^  is  a  unit  vector  tangent  to 


the  cable, 


-  {£) 
<s  KzsJ 


transforms  to 


1+6  \ 

coordinate  sj'stem.  This  transformation  is 


in  the  double-primed 


/ 

1+6 


(Mho'll]  , 


or 


(ih) 


5X  ^ 

r 

dse 

C 

ay 

3* 

3z 

to. 

SlN  <p 
C OS(p  SIN  9 


(143) 


The  components  of  equation  (143)  are 


3X 


“  (1  +  COS  (j)  COS&  ) 


§|  a  (1  +  6)  SIN  , 


(144A) 


(144B) 


I j  ss  (i+£)  cos  (p  SIN  e  , 


(144C) 


Lot  H,  G,  and  I  be  forces  per  unit  length  acting  on  the  cable  along  the  x”,  y", 
and  ■/."  axes,  respectively.  In  earth  coordinates,  the  cable  forces  per  unit 
length  are 


xl 

»  - 

H 

y  =  a" 

Q 

2 

U  J 

I 

a  - 

(145) 


Writing  equations  (139A),  (139B),  and  (139C)  in  vector  form  (in  the  unprimed 
coordinate  system),  we  see  that 


(T cos (p  cose ) 


'a  / 


(tsincp)  +/\"  G  -  yuCs,) 


Its  (TcostsiNe) 


(146) 


where  u,  v,  and  w  are  velocities  in  the  x.  y,  and  z  directions,  respectively. 
Transforming  to  the  double-primed  coordinate  system,  equation  (146)  becomes 


f^(T  cos <p cose)  H 

A  | ^(TS)N(p)  -h  Q 

(Tcospme)  i 


-  /a(&)A  % 


The  first  term  on  the  left  is 


|I  (cos*#cos*e  +Sin*$+cos*<psiNze ) 

+T||  (-CO$*0COS<pSlN<p  +SlN4cO5(j>~StH*eSlN<pC0Slp) 
+rp||  (-COS^SlN&COSe*  C OS*<f>SlU0  cose) 

^£(-cosa©  siufcosp  +  COSpSIN#  -  Sin*S  sin#  cos  #) 
+T  H  (  COS2©  SIN*#  +  COS*#  +  SINZ9  SIN11#) 

+T i&(siN<p cose swe cosep  - sin#  cose smd  cos#') 

(-sihg  cos#cose+siNd  cose  cos#) 

+Tl|  (cose  sinosin#-  cose  sin 0  s/N#) 

vSq 

+rp||  (c os#siN*e  +  cos#  cos* o) 


which  can  be  reduced  to 


f  “*  %  • 

Now  consider  the  term  on  the  right  side  of  equation  (147), 


yU(S.)A 


89 


If  U,  V,  and  W  are  the  velocity  components  along  the  x",  y”,  and  z"  axes, 
respectively,  the  velocity  components  in  inertial  coordinates  are 


’ttl  \v 

V  =  A~'  v 

w  w 


then,  the  acceleration  components  are 


Multiply  through  by  the  mass  density  and  the  transformation  matrix  (A), 


3U 

V 

+/<(s,)AA’' 

yU(S.)A 

5t 

=  /t6s.)A^M") 

v 

i  J 

v 

. 

(149) 

Considering  the  first  term  on  the  right  of  equation  (149),  we  see  that 

C0S<p  cose  -Siscos  9  •‘SIN  9 

ft  (A  )  =■  ft  394 1  0 

CCS(f>SINd  -SlN<pSlN9  COSO 


or 


■  — '  V*  3  ‘  'J>K  "5  J--  “ 
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f  .jse&ttitty  \ 
\  -cospsm^J 


( 


-cose  cost  ||  'N 

+S(M0Wg  J 


(- «*»  n) 


(coi^)  (-*»*%)  0 

(-sine sint&  \  /-sine  corfU  'N  /  Km *  &) 

[  -¥COStcOJd^)  \  -SIH+Coselg  )  \sm6&) 


A  further  expansion  of 


A  ft  (A-1) 


gives 


(^C0S*eCOS<)S/M^  f-coshcos*#^ 

j  -cos*toosesine\l^  J+SlNpSINBCOSt003^^ 

■<  +StH<f)COSp  |f 

|  -SJW*0CO$*<£  || 

(_  +COS*<flCOSd  SlNd  ^  ^-SW$*W0  WJ$coj0  ^| 


C0Sa<9  swa(^  ^  r CO$*0  SIA/f  ««$>  ||  'I 

+SMtSlN9GOs9COst^i  I  J-5IN*PC4T0  SlN$ 

<  +655*  ♦  $  c>< -cos^sw*^ 

I+Scwtyjw^ll  ( 


\j-Smtsin9cosdcosp  ^j|  j  ^+$lNg4s/U$C0se  ^ J 


I -costcos*^  I 
)  ~a>stsiNe0$r 


coszesnj(f  U 

+Slti*d3)Ht\§. 

3C 


rSIMdC0S,95M^^' 

J  +'swa0  COJ^  i| 

-Sintco3$sin9  j 

+CO3*0Q>54>  | 


r SIN0  cosecos^ 

-sin2©s;n^^  ‘ 
-SiH'JccSGcospl 
- C0Sa6  sin<j>}£ i 


S|N0COS0  3| 
}p 

-sine  coret* 

it 


The  above  matrix  can  be  reduced,  and  we  find  that 


91 


Substituting  into  equation  (149),  the  accelerations  become 


(151) 


Substituting  equations  (148)  and  (151)  into  equation  (147)  and  writing  the  compo¬ 
nent  equations,  we  obtain  the  equations  of  motion: 


S  +H  -  -  V*4-Wcos4  If)  (152) 

T||e-t-Q  -  /«&)(!£  +  1/ II  +Vswt  n) 


(153) 


92 


Tcos^H  *  I  **/*&)($  +l/cos0|f  ~V’sw0|f)  ^ 


In  addition,  we  have  the  strain  definition, 


H  ~l+e 


and  the  constitutive  relation, 


T  =  TY«) 


(156) 


Equations  (152)  through  (156)  give  five  equations  in  eight  unknowns 
(T}4>,6>U')V)W’i3l€)t  Three  more  geometrical  equations  can  be  developed  by 
considering  the  velocities.  The  velocities  are 


V  -  A  V 

w  V 

L  J  L  J  ) 

and  the  space  derivative  of  the  velocity  vector  is 


i-  V 


|  (A"')  V  +  A'  % 
V  w 


Consider  the  left-hand  term  of  equation  (157).  The  space  derivative  of  the 
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velocity,  vector  becomes 
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From  equations  (144A),  (144B),  and  (1440,  we  see  that 
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(l+e)  cosip  sine 
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Rewrite  equation  (157)  and  transform  with  A;  thus 
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(159) 


By  analcgy  with  the  development  of  equation  (148),  the  first  term  on  the  right- 


hand  side  of  the  above  equation  is 


'  'tat 


-  vi  W-  ^■fcTffg*g*qS»  g*  * 


&*)  u 


(It 6) COS <p  If 

OS’ 


By  analogy  with  the  development  of  equation  (150),  the  second  term  on  the  right- 
hand  of  equation  (139)  is 

-&  -cos*!i 

fs.  0  sm*fs. 


[co^H  -sw*33. 


Substituting  into  equation  (159),  and  writing  the  component  equations  gives 


(160A) 


(160B) 


W  =  O+£)cos<^5l“(u'C0^li'V's'Nli’i,) 


J  ,  (xdOC> 


Rewriting  the  constitutive  relation,  we  obtain  the  time  derivative  of  tension: 

ar  __  ctr  36_ 
at  9*t  . 


Finally,  the  equations  are  summarized  as  follows: 


g  -  -  V  -  W  cost ■£)+}{  =  0  t  (ici) 

T||o-^&)(i+l/«+v^w^)+(J=0  >  <162) 

T  cos*  ^ t  t/cos*  -  V W|§)+ 1=0,  ,i63) 

+  I)  =  °> 

g  -  «♦«)$  +(V||  +W*«*  ||3=  0  ,  (165) 


3W 

33. 


(l+6)co.-.^H  +(I/GW$||- Vsw*|§3  =0 


(166) 


3T  _  dl  M  «  0 

Jt  Me  H  u  <167> 

We  have  seven,  nonlinear,  partial  differential  equations  in  seven  unknowns. 

Using  the  constitutive  relation  equation  (167),  those  reduce  to  six  equations  with 

G 1 

six  unknowns.  Jeffrey  and  Taniuti  show  that  an  analytic  solution  in  the  form 
of  a  power  series  for  a  single,  nonlinear,  hyperbolic,  partial  differential  equa¬ 
tion  can  only  be  obtained  locally,  or  "in  the  small"  for  a  point  P  on  a  nonchar¬ 
acteristic  curve  T  in  the  time-space  plane.  Physically,  it  is  known  that  the 


it,'  vJ^V^-V^  :C4Ji2£vS 


»-vnfcyM.  — -II  v,'«m»^ 
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flexible  cable  can  propagate  both  tensile  stress  waves  and  transverse  flexural 
waves.  It  is  also  known  that  the  characteristic  velocity  of  the  transverse  waves 
is  a  function  of  the  state  of  stress  of  the  cable.  Thus,  the  usual  method  of 
characteristics  approach  to  solve  nonlinear  wave  equations  is  not  tractable 
because  the  characteristics  diverge  in  the  time-space  domain.  Solutions  m  ist 
be  obtained  simultaneously  at  the  same  locations  m  the  time-space  domain  in 
order  to  continually  update  the  transverse  wave  characteristic.  A  numerical 
method,  which  is  an  extension  of  Hartree’s  "hybrid"  method,  '  will  be 
employed  to  conduct  simultaneous  integrations  at  nodal  points  of  a  rectangular 
grid  in  the  time-space  domain.  Because  this  method  utilizes  integration  along 
characteristics  in  the  immediate  vicinity  of  the  time-space  grid  nodal  points,  it 
is  necessary  to  rewrite  the  cable  equations  in  their  "normal"  form,  i.e. , 
characteristic  equations. 

V’irst,  the  characteristics  must  be  found.  Assume  a  linear  variation  of 
tension  with  strain.  Equation  (167)  reduces  to 

r-f  e  -ks  . 

Rewriting  equations  (161)  through  (166)  in  the  form  +  +  C  ~  0  0-68) 
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Rewriting  equation  (168),  we  see  that 

-I  .  -r  r  A-*  T~N  T  r  1  ~l 


or 


A  AUt+ABUs+A  C  -0  > 


14  -M~'Bl/s  +  A"C  =  0 


1 


(109) 


The  six  characteristics  are  derived  from  equation  (169).  This  derivation  is 
shown  in  detail  in  appendix  C.  From  appendix  C  the  six  equations  (equations 
(161)  through  (166))  yield  six  characteristics; 


^  =r  +  /_J _  dl. 

'  de 


(170) 


»  =  _/_i_  dr 

2  d£  , 


(171) 


A.  = 


+  n"  e  ~7r 

>1/1(3,)  (1+6)  ^6  , 


(172) 


Rewrite  the  cable  equations  in  their  normal  form  using  the  characteristics.  To 
find  the  characteristic  equations,  we  can  write  equation  (161)  in  terms  of  € 
using  the  constitutive  relation  (equation  (167))  multiplied  by  ^  and  subtract 
from  equation  (164)  multiplied  by  dS.  This  operation  is 

MU'* 

+  XHJt] 

+  [  ^  -  if  ‘ ^ +  (-V  H ^  n  <^s)]  *  0 

o 

(176) 


Since  there  are  only  two  independent  variables,  s  and  t  ,  the  total  differential 
of  any  dependent  variable,  for  example  V,  is 


The  variables  in  equation  (176)  are  expanded  as  follows: 
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-vy -WaHtgJt-yrastgJs-jk H  Jt-o}  (178) 
-ja  §  H'1*  -  +M-  V44-  WccsMe-iKJt*oi(m) 


In  terms  of  the  characteristic  values.  .  v/e  see  that 


-jiH)jt  -  fLjs+JV-VH-WcnW-O,  (180) 


where 
in  the  t-s 


=  ±  /nrr~1 

/><  5? 


are  the  axial  wave  characteristics 


plane;  hence, 


4S  _  ±  f±~£t' 

cAt  \J els  , 


The  characteristic  equation  becomes 


~JjT%  -o(€  +  JV~VJ<p-Wcos<t>^-±HJtt~Ot(m) 

This  equation  represents  the  motion  of  the  cable  in  the  axial  direction  (stretching). 
If  we  let  the  <x  characteristic  be  associated  with  -f-  / -L  ■dll”1  and 
represent  disturbances  traveling  down  the  cable,  we  see  that 

d.V.  -  /TIT  ds.  _  v  ■urcosS'l6  ~  J-uit  -  n 
«l«  «(6  V  (182) 

f 

Similarly,  let  xS  be  associated  with  disturbances  traveling  up  the  cable, 
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;  then. 


-t  /Z"3?  ^  -  y  -V/cos  _ t  K  2^ 

X#  //U  *<€  *08  V  <*08  W  dJS 


O 


(183) 


Again,  equation  (162)  multiplied  by  (j^itf)dt  is  subtracted  from  equation  (165) 
multiplied  by  ds  to  give 

15' °<s  +  %■ At  + v Ts  0,3  + u H M  +' WS'N *  H 0,5 
+¥ -jx  ££ ||c#  -j^Gfclt  =»  O , 
Dividing  through  by  ( 1  +€  ),  we  find  that 

-(/» 77+i)  !j?  M  ^  ^  ft  ^  +  oT«T 17  '*'(rfe)cnr 


+  rr^  ¥  5JN</)  $  clt  -O  9 


x_ 

(ii-6) 


yU(l  +  6) 

The  transverse  wave  characteristics  in  the  t-s  plane  are 


hence, 


and 


dt 


n 


X 

6 

dr' 

X 

(1  +  6) 

J6  ) 

Ijl 

6 

c it’ 

!  M 

(1  +  6) 

de  y 

e 

(i+e) 

+  1/ d<P  +~W sin  cf>  d0 

Let  the  0<  characteristic  be  associated  with  -f- 

i 

represent  disturbances  traveling  down  the  cable  in  the  y"  -  x"  plane: 


a 


6  dX 

im  St 


r~" 


0 

4 

and 
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+w -ji^as  *  °,  om) 

Let  the  j8  characteristic  be  associated  with  —  /-i-  7<— -  gf-X  and 

/>*  (7+£)  ^ 

represent  disturbances  traveling  up  the  cable  in  the  y"  -  x"  plane: 

S+hlM^-v)M  +^j-i  $$=<?.  {i85, 

Finally,  equation  (163)  multiplied  by^fc)dt  is  subtracted  from  equation  (166) 
multiplied  by  ds  to  give 

-(j+6)cos$  (!§«<$+ jx(j~§  jz  ffcte)  +  <>iw 


+  (1/  COS<{)~  ySlN<{>)  AO  ~  I  At  ~  0  t 

The  transverse  wave  characteristic  in  the  t-s  plane  is 

- - - 1 


K,.  «  ±U*£ 


'5,6 


hence 


(It  6)  d6  > 


and 


4f  =a  +  /Z  ^ 

^  JyM(l+e)  Me 


</W~ (I +e)cos (p * c(e  +  0 7  cos<$>-Vs\n<$)M6  -jjl&b  -Ot 

Let  the  <*  characteristic  be  associated  with  -f  /  /  6  gfp  and 

\/^u  (t+6)Ste 

represent  disturbances  traveling  down  the  cable  in  the  z"  -  x"  plane: 


Let  thefi  characteristic  be  associated  with 


-  /zur  xf1 

v//«  f/+6) 


and 
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represent  disturbances  traveling  up  the  cable  in  the  z"  -  x"  plane: 

M + (<  Htyaaifi  fa  g +(V 


The  characteristic  equations  are  now  summarized.  For  tensile  waves 


traveling  down  the  cable  with  velocity  •*• /_ g(2 I 


,  we  have 


\&'ir 


tensile  waves  traveling  up  the  cable  with  velocity  i?  >  we 


For  transverse  waves  in  the  y"  -  x”  plane  traveling  down  the  cable  with 


,  we  have 


velocity  -+  /-L.  -£ _  dT 

TiH-e)  5T6 


sW +u)^  15 ^ ■= 0 .  <190> 

For  transverse  waves  in  the  y"  -  x"  plane  traveling  up  the  cable  with  velocity 


we  have 


~U)$  191) 


For  transverse  waves  in  the  z"  -  x"  plane  traveling  down  the  cable  with 


velocity  -4-  /_(. _ 

VMS&M iae 


,  we  have 


smrwmaammtaet 


! 


I- 


’O' 


m 


% 


1 
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For  transverse  waves  in  the  z"  -  x"  plane  traveling  up  the  cable  with  velocity 

,  we  have 

AJ+(('+6)^*Jm(m)m' +(V»^-V«n*)^|  -^1$  =  0  ( 193) 

Equations  (188)  through  (193)  form  the  basis  for  a  numerical  solution  of 
the  three-dimensional  equations  of  motion  for  a  buoy  mooring  cable.  The 
original  set  of  six  partial  differential  equations  with  six  unknowns  has  been 
transformed  to  a  set  of  six  ordinary  differential  equations  with  six  unknowns 
with  the  restriction  that  integration  operations  must  be  carried  out  along 
characteristic  curves. 

The  solution  of  sets,  of  coupled,  nonlinear,  partial  differential  equations 
with  different  characteristic  velocities  is  usually  accomplished  by  assuming  that 
the  equations  can  be  decoupled.  At  best,  they  can  then  be  linearized  and  a 
separation  of  variables  method  can  be  used  to  obtain  an  analytical  solution.  At 
worst,  the  normal  fc^ms  can  be  integrated  numerically  along  characteristics. 

Examination  of  the  functional  form  of  the  transverse  wave  characteristics 
show  that  they  are  dependent  on  the  state  of  stress  in  the  cable.  Thus,  the  cable 
equations  can  not  be  decoupled  if  they  are  subject  to  time  or  space  varying  forces. 
Also,  the  particular  cable  system  studied  here  is  a  combined  initial  value  — 


boundary  value  system: 


jBg 


gg*yaaaaa3»aj 


gvrTOTjj.-  <iii^'Ji^»jg/».*»^iay5iBB!ESBSS 
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1.  All  the  dependent  parameters  are  known  at  time  =  0. 

2.  The  displacements  at  the  lower  bound  (anchor)  are  constants  for 
all  time  and  their  time  derivations  are  zero. 

3.  The  displacements  at  the  upper  bound  (buoy  motions)  must  be  solved 
for  simultaneously  as  are  cable  motions. 

Buoy  system  cable  motions,  as  posed  here,  can  best  be  handled  numeri¬ 
cally  with  a  modification  of  Hartree's  method,  i.  e. ,  solving  for  the  values  of 

the  dependent  variables  at  rectangular  grid  nodes  in  the  time-space  domain. 

66 

Hartree's  method,  as  described  by  Ames,  was  developed  for  a  single  hyper¬ 
bolic  equation  whose  coefficients  may  be  time  or  space  dependent.  In  this  study, 
the  basic  method  will  be  extended  for  a  set  of  coupled  hyperbolic  equations  with 
variable  coefficients. 

3.2.2  Finite- Difference  Methods 

Consider  a  rectangular  grid  in  the  time-space  plane  (figure  15).  Assume 
that  the  values  of  the  six  independent  parameters  are  known  at  the  nodal  points 
on  the  jth  time  line.  It  is  desired  to  advance  the  solution  to  the  nodal  points  on 
the  (j  +  l)th  time  line;  specifically,  for  the  ith  point  R.  If  the  characteristics 
are  known  at  point  R ,  characteristic  lines  can  be  drawn  back  from  R  to  inter¬ 
sect  the  jth  time  line.  The  six  parameter  values  at  points  B,  A,  and  D  can  be 
linearly  interpolated  to  find  the  values  at  each  intersection  point,  P^  ,  P  ,  , 

and  Q  .  Now,  the  characteristic  equations,  (188)  through  (193),  can  be  used  to 
£ * 

advance  the  solution  to  point  R  .  For  tensile  waves,  the  length  between  point  A 


4  ' 


and  point  P  in  the  grid  is 


1 


*"Si  -  k  (  ChiR  +  ChJri )  •  fc  * 

If  constant  cable  mass  and  modulus  is  assumed,  the  tensile  characteristic  is 
independent  of  the  tension;  thus. 


Ch„-  0>1P|  «  Ch,  . 

The  location  of  point  P  is  given  by  S.-Sr.  =  Ch,-k  ;  in  a  similar 

1  » I  1 

manner,  point  (associated  with  tensile  waves  propagating  up  the  cable)  is 
given  by  .  For  tensile  waves  traveling  down 

the  cable,  the  finite-difference  equation  is 

^R-ty  "  g(c*'ifi+Chlf|)(€R-€ri)  -jr(VK  +Vp)(<Pit~,Pf‘) 

+wR)[cos(£  Wfe+&,)|(e„-  ePl) 

~Jn'i (Ht+HrJ'k  =0  _ 

For  tensile  waves  traveling  up  the  cable,  the  finite-difference  equation  is 

U*~  V<a, +  2  (CK+  Chss)(e R-  £>)  -  k 

-  k  (w*  •*%,)[  cos& <  *«,))](  e*  -  %) 

=  o  , 

For  transverse  waves  in  the  y"  -  x"  plane  traveling  down  the  cable,  point  P 


(194) 


(195) 


is  located  by 


where 


Sr  ~  “  a  Chapa)  'k  > 


c^r  -J: 


j_€  dr 
M  (H-I)  Ji 


The  finite -difference  equation  is 


V*  -  V  Kck2;,+csJa-4) 

+  i  ( V,  -  4))]  C  0*  -  %) 

-jri(Q«  +  £j-k  =0  .  <196> 

For  transverse  waves  in  the  y"  -  x"  plane  traveling  up  the  cable,  point  Q  is 

t* 

located  by 

•S’r-5^  +  Cli2<?a)*k  , 

The  finite-difference  equation  is 

\-V^-k(c^ChzJ(<prfrJ 

->-*0VGJ*k«o  .  <197> 

For  transverse  waves  in  the  z"  -  x"  plane  traveling  down  the  cable  (note  that 
the  characteristics  for  the  transverse  waves  are  the  same), the  finite-difference 
equation  is 

WR  -  WPz  -  +  fl+%)  C/i2Pt]  cos 

- i(^pz)cos(i((p^(pJ-p  i(vR^)siN(M^(e R-0p2) 

~0.  (198) 

lor  transverse  waves  in  the  z"  -  x"  plane  traveling  up  the  cable,  the  finite- 
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difference  equation  is 

WR  -  W9l+h[(i*^)C^(i*^)Ch^cos(i(^*0 
+  i(Vt+v^)cos(M+^-£(v^)siN(i( 

"  Jji'  2  (  +  —0  ,  (199) 

Rewriting  the  six  finite-difference  equations,  we  see  that 

v*  -  4  (chiR*  <Vv  *v^-i(K^,psm^i]  e* 

+f  V  £  ^V^S'6'-  *  ± 

-jx-  £  f  H^+Hj-k}  =  0  ,  (200) 

14 + Kc^+ct-1Q>R  -  4  (WJ&  -  4fws  ^(cosdi^ijj  eR 

+ a  C^iu)€? + g  2  (WR  i-Wg^cos  (i  f  4))J 

~  jui’  2  (Hr^HsJ’KJ  ?  (201) 

VR  -  g-(o+6s‘)  CkR  +  ( ,+£a)CliZ(i)(j)/i  +  g^+fo^jsiN (£(V^)J  Or 

+ {~\+  £  (d+£«)Ch2j,+ O+fi/j)  C  (iZ)>2) 

j;  g(Gfi+  %)'k :]=o  I 
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¥ 


I 


R 


Vr  +  5  (d+eJChtt+O+eddifJ  <pK + n(i(^<ljj]  dR 
-g(w,-*wai)[5iN(i(i/is+^2|]0S£  -jii(GR+Ga!)-k]-  = 


(203) 


(7 


K~  ^ (0+4)  ChZfji)coS (2  (fa+tffk'j) 

-  i  M+l4)  COS (£( +&J  +  £  (V*  +Vpe)  swfety+fy))}  Or 

+  ["{wrneiwr  j  ~  /a' ^P/J'kJ  ~  0^  (2Q4) 

and 

+  [z  ({,+€,)  C^#+ (|+€<u)  Cli20i)  cos(i(fat  <paj) 

+  £(v,<+Uic)cos(i(h*/pJ-  £(V*+\)$w(£  WW<0)}  6>K 

44r;*4 v jr  ^wj-k]=a.  (205) 

In  determinant  form,  they  become 


yR  to 
uR  to 

0  tVR 

0  -vR 
0+0 
0+0 


+  0  •*•  Ai  €r  +  A, ,  +  Al6oK  +  A,7  ~  o 

+  0  4  A246r  t  Ars^r  t  Az<  4  A27  “  0 

t  0  +  0  +  A3s(pR  +  A3ff0*  •*■  A3r  =  0 
+  0+0  t avj0r  +  Aw  +  AH7  =  o 
+  H+  0  +  0  +  A56Qr+AS7  =  0 
+  WR  +  0  t  0  +  oR  t  A<7  -  0 


) 
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where  the  A  are  the  coefficients  in  the  preceding  set  of  equations,  (200)  through 
ij 

(205).  The  six  simultaneous  equations  are  solved  to  give 

W*  +  A Si&h  ~  -A  57 

-(Wh  t  A  g  Qt  -Ag ) 

( -  A6e)  0*=  A(r  -  A57 

•  A  —  ( At7~~Asr) 

,=2Z  (A „-A,)  } 

WR  +  A5<  0*  -  -A57  **«  Wjf  ~  ~  A g7  -  A5<  0* 


"Vr  +  Ajs  (j)R  —  -A3g  Or  ~  A3 7 
"•(Vr  A qy  (pR  ~  -  Aifc  0R  -  Aht  ) 
(Ajs-A^)^-  (AH6 -A36)0  +Aur-A37 

•  d)  —  (Ahs~Asc)&  +  Aj1~A37 

(A  „-Aw) 

= ■— - ) 


VR  +  ^35-  ^R  —  A3(;  A  37  ***  ~  A35r  ^r  ~ A2c  Or  "  A37 

^R  +  Ah  ~  “  A)5  -  Aj$0R  ~  A17 

"(Vr  +AZh  £r~~  Azs  fit "  At(  Or-Am) 

~~JAn~A^J^^TA^^  As)(pR  +  (A2^“  A|f)0R  +  A27“A/7 

.  £  —  (A2ff  —/A<s)  (j)R  T  ('Aag  ~Aic)0r  +  Azr-Ag 

(Aw-A„)  _  % 


~  “A){<  £g  ~  A (5  -  Ajg  0g  ~  A(7  , 


Ill 


The  numerical  procedure  within  the  time-space  plane  grid  is  as  follows: 

Step  1 

Compute  the  coefficients  of  the  six  simultaneous  equations  using  the 
values  previously  computed  along  the  two  preceding  iso-time  lines. 

Parameters  between  grid  points  are  estimated  by  linear  interpolation 
along  an  iso-time  line;  for  example, 

t'n  =tWhp,/h)(I/,-t/s) , 

where 


Accelerations  are  needed  to  compute  the  hydrodynamic  mass  terms  in 
the  loading  functions.  Accelerations  are  estimated  by  using  an  Euler 
numerical  approximation  along  an  iso-space  line,  for  example, 

Ur  =  (V«-Va)A  . 

In  order  to  begin  the  iteration,  assume  that  the  values  of  the  para¬ 
meters  at  point  R  are  equal  to  the  values  at  point  A. 

Step  2 

Solve  the  simultaneous  equations  for  the  six  parameters  at  point  R. 
Step  3 

Go  back  to  stepl  using  the  new  values  at  R  in  the  coefficients  and 
recompute  values  atR. 

Step  4 

Repeat  this  procedure  three  times.  * 

*A  more  efficient  method  would  be  to  specify  an  error  and  iterate  until 
the  computed  parameter  values  converge  within  the  eiror  band.  However, 
three  iterations  were  found  to  give  good  convergence. 
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From  the  above  scheme,  it  is  obvious  that  solutions  along  two  iso-time 
lines  are  needed  to  start  the  iteration.  Use  the  initial  conditions  along  the 
cable  at  t  =  0  to  fill  in  the  parameter  values  at  points  C  ,  E ,  and  F.  Assign 
these  same  values  to  points  along  the  cable  at  t  -  k  to  obtain  parameter 
values  at  points  B,  A  ,  and  D.  Begin  iteration  at  grid  point  (1,  2);  S  =  h, 
t  -  2k,  and  iterate  down  the  cable.  At  the  lower  end  of  the  cable,  set  S  =  h 
and  t  =  3k  and  repeat. 

Velocities  along  the  upper  boundary  are  described  by  the  motion  of  the 
cable  end.  However,  at  the  beginning  of  each  new  iso-time  line,  strain  and 
angles  must  be  computed  at  S  =  0  for  the  previous  time  increment.  Hartree's 
method  can  not  be  used  here  because  parameter  values  at  times  less  than  0  are 
not  available.  A  linear  extrapolation  along  the  previous  iso-time  line  is  used 
to  obtain  the  required  values;  for  example, 


®S  =  o  ;  6„  =  6, 

The  same  problem  exists  along  the  lower  boundary;  only  here  the  velocities 
are  set  to  0.  Once  again  the  values  of  strain  and  angles  are  determined  by 
linear  extrapolation  along  the  previous  iso-time  line. 

This  method  fails  to  converge  if  points  P  or  P  fall  outside  the  space 
interval  A-D.  Since  the  locations  of  these  points  are  determined  by  the 
characteristics,  the  minimum  relative  size  of  the  space-to-time  increment  is 
equal  to  the  value  of  the  greatest  characteristic.  The  tensile  wave  character¬ 
istic  is  always  greater  than  the  transverse  wave  characteristic;  thus, 


« 


1 
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The  maximum  space  increment  is  determined  by  the  size  of  the  smallest 
wavelength  in  the  cable.  The  transverse  wave  characteristic  will  determine 
the  smallest  wavelength.  If  the  transverse  wave  characteristic  is  known  to  be 
250  ft/sec  and  the  cable  is  being  excited  by  a  sinusoid  with  a  2-sec  period,  the 
wavelength  will  be  500  ft  in  length.  At  least  ten  points  are  needed  to  describe 
a  sinusoid;  thus,  the  maximum  size  of  the  space  increment  is  50  ft.  The  time 
increment  for  a  3,000-ft/scc  tensile  wave  characteristic  is 


h 


50  ft. 


=  0.01  a  7  sec 


ft**.  Lb.  3000  Ft./sec 

Note  that  the  characteristics  must  always  be  real,  finite  values.  If  the 

characteristics  are  zero  or  imaginary,  the  equations  become  ultrahyperbolic 
with  multiple  solutions  for  a  given  set  of  initial  conditions.  In  this  system, 
this  is  possible  if  the  tensions  are  less  than  zero.  Transverse  wave  character¬ 
istics  are  then  imaginary  and  the  problem  is  indeterminate.  This  study  is 
further  constrained  by  the  requirement  that  the  tension  must  be  greater  than 
zero  at  all  times  over  the  whole  length  of  cable. 

3.2.3  Cable  Lo^  ang  Functions 

The  loading  functions  H,  G,  and  I  (figure  16)  in  the  cable  equations  are 
composed  of  the  weight  components  in  the  double-primed  coordinate  system, 
the  steady-state  and  dynamic  drag  components  (viscous  forces),  and  hydro- 
dynamic  inertia  components  (acceleration  proportional  forces). 

The  weight  components  are  found  by  simply  transforming  the  weight  per 
unit  length  from  the  inertial  coordinate  system  to  the  cable  coordinates.  The 
weight  components  are  ~Wc  =  ,  or 


Figure  16.  Cable  Loading  Functions 


# 


(206) 


'  wc  cos  <j>  cose 
“ wc  siN(j>  cose 
I"  We  SIN  6 


Drag  forces  are  assumed  to  follow  a  velocity- squared  drag  law  of  the  form: 


D-  ^4AV2. 

On  a  unit  length  basis,  they  become 

J)^±J>CpdVM) 

where  d  is  the  cable  diameter.  In  order  to  maintain  a  sign  convention, 
equation  (208)  is  rewritten  as 


ivl. 
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Casarella  and  Parsons  have  reviewed  the  state  of  the  art  for  analysis 
of  hydrodynamic  forces  on  cable  systems .  Two  approaches  are  described:  the 
use  of  loading  functions*  to  compute  the  normal  and  tangential  drag  force  compo¬ 
nents,  and  the  direct  computation  of  normal  and  tangential  drag  using  normal 

11  12 

and  tangential  drag  coefficients.  The  latter  method,  used  by  Wilson  *  for 
mooring  problems,  is  employed  in  this  study  because  it  lends  itself  to  three- 
dimensional  cable  problems  and  only  requires  a  simple  transformation  of  the 
relative  velocity  components  to  cable  coordinates.  The  first  method  requires 

*Loading  functions  are  defined  as  the  ratio  of  the  hydrodynamic  force 
component  of  interest  to  the  drag  force  when  the  cable  is  oriented  normal  to  the 
flow. 


multiple  transformations  to  define  the  normal  and  tangential  forces  and  to 
resolve  them  into  three  force  components. 


If  the  current  velocity  components  are  given  by 


(no  vertical  currents)  and  the  velocity  components  of  the  cable  element  are 


then  the  velocity  components  of  the  water  relative  to  the  cable  are 

VJ  \o~k' 

V*  = 

w«  W3-z 

Transforming  to  cable  coordinates  gives 

I/r  -XC0S<f>C05e+(Vs-?)SlN<p  +  (Ws-£)  cos  <f>  SING 
VR  -  X  SIN(j)  cose  +(VS-?)  COS<}>  -  ('Ws-i)siN<)>SlN$ 

W*J  L  *  5/n©  +  (ws-z)  cose 
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The  drag  force  components  are  now  written: 


d„.  =  ij>c„dv:m 


(210A) 


Dy.« 


(210B) 


v,„  =  is  cMdw;\w;  i. 


(210C) 


Hoerner,^8  Wilson,  ^  and  others  show  plots  of  normal  and  tangential 
drag  coefficients  versus  Reynolds  number  for  stranded  cables.  Wilson's  data 
were  used  in  the  G.  S.  A.  program  "CURFIT, "  and  the  following  expressions  for 
normal  and  tangential  drag  coefficients  were  developed: 


CDN  =  1.3Z027  +  10,Gd$Z/Rz 


0<  Re  <  100 


-  l.H 


IOO<Re<5*lo% 


CDr  =  O.G05H6-(Rz 


-O.H75t\ 


0<  Rz<5xlos 


(211B) 

where  R.-TC4  ,»  is  the  kinematic  viscosity.  These  expressions  are 

5 

valid  only  in  the  subcritical  Reynolds  number  region  (Re  5  X  10  ).  However, 
since  this  Reynolds  number  would  apply  to  a  2-in.-<diameter  cable  moving  at 
35  knots,  it  is  felt  that  the  range  of  Reynolds  numbers  covered  is  adequate. 

Hydrodynamic  inertia  forces  are  defined  normal  and  tangential  to  the 
cable.  Transforming  the  accelerations  to  cable  coordinates  and  assuming 
steady  ocean  currents,  we  find  that  the  accelerations  become 
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j(//  X  C0S<p  C0s6  +  ¥  SIN  *  *+  H  cos*  57A/  0 
y"  ta  - X  SIN  Ip  cose  -»•  ycostp  -  ji  5/a/*  sva/<9 
z"  -x  SlMd  +  2  COS  6 

J  W  9 

The  hydrodynamic  inertia  force  vector  is 

>U  K,-  o  O  irr 

0  nv  o  •  /" 

[  0  O  hv-J  [  z"J . 

The  hydrodynamic  mass  dyadic  is  a  diagonal  matrix  due  to  the  axisymmetry  of 

the  cable.  Also,  for  a  cylindrical  object  of  infinite  length,  the  tangential  hydro- 

dynamic  mass  (  Whxx«)  is  zero.  For  a  smooth,  constant  diameter  cable,  the 

normal  hydrodynamic  masses  <  mKr-,  Wh„  )  are  equal. 

45  72 

Lamb,  Basset,  and  others  have  used  potential  flow  theory  to  com¬ 
pute  the  normal  hydrodynamic  mass  of  a  circular  cylinder: 


mhM  =  tt/1  .  (a*) 
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Miller  and  Miller  and  Hagist  have  investigated  the  frequency  dependence  of 
hydrodynamic  mass  for  various  bodies.  Their  data  show  a  linear  decrease  in 

hydrodynamic  mass  with  increasing  frequency  in  the  Stokes  number  region 

5  74 

0  *<  St.  <  3  X  10  .  By  using  Miller's  data  for  a  5:1  cylinder,  we  can  see 

__  0 

that  the  slope  is  -1.62X10  .  Equation  (212)  is  modified  to  include  the 

frequency  dependence  as  follows: 

=  (l0~ (L 62 xW'S]‘  17 ) ^/H  5  0<^<3\m) • 

The  components  of  the  hydrodynamic  inertia  force  are 


FT  -  0 


(214A) 
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(2)4  B) 


fL.  =  mh  •  z" . 


(214C) 


The  loading  functions  are  summarized  by  using  equations  (206),  (210),  and 


(214): 


H  —  Wc  COS  (/)  COS  6  -t-  (0.605H6>(R<Zv„  Vr  /l/fi  I  (215A) 


-Wc  SIN <j>  COSO +  ^  dVft  IVj 

+  (lo-(i&Aio-‘)-qf)'Trj>£'  y"  > 


(215B) 


=  -WcSlNB  +  C^,,  ^  dW/IWfil 

+  ( lo  -  (l  62  *10*0'  ^Z)tTJ>  fz" 


(215C) 


Re  =  ^ 

r'eU"  -y 


Rev, 


R  »  — 


=  Y^< 

'V 

_  W'W 


where 


CM  =  1.3&0Z7  +  i0.616Z/fov. 
=  1.H 


0<  fiev.<  100 
I00<  fay  <  5  *los 


CL  ..  =  152027+  iO.  6K£/ Rtw* 

uriVt 

a  1.*/ 


0<Riwh<  l 00 
)00<  <5*  10s 


3.2.4  Lumped-Mass  Model 

The  finite-difference  analysis  for  cable  dynamics  can  become  very 
expensive  in  computer  time  and  can  require  large  amounts  of  computer 
"storage"  because  of  the  large  arrays.  A  lumped- mass  analysis  can  offer 
significant  savings  in  computational  time  at  the  expense  of  simulation  accuracy. 
In  general,  the  lumped-mass  analysis  will  truncate  the  high-frequency  response 
of  the  system.  However,  for  many  engineering  applications,  the  high-frequency, 

V 

low-amplitude  response  is  not  of  interest,  and  the  cable  can  be  represented  as 

a  small  number  of  lumped  masses. 

% 

Assume  that  a  uniform  cable  of  length  L  can  be  broken  up  into  N 
segments  of  equal  length  (figure  17).  The  length  of  each  segment  is 

aL  *  L/n 

t 

Assume  that  the  properties  (weight,  mass,  hydrodynamic  forces,  etc.)  of  the 
cable  can  be  concentrated  at  points  2,  3,  . . . ,  N,  which  are  located  A  L, 


dKbWKMMM a 
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2AL. (N-l)AL  from  the  upper  end  of  the  cable.  All  the  forces  acting 

AL 

on  the  cable  span  from  (N  -  1)  4  L  _+  *—  are  concentrated  at  the  nth  mass 

point.  Cristecu's  cable  equations,  (135A),  (135B),  and  (135C),  are  written  for 
each  mass  point: 


yU&UL  0  *XnAL  +  (f.-1-Z,-f)  , 

/a(s.)aL  0  «  Ym  aL+  (%3 -  ,  <216B> 

and 

/j.(S,)aL  <0  =  ^nAL  +  (T„4-Tn.,4)  .  <2i6C) 

This  is  equivalent  to  modeling  the  cable  with  a  system  of  spring-mass  elements 
as  shown  in  figure  18.  In  order  to  compute  the  forces  acting  on  each  mass 
element,  the  cable  angles  (  0  and  9  )  for  each  cable  element  must  be  defined. 
From  the  geometry  between  the  nth  and  (N  +  1  )th  mass  element,  we  see  that 


e„  =  sin 


r1/  ~Zn) 


J~CXt\+i~Xn)  t 


and 


(217) 


<j K  =  SIN 


~‘f  (y.„-y„) 


(218) 


Tensions  between  the  lumped  masses  are  computed  by  using  the  elastic 


properties  of  the  cable  and  the  deformation  of  the  cable.  If  the  effective  cable 


modulus  is  E  and  the  deformation  of  the  cable  of  length  A  L  is  S  the 
c 

spring  constant  along  the  cable  is 


-n-  $r<r_  _  /  Be 

8  "  *f  AL  - 


(219) 


The  stretched  length  between  the  nth  and  (N  +  l)th  mass  element  is 

If  the  difference  between  the  stretched  length  of  the  cable  and  the  unstretched 
length  (AL)  is  less  than  zero,  the  tension  is  zero  since  the  cable  can  not  sup¬ 
port  compression.  (This  is  analogous  to  the  ultrahyperbolic  equations  that 
occur  in  the  finite-difference  analysis  if  tensions  go  to  zero. )  Otherwise,  the 
tension  in  the  nth  cable  segment  is  defined  as 


T^„  -  (j (xn+1  -xnf -t(yriH~Y„t+  (z«*r 1 zi>)* '  -  a  l)  # 

Transforming  to  inertial  coordinates,  the  tension  becomes 


(220) 


T„x„  cos  <fi„  cos  e„ 
sin  <t>„ 

T  cost  SIN  e„ 

0 

(221) 

These  tension  components  are  used  in  equations  (216)  to  compute  the  tension 
difference  across  the  mass  element. 

As  before,  the  forces  X„aL,Y„aL  and  aL  acting  on 


each  mass  element  consist  of  weight,  viscous  drag,  and  hydrodynamic  inertia 
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forces.  The  weight  force  vector  per  unit  length  is 

where  Wc  is  the  weight  per  unit  length  in  water  of  the  cable.  Again,  if  the 
current  velocity  components  at  the  nth  element  are 


fo 


and  the  velocity  components  of  the  nth  element  are 


% 

L^J  ) 


then  the  relative  velocity  components  are 


7? 

3 

_ l 

i - 

£ 

i 

o 

■ 

11 

*c 

~  i 

5<* 

I 

C 

'  (222) 
If  the  mean  cable  angles  at  tue  nth  mass  element  are  computed,  we  find 


<t>n  =  ±(<fc,  +  0 


(223A) 


and 


126 


©„  =  f  ( e„  +  e„  J  . 


(223  B) 


The  relative  velocity  components  are  transformed  to  cable  coordinates  as 
follows : 


I cos  ^  cos  e„  +  (\kB-  y„)  SIN  +(WJn  -  H„)  cos  <j>„  .sin  0„  , 


V,"  =  X„  sin  ^  cos s„  +  (\&B-  X„)  cos  <J>„  —  (V&„' i»)  sin  <j>„  Sin  e„  , 

and 

W„"  =  An  SIN  0„  +  (w5n  -  H,,)  COS  e„  . 

The  drag  forces  per  unit  length  are 

i,;  „  j  /.  c„  d  K I UC  I  , 

K -  is  c„  J «. lull , 

and 

where  the  drag  coefficients  are  computed  from  equations  (211A)  and  (211B). 

The  accelerations  of  the  nth  mass  element  in  cable  coordinates  are 


K  =  x„  cos  <#>»  cos  e„  +  x,  sin  +  z„  cosfy  sin  %  j 
y"  =  -X„  SIN  $„  COS0n  +  X,  CCS&  -  S'H  SIN  §„  j 

iC  = -X,  SIN  e„  +  2„  C os§„  . 
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The  hydrodynamic  inertia  force  components  per  unit  length  are 


r-  “ 

fL  =  o 


nr  ^  **  ■/ 

'Hriy  ) 


Pr  11  •’ » 

►"»*„  =  -nVx»  , 


(224A) 


(224B) 


(224C) 


where  JT).  is  computed  from  equation  (213).  The  drag  and  hydrodynamic 
inertia  forces  per  unit  length  are  summed  and  transformed  back  to  inertial 
coordinates  and  added  to  the  weight  force;  thus, 


X]  M  br+o  " 

r„  =  o  -  -4  •  d;  +  c* 

XJ  M  K-  X,„ 


(225) 


All  terms  have  been  defined  in  the  equations  of  motion  for  the  cable 
elements.  The  three  equations,  (21GA),  (216B),  and  (216C),  must  be  integrated 
simultaneously  for  each  mass  element.  Thus,  if  the  cable  is  broken  up  into  ten 
segments,  there  will  be  nine  mass  elements,  each  having  three  degrees  of 
freedom.  Therefore,  there  will  be  9  x  3  x  2  =  54  simultaneous  first-order 
equations.  Note  that  the  upper  half  of  the  first  cable  length  and  the  lower  half 
of  the  last  cable  length  are  not  included  in  these  equations  of  motion.  These 
cable  segments  are  assumed  to  be  moving  with  the  upper  and  lower  boundaries 


i 
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and  their  properties  should  be  lumped  in  with  the  properties  of  the  buoy  and 
anchor.  The  computer  program  developed  for  the  lumped- mass  simulation  of 
cable  dynamics  will  be  discussed  in  a  later  section. 

3.3  Steady-State  Buoy  System  Configurations 

As  an  introduction  to  the  coupling  of  the  buoy  and  cable  equations  of 
motion  and  the  resulting  computer  programs,  the  steady-state  buoy  system 
configuration  (zeroth-order  case  of  buoy  system  dynamics)  will  be  investigated. 

The  specification  of  the  proper  mooring  line  length  for  a  moored  buoy  system 
is  critical  in  the  design  of  the  system  in  order  to  avoid  tow-under  of  the  buoy, 
minimize  the  "watch  circle"  of  the  buoy,  reduce  steady-state  tensions  in  the 
moor,  etc.  This  analysis  offers  a  method  to  select  the  proper  mooring  line 
length  for  a  given  buoy,  water  depth,  and  current. 

If  the  time  dependent  terms  in  the  cable  equations  (equations  (161)  through 
(167))  are  allowed  to  go  to  zero,  we  find  that 

+  H  -  0  ,  (326) 

T  G  =  0  >  (227) 

T  COS  (f>  ||  +  I  «  0  ^  (228) 

and 


(229) 
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Lotting  cable  velocities  go  to  zero  in  the  loading  functions,  we  find  that 

H=wc  COS  if)  cose  +  (o.60Si6-(f^^Srj^  dl/„  jUfl  J 

G  =  ~wc  siH<t>cose  +  CWr,  ?  ol  VijV'"}  , 

and 

I  =  -Wfc  SIN€>  +  C„Nv„  £  dWfi  l  W„"l  y 

where 


SlN<p+ ys  COS  (j*  SIN  6 

VR"  »  X  ZOS^-Ws  SIN <p  SINS 
W«  =  W5  COS0 

ISZ0Z7  +  iO.CHZ/Rty* 
CPN . =  1.3Z0Z7  +  1ft  6%z/tew, 

n 

»!// 


0<  Rev»  <  I OO 
/0O^  RZy'<5*tvs 

0<  Rev*<  j0O 
100  <  Re*,”  <5 *I0S 


Tf"d 

v 


V^d 

v 


Expanding  equations  (226),  (227),  and  (228)  to  include  the  loading  functions  and 
taking  total  derivathes,  we  find  that  the  steady-state  cable  equations  become 
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—  =  -wc  cos  <f>  cose  -  (o,605°t6  (Kc^75  ))^T  d  V#  )Ur  I  >  (230) 


j|-  f(t*  *»♦«*•  -  C„Uv„ 4.  ^''|V/f)  , 

3i  =  fa&(w‘sme~  WW), 


(231) 


(232) 


(233) 


Before  integrating  these  equations,  the  upper  boundary  conditions  (at  the 
buoy)  must  be  determined.  The  tension  at  the  buoy  is 


V<  ^,SP  “  ^8  T *  ^r«S2  1  >  (2J 

where 

Disp  is  the  buoy  displacement 
Drag  is  the  buoy  drag  force  due  to  surface  currents 
is  the  buoy  weight. 

For  convenience,  let  the  z  axis  of  the  inertial  coordinate  system  be  aligned 
with  the  surface  current.  Then,  the  initial  buoy  moor  angle  is 


a  - 


(P»5p  J  ) 


and  the  initial  moor  transverse  angle  is 


t-0 


(236) 


Unfortunately,  the  displacement  and  drag  of  the  buoy  arc  functions  of  the  cable 

tension  at  the  buoy.  Thus,  a  trial  and  error  solution  must  be  used.  For  a 

given  current  profile  ^  ys  =  %(*)  ;  wt  =£(*)  ,  the  depth  of  water  is 

the  controlling  parameter.  When  the  vertical  projection  of  the  cable  is  equal 

to  the  depth  of  water,  the  correct  solution  has  been  obtained.  The  differential 

equations  defining  the  cable  shape  (equations  (230)  to  (233))  are  controlled  by 

the  unstressed  cable  length  Sq  .  Each  iteration  must  take  into  account  the 

stretch  in  the  segment  before  the  x  component  of  segment  is  computed.  Thus, 

if  the  integration  step  size  is  dS  and  the  tension  and  angles  for  that  segment 

o 

are  T ,  Q  ,  and  <f>  ,  the  stretched  length  is 


dS  =  (l+£)420  0 

From  the  cable  properties,  we  see  that 

c  _  <T  IT 

fc  -  £c  IT  olsEo  . 

Rewriting  equation  (237),  we  see  that  the  stretched  length  increment  is 


dS  ~(i  +  fc )  dSo  # 

The  components  of  the  cable  incremental  length  dS  are 


(238) 


dx  =  d$  COS<j>  COS  $  J 

dy  -  dS  Sin  <p  ) 
and 

dz  =  dS  cos  sin  &  % 


Since  both  the  displacement  and  drag  of  the  buoy  are  functions  of  the  draft  of 
the  buoy  (equations  (52)  and  (95)  developed  for  an  oblate  spheroid  —  similar 
equations  can  be  developed  for  buoys  of  different  shape),  the  draft  H  is  incre¬ 
mented  upward  from  its  "free-floating"  value,  i.e. ,  as  if  the  buoy  were  floating 
on  a  calm  surface  with  no  mooring.  The  configuration  of  the  cable  is  then 
computed  using  a  fourth-order,  Iiunge-Kutta  numerical  integration  algorithm 
for  the  given  current  structure,  which  may  vary  in  magnitude  and  direction  as 
a  function  of  depth.  The  components  of  the  incremental  length  dS  are  com¬ 
puted  and  summed.  The  x  component  of  the  end  of  the  cable  is  tested  logically 
to  see  if  it  falls  within  a  specified  error  band  about  the  water  depth.  If  the 
computed  vertical  projection  is  less  than  the  water  depth,  the  buoy  draft  is 
incremented  upward  and  the  process  is  repeated.  If  the  computed  vertical 
projection  falls  within  the  error  band,  the  solution  can  be  accepted  or  the 
width  of  the  error  band  can  bo  reduced  and  the  process  repeated  until  the 
solution  achieves  the  desired  accuracy.  This  process  is  shown  schematically 
in  figure  19. 

This  method  has  been  programmed  in  FORTRAN  IV  and  is  shown  in 
appendix  B.  Subroutines  arc  included  for  the  displacement  and  drag  of  oblate 
spheroids,  spherical  buoys,  cylindrical  buoys  (cans  or  spars),  torroidal  buovs 
("donut"  buoys),  and  discus  buoys  ("monster"  buoys).  The  user  must  input 
the  buoy  characteristics  (weight,  dimensions,  etc.)  the  cable  characteristics 
(diameter,  length,  cable  modulus,  and  weight  per  foot  in  water  -  negative  if 
buoyant),  and  the  environmental  characteristics  (water  depth,  current  struc¬ 
ture,  etc.).  The  output  includes  unstretched  length,  stretched  length,  tensions. 
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the  cable  angles,  and  the  x,  y,  z  coordinates  of  each  integration  step.  The 
step  size  of  the  buoy  draft  is  incremented  upward  as  some  percentage  of  the 
total  buoy  draft  (as  if  it  were  fully  submerged).  Downward  buoy  draft  incre¬ 
ments  are  taken  as  half  the  prior  upward  increments  to  avoid  a  lock-step 
situation  where  the  solution  swings  between  two  values  and  never  converges. 
Subsequent  draft  increments  are  halved  to  avoid  the  same  situation.  If  the 
initial  error  band  specified  is  too  narrow  and  the  solution  overshoots  the  error 
band,  the  solution  may  converge  too  slowly.  A  limiter  is  built  in  to  stop  the 
program  after  15  configurations  are  computed.  The  user  should  then  open  up 
the  initial  error  band.  If  the  buoy  draft  is  increased  to  its  maximum  value, 
the  program  prints  out  a  statement  that  the  buoy  sinks  and  the  computation  stops. 
This  indicates  that  the  buoy  tows  under  either  because  it  has  insufficient  excess 
buoyancy  or  the  system  drag  is  too  high. 

.  If  the  mooring  line  is  made  up  of  more  than  one  type  of  cable,  rope,  or 
chain,  logical  "IF"  statements  are  used  to  change  the  cable  properties  at  the 
proper  cable  lengths.  The  writer  has  found  that  at  the  transition  of  cable  to 
chain,  the  integration  step  size  must  be  reduced  to  ensure  numerical  stability. 

Also,  if  the  mooring  line  has  objects  (instruments,  buoys,  weights,  etc.) 
attached  to  it,  the  change  in  tension  and  angles  across  the  discontinuity  must  be 
computed  from  the  free  body  of  the  object  (figure  20).  If  the  cable  parameters 
just  above  the  discontinuity  are  T#  >  ( PH  and  0H  ,  the  tension  compo¬ 

nents  at  that  point  are 

~  %  COS  </>H  COS  Oh 


) 


(239A) 
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Th„  =  T„  SIN  (pH  , 


(239B) 


—  Th  COS  <pH  SIN  &H  ' 


(239C) 


By  summing  the  forces  on  the  free  body,  we  can  determine  the  tension  compo¬ 
nents  below  the  discontinuity.  The  tensio-i  components  are 


where 


\  -  TNy  +  D0/  , 
\  -  THe-D„t  , 


K  is  the  in-water  weight  of  the  object 

is  the  drug  of  the  object  in  the  minus  y  direction 
Dtf-  is  the  drag  of  the  object  in  the  minus  z  direction  # 
The  cable  tension  just  below  the  discontinuity  is 


^  -yxSiyc  } 


and  the  cable  angles  become 


eL=  TAN-‘(rLyrJ 
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■V 

and 

A  =  tan~1(  %yjrJ7^) 

TL  N  7  '  *  7  •  (242) 

The  integration  down  the  cable  can  resume  using  the  new  tension  and  angles. 
The  process  can  be  repeated  for  other  objects  along  the  mooring  line. 

Computer  programs  developed  for  buoy  system  statics  and  dynamics  are 
shown  and  discussed  in  appendix  B.  The  experimental  data  taken  and  the 
validation  of  the  analytical  models  will  now  be  discussed. 


IV.  EXPERIMENTAL  MEASUREMENTS 
AND  COMPARISON  WITH  MODELS 

4. 1  Steady-State  Buoy  System 

Configurations 

The  steady-state  analysis  of  the  buoy  system  configuration  is  important 
since  it  serves  as  the  set  of  initial  conditions  for  the  dynamic  analysis.  Also, 
the  steady-state  analysis  is  the  first  step  for  the  buoy  system  designer  to  ensure 
that  his  system  will  not  tow-under  or  have  other  undesirable  static  characteris¬ 
tics.  The  analysis  for  both  shallow  and  deep  water  oceanographic  buoy  systems 
are  compared  with  data  taken  at  sea  to  validate  the  steady-state  analysis. 

4. 1. 1  Torroid  and  Current  Meter 

Array  at  Station  BRAVO 

On  22  August  1967,  the  writer  installed  a  buoy-supported,  current  meter 
array  at  station  BRAVO  (41°  51. 15‘N,  71°  46. 50'W)  at  the  Block  Island  - 
Fishers  Island  (BIFI)  Oceanographic  and  Acoustic  Range  in  Block  Island  Sound. 
This  array  was  recovered  on  19  September  1967,  and  the  current  meter  data 
were  analyzed.  The  components  of  the  array  are  shown  in  figure  21.  The  three 
Braincon  type  31G  current  meters  were  suspended  below  the  buoy  at  cable  lengths 
of  15,  60,  and  105  ft.  The  buoy  and  array  were  moored  with  a  70-lb  Danforth 
anchor  and  a  100-lb  lead  weight  to  keep  the  line  pull  horizontal.  The  buoy  used 
in  this  experiment  is  shown  in  figure  22. 
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Current  Meter  Array  at  Station  BRAVO 
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Twenty-eight  samples  of  the  current  meter  data  were  selected  and  are 

shown  in  table  1.  Data  from  the  bottom  current  meter  were  highly  variable  and 

erratic  and  were  not  regarded  as  reliable.  The  data  shown  are  10-min  time 

averages  due  to  the  photographic  method  used  to  record  data.  The  instruction 

75 

manual  for  the  type  316  current  meter  gives  the  following  ranges  and 


accuracies: 

/ 

Range 

Accuracy 

Current  speed 

0.5-5  knots 

+  0. 15  knot 

Current  direction 

0  -  360  deg 

^  5  deg 

Current  meter  tilt 

0  40  deg 

,  +  1  deg 

Timing  mechanism 

5  months 

+  10  scc/day . 

The  weight  in  sea  water  of  the  current  meter  is  given  in  the  instruction  manual 

as  67  lb.  The  cylindrical  body  of  the  current  meter  is  8k  in.  in  diameter  by 

38|  in.  long  with  a  36-in.  vertical  vane  attached  to  al'gn  it  with  the  flow. 

76 

Sunblad  gives  a  normal  drag  coefficient  of  0.59  based  on  frontal  area  derived 
from  tow  tank  data  for  the  type  316  current  meter. 

The  torroidal  buoy  has  an  8-ft  outer  diameter  and  a  3- ft  hole  through  the 
center.  The  buoy  weighed  1200  lb  in  air.  The  mooring  line  was  5/8- in. -diam¬ 
eter  polypropylene  rope  having  0. 02  Ib/ft  buoyancy.  The  3/8-in.-diameter 

wire  rope  at  the  anchor  weighed  0.2  lb/ft.  The  cable  modulus  of  elasticity  for 

5  2 

the  polypropylene  rope  was  taken  as  1.67  x  10’  lb/in.  (50°  F)  and  for  the  steel 

cable  as  12.0  x  10^  lb/in.“  . 

77  78 

Williams  and  Nalwalk  el  al.  have  made  current  measurements  at 
station  BRAVO  and  have  found  a  two-layer  current  structure.  On  an  ebb  tide 
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<0 


(current  setting  to  the  easi)  the  bottom  layer  appears  to  set  to  the  northeast. 

On  a  flood  tide  (current  setting  to  the  west)  the  bottom  layer  sets  to  the  north¬ 
west.  Williams  and  Nalwalk  never  observed  bottom  currents  greater  than 
1.5  knots,  which  make  the  data  from  the  bottom  current  meter  questionable. 

Also,  the  high  tilt  angles  of  the  bottom  current  meter  preclude  proper  response 
of  the  Savonius  rotor  to  the  ambient  currents.  In  the  following  studv,  which  was 
made  to  compare  predicted  current  meier  tilt  angles  with  observed  tilt  angles, 
the  current  is  modeled  as  an  upper  layer  having  a  thickness  of  70  ft  and  a 
bottom  layer  with  a  thickness  of  50  ft.  The  current  speed  and  direction  of  the 
upper  layer  is  approximated  by  the  mean  of  the  speeds  and  directions  from  the 
two  upper  current  meters.  The  speed  and  direction  of  the  lower  layer  is 
assumed  to  be  equal  to  the  speed  and  direction  from  the  bottom  current  meter. 
Bottom  curi’ents  greater  than  1.5  knots  were  set  equal  to  the  value  of  the  upper 
layer  current. 

The  computer  program  for  steady-state  buoy  configurations  shown  in 
appendix  B  was  modified  to  include  the  effects  of  the  current  meters  and  lead 
weight  by  solving  the  free  body  at  each  object  (equations (239A)  through  (242)).  The 
torroidal  buoy  subroutine  was  used  and  cable  configurations  were  integrated 
with  a  fourth- order,  Runge-Kutta  algorithm  using  a  1— ft  step  size.  Current 
meter  tilt  angles  were  computed  by  balancing  moments  on  each  current  meter. 
The  program  was  run  on  the  CiSA-3fiO  time  sharing  computer  and  results  for  the 
first  ten  eases  are  shown  in  table  2  and  in  figure  23.  Average  tilt  angle  errors 
for  the  ten  cases  are  as  follows: 
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Figure  23.  Current  Meter  Tilt  Angles  Computed  by  Integrating 
Down  the  Cable  Compared  With  Observed  Data 


Upper  current  meter 


+  0. 22  deg 


Middle  current  meter  +0.2  deg 

Lower  current  meter*  +  1.06  deg 

Overall  average  error  +0.49  deg. 

With  the  exception  of  the  lower  current  meter,  these  errors  fall  within  the 
+  1  deg  accuracy  of  the  tilt  indicator  indicating  good  agreement  of  the  computer 
model  with  observed  data.  The  computer  study  indicated  that  the  100-lb  weight 
was  never  picked  up  off  the  bottom  by  the  strongest  currents  and  that  the  weight 
and  drag-  of  the  current  meters  control  the  buoy  system  configuration  to  a 
greater  degree  than  the  weight  and  drag  of  the  polypropylene  rope.  With  this  in 
mind,  a  simple  statics  model  of  the  buoy  system  was  developed;  it  was  assumed 
that  the  rope  was  not  deflected  between  current  meters  and  that  half  of  the  drag 
force  acting  on  each  rope  span  could  be  assumed  to  be  concentrated  at  the  end 
of  the  current  meter  to  which  it  was  attached.  Since  the  tensions  in  the  inte¬ 
grated  cable  configurations  were  observed  to  be  very  small,  the  buoy  draft  was 
computed  for  the  buoy  weight  and  the  weight  of  the  current  meters  only.  Thus, 
the  vertical  force  component  is  equal  to  the  weight  of  the  current  meters.  The 
simplified  model  is  shown  on  figure  24.  Lateral  deflections  computed  by  the 
integrated  configurations  were  very  small;  thus,  the  simplified  model  was 
restricted  to  two  dimensions.  Drag  forces  on  the  rope  spans  and  on  the  current 
meters  were  computed  as  if  the  currents  were  acting  normal  to  the  rope  or 
current  meter.  The  small  buoyancy  and  the  stretch  of  the  rope  were  neglected. 

*The  errors  of  the  lowest  current  meter  were  not  computed  for  data  where 
the  observed  tilt  angle  was  40  deg. 


Figure  24.  A  Simple  Statics  Model  of  the  Current  Meter  Array 


Buoy  system  configurations  for  the  case  at  current  meter  reference  time 
of  0.333  hr  are  shown  in  figure  25  and  are  quite  similar.  Computed  current 
meter  tilt  angles  were  compared  with  the  observed  tilt  angles.  The  angles 
were  computed  by  using  a  uniform  current  equal  to  the  mean  value  of  the  current 
speeds  from  the  upper  two  current  meters.  The  angles  are  shown  in  table  3  and 
on  figure  26.  The  average  errors  for  the  first  ten  data  sets  are  as  follows: 

Upper  current  meter  -  2. 15  deg 

Middle  current  meter  -  0.26  deg 

Lower  current  meter  +7.64  deg 

Overall  average  error  +  1. 74  deg. 

These  errors  indicate  that  the  simple  statics  model  is  about  three  times  less 
accurate  than  integration  down  the  cable  but  may  be  adequate  for  engineering 
applications.  The  average  errors  for  28  data  sets  are  as  follows: 

Upper  current  meter  -  3.77  deg 

Middle  current  meter  -  2. 69  deg 

Lower  current  meter  16.66  deg 

Overall  average  error  3. 4  deg. 

The  preceding  study  indicates  that  the  steady-state  buoy  system  con¬ 
figuration  model  can  predict  current  meter  inclination  angles  to  within  \  deg  on 
the  average.  No  tension  data  were  recorded;  thus,  the  steady-state  tension 
errors  were  not  computed.  The  shallow  water  buoy  system  described  above  is 
most  heavily  influenced  by  the  weight  and  drag  of  the  current  meters.  In  most 
deep  water  buoy  systems,  the  weight  and  drag  of  the  mooring  cables  are  the 


predominant  forces. 


Comparison  of  Computed  Buoy  System  Configurations 
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TABLE  3.  COMPARISON  OF  OBSERVED  CURRENT  METER  TILT  ANGLES 
WITH  DATA  COMPUTED  BY  ASSUMING  NO  CABLE  CURVATURE 


Time 

Speed  (knots) 

Observed  Tilt  Angle 
(deg  +  1  deg) 

Computed  Tilt  Angle  (deg) 

Tilt  Angle  Error  (deg) 

(hr) 

LJfJ 

60 

15 

[6° 

j  105 

15 

60 

105 

15 

60 

105 

.33 

1.177 

1.27 

10 

22 

40 

6.90 

20.77 

62.53 

-3.10 

-1.12 

- 

.66 

1.205 

1.16 

10 

20 

40 

6.45 

19.50 

60.90 

-3.55 

-  .50 

- 

1.00 

1.079 

1.08 

7 

18 

40 

5.38 

16.45 

56.27 

-1.62 

-1.55 

- 

1.33 

1.099 

1.07 

6 

14 

40 

5.43 

16.59 

56.51 

-  .57 

2.59 

- 

1.66 

.945 

.928 

6 

11 

18 

4.06 

12.53 

48.42 

-1.94 

1.53 

30.42 

2. 00 

.836 

.713 

5 

7 

40 

2.78 

8.64 

37.62 

-2.22 

1.64 

- 

2.33 

.717 

.578 

4 

7 

40 

1.94 

6.06 

28.31 

-2.06 

-.94 

- 

2. 66 

.538 

.497 

4 

5 

40 

1.24 

3.88 

18.99 

-2.76 

-1.12 

- 

3.00 

.509 

.341 

3 

3 

10 

.84 

2.62 

13.07 

-2.16 

-  .38 

3.07 

3.33 

.378 

.214 

2 

4 

17 

.41 

1.27 

6.42 

-1.59 

-2.73 

10.58 

3.66 

.420 

.159 

3 

1 

7 

.39 

1.22 

6.15 

-2.61 

.22 

-.85 

49.0 

1.253 

1.086 

15 

25 

18 

6.31 

19.11 

60.36 

-8.69 

-5.89 

42.36 

50.0 

1.169 

1.088 

15 

25 

40 

5.88 

17.88 

58.57 

-9.12 

-7.12 

- 

51.0 

.940 

.943 

10 

20 

40 

4.10 

12.66 

48.72 

-5.90 

-7.34 

- 

52.0 

.748 

.516 

6 

11 

40 

1.85 

5.78 

27.17 

-4.15 

-5.22 

- 

53.0 

.617 

.265 

5 

5 

9 

.90 

2.82 

14.03 

-4.10 

-2.18 

5.03 

54.0 

.790 

.  576 

5 

5 

10 

2.17 

6.76 

31.01 

-2.83 

1.76 

21.01 

55.0 

.908 

.666 

4 

10 

40 

2.87 

8.92 

38.51 

-1.13 

-1.08 

- 

56.0 

.824 

.814 

5 

5 

19 

3.10 

9.64 

40.  76 

-1.90 

4.64 

21.76 

20.0 

.632 

.699 

5 

5 

40 

2.05 

6.40 

29.64 

-2.95 

1.40 

- 

21.0 

.578 

.345 

5 

5 

5 

.99 

3.09 

15.31 

-4.01 

-1.91 

10.31 

22.0 

.630 

.363 

6 

15 

40 

1.14 

3.57 

17.58 

-4.86  - 

11.43 

- 

23.0 

.910 

.845 

10 

23 

40 

3.56 

11.04 

44.70 

-6.44 

-11.96 

- 

24.0 

1.178 

1.108 

16 

30 

41 

6.03 

18.31 

59.22 

-9.97 

-11.69 

- 

■ .  .0 

1.219 

.985 

14 

30 

40 

5.61 

17.10 

57.35 

-8.39 

-12.90 

- 

£,*>•  0 

1.158 

1.056 

15 

26 

18 

5.66 

17.25 

57.58 

-9.34 

-8.75 

39.58 

27.0 

.893 

699 

9 

20 

18 

2.  93 

9.12 

39.15 

-6.07 

-10.88 

21.15 

28.0 

.691 

.342 

6 

6 

40 

1.24 

3.87 

18.92 

-4.76 

-2.13 

20  Set  average  error  -3.77  -2.69  16.66 

Overall  average  error  3.4 
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Figure  26.  Current  Meter  Tilt  Angles  Computed  bv  Simple 
Statistics  Model  Compared  With  Observed  Data 
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4.1.2  WHOI  Mooring  No.  279 
37 

Millard  describes  tension  measurements  made  on  a  taut-moored  buoy 

system.  The  buoy  system  (figure  27)  was  installed  in  water  2685  m  deep  at 

Woods  Hole  Oceanographic  Institution  Site  D.  Tensions  were  recorded  at  four 

locations  along  the  mooring  line,  and  currents  were  recorded  at  a  depth  of  12  m 

79 

for  the  2h  months  that  the  buoy  system  was  on  station.  Berteaux  and  Walden 
describe  the  properties  of  the  wire  rope  and  plaited  nylon  rope  used  in  this  buoy 
system  as  follows: 

Diameter  (in. )  Weight/ft  in  Sea  Water  (lb/ ft) 

l/4-in.  1  x  50  wire  0.25  0.090 

rope 

5/8-in.  plaited  nylon  0.625  0.01047. 

7  2 

The  cable  modulus  for  the  wire  rope  was  taken  as  1.682  x  10  Ib/in.  . 

Nylon  rope  is  subject  to  both  elastic  and  inelastic  deformation  when  loaded. 

New  rope,  when  first  loaded,  will  acquire  a  permanent  deformation,  the  amount 

of  which  depends  on  the  initial  load.  Furthermore,  if  the  load  is  left  on  the 

rope,  the  rope  is  subject  to  creep  and  the  permanent  deformation  increases 
14 

with  time.  Martin  discusses  the  various  mechanisms  for  the  deformation  of 

nylon  rope.  Using  Martin's  curve  for  the  percent  stretch  versus  load  for  the 

5/8-in. -diameter  plaited  nylon  rope,  the  rope  modulus  is  computed  as  follows: 

E  =  3.52  x  105  lb/in.2  ;  0  <  T  <  1000  lb 
5  2 

E  ^  6.  79  x  10  lb/in.  :  1000  <  T  <  2000  lb 
6  2 

E  =  1.041  x  10  Ib/in.  ;  2000  <  T  . 
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PURPOSE  OF  TEST 


PROCEDURE 

EQUIPMENT 


STA.  279 
SITE  "D" 

SET:  OCT.  1, 1968 
RET:  DEC.  11, 1968 
DEPTH:  2685  M. 


-  EVALUATION  OF  MOORING  CONFIGURATION  AS  SHOWN  OVER 
A  TWO  MONTH  PERIOD  -  MEASUREMENT  OF  MOORING  TENSION 
STRETCH  13% 


-  LAUNCH  BUOY,  PAY  OUT  MOORING  LINE,  ATTACH  BALLS  , 
LAUNCH  ANCHOR,  CHECK  ANCHORING,  RETRIEVE  NEXT  CRUISE 

-  AS  SHOWN 


Figure  27,  Woods  Hole  Oceanographic  Institution  Mooring  No.  279 
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Webster  measured  currents  for  a  2-month  period  at  site  D  ,  Webster's 
data  were  curve-fitted  using  the  program  CURFIT  on  the  GSA  time-sharing 
computer,  and  the  following  function  for  the  current  strength  as  a  function  of 
depth  was  developed: 


C  =  2.6 


Os  .  D-0-418 


where 


C  is  the  current  (ft/sec) 

Cs  is  the  surface  current  (ft/sec) 

D  is  the  depth  (meters). 

The  steady-state  buoy  system  configuration  computer  program  was 
modified  to  include  Webster's  current  profile,  Martin's  elastic  properties  for 
the  nylon  rope,  and  the  cable  properties  given  by  Berteaux  and  Walden.  Since 
information  on  the  variations  of  current  direction  with  depth  was  not  available,  the 
currents  were  assumed  to  be  acting  in  the  same  direction  at  all  depths.  Further¬ 
more,  the  initial  inelastic  stretch  due  to  the  emplantment  and  the  dynamic  wave 
loads  is  not  known.  The  no-current  elastic  stretch  of  the  nylon  rope  was 
assumed,  arid  the  tensions  in  the  system  were  computed  while  currents  acted  on 
the  system.  Tensions  at  the  junction  of  the  wire  rope  and  nylon  rope  are 
shown  in  figure  28  as  functions  of  the  surface  current  and  of  the  no-current 
olastic  stretch.  Data  taken  by  Millard  at  the  same  location  on  the  mooring  line 
are  also  shown  on  figure  28.  Because  of  the  creep  properties  ot  the  nylon  rope, 
the  inelas'.ic  stretch  will  increase  ar  J  the  no-current  elastic  stretch  will 
decrease  as  time  increases.  The  shape  of  the  tension  curve  will  remain 
roughly  the  same,  but  the  tension  bias  will  decrease  with  time.  Since  the 


creep  properties  of  the  nylon  rope  are  not  known,  computation  of  errors 
between  observed  and  computed  tensions  in  this  case  is  of  little  value. 

Millard  analyzed  the  tension  versus  surface  current  data  shown  on  ■ 
figure  28  and  least-squares  fitted  the  linear  function: 

T  =  12.3  •  Cs  +  369  . 

Least-squares  linear  fits  to  the  computed  curves  A,  B,  and  C  are 
T  a  12.59  •  Cx  +  608  , 

T  =  12.25  ♦  Cs  +  381  , 
and 

T  =  12. 10  •  Cs  +  203  . 

Slopes  of  the  computed  linearized  functions  are  in  very  close  agreement,  which 
indicates  that  the  functional  form  of  the  computed  tensions  is  accurate. 

From  this  comparison,  it  is  obvious  that  more  experimentation  on  the 
elastoplastic  properties  of  nylon  rope  is  needed  in  order  to  predict  the  steady- 
state  configurations  of  deep  ocean,  taut- moored  buoy  systems.  Also,  a  deep 
ocean  buoy  system  should  be  installed  with  both  recording  tensiometers  and 
inclinometers  to  belter  validate  the  steady-state  computer  model. 


32-ft'diameter  sphere,  was  installed  in  March  1970  off  Great  Salt  Pond 
entrance,  Block  Island,  Rhode  Island.  The  larger  buoy,  an  8-ft-diameter 
torroid  was  installed  at  station  BRAVO  during  May  1970.  A  description  of  the 
instrumentation  and  a  discussion  of  the  measurements  taken  are  presented  in 
reference  81. 

4.2.1  Spherical  Buoy  at  Station  D 

The  3|-ft-diameter  spherical  buoy  is  shown  in  figure  29.  This  buoy  was 
loaned  to  the  writer  by  Dr.  A.  Nalwalk,  of  the  Marine  Sciences  Institute  of  the 
University  of  Connecticut.  The  buoy  was  equipped  with  the  following  instru¬ 
mentation: 

Heave  motion  statistical  accelerometers 

Current  meter 

Heave  accelerometer 

Surge  Accelerometer 

Sway  accelerometer 

Pitch  pendulous  potentiometer 

Roll  pendulous  potentiometer 

Cable  pitch  pendulous  potentiometer 

Cable  roll  pendulous  potentiometer 

Cable  tension  gage 

Thermistor. 

The  buoy  system  was  installed  in  62  ft  of  water,  west  of  Great  Salt  Pond 
entrance  at  Block  Island  by  the  Research  Vessel,  UCONN,  on  2  March  1970. 
The  lg-in.-diameter,  14-conductor  armored  cable  was  laid  along  the  bottom  to 
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the  beach  and  was  carried  over  the  beach  to  the  BIFI  field  station:  A  Snodgrass 
wave  sensor  was  also  installed  off  the  beach  in  25  ft  of  water.  The  arrangement 
of  cables  and  instruments  is  shown  schematically  in  figure  30. 

Two  weeks  of  statistical  heave  accelerometer  data  were  collected.  Some 
buoy  motion  data  were  also  recorded.  The  electrical  conductors  in  the  cable 
began  to  fail  after  a  month  of  use  and  a  number  of  attempts  were  made  to  repair 
the  cable  in  order  to  continue  collecting  data.  The  buoy  broke  loose  during  a 
storm  in  early  November  1970  and  is  missing. 

The  statistical  accelerometers  were  designed  to  count  at  0.35-,  0.50-,  0.65-, 
0.  80- ,  1.20-,  1.35“,  1.50-,  and  1. 65-g  heave  acceleration  levels.  The  pulses 
from  the  accelerometers  actuated  counters  in  the  van  at  the  BIFI  Field  Station, 
and  the  total  count  was  recorded  daily.  Sea  state  data  were  based  on  estimates 
by  the  resident  engineer  (Carl  T.  Milner)  and  by  Coast  Guard  observations 
reported  by  the  ESSA  Marine  Weather  Service.  Data  were  recorded  during  the 
period  4  March  to  12  March  1970. 

Data  from  the  22  observations  of  positive  acceleration  counts  were  used 
to  generate  figure  31,  a  plot  of  the  counts  per  hour  versus  acceleration  level 
for  various  sea  states.  The  negative  acceleration  counters  did  not  work,  because 
of  the  failure  of  leads  in  the  cable.  Figure  31  indicates  that  for  any  sea  state 
greater  than  sea  state  0,  the  buoy  will  always  undergo  1.2-g  accelerations  at  a 
rate  of  1,000/hr.  The  number  of  cycles  per  hour  for  higher  acceleration  levels 
will  increase  with  increasing  sea  state.  Figure  32  is  the  conditional  probability 
of  the  buoy  exceeding  various  positive  acceleration  levels  given  that  the  buoy 
exceeds  1.2-g  level  accelerations.  Figure  33  shows  the  buoy  heave  acceleration 
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Figure  30.  Spherical  Buoy  System  at  Block  Island 
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amplitude  histogram  for  various  sea  state  conditions.  Rayleigh  distributions 

are  plotted  over  the  histograms  and  were  computed  by  using  the  mean  of  the 

55  51 

histograms.  Longuet- Higgins  and  Bretschneider  have  shown  that  the 
distribution  of  wave  heights  is  given  by  the  Rayleigh  distribution. 

Analog  data  of  the  spherical  buoy  motions  were  also  recorded.  On 
16  March  1970,  the  following  records  were  obtained  on  a  two-channel  strip 
chart  recorder: 

Wave  height  (20  min) 

Buoy  heave  and  buoy  surge  (15  min) 

Buoy  heave  and  buoy  pitch  (10  min) 

Buoy  heave  and  buoy  roll  (10  min) 

Buoy  heave  and  cable  pitch  (10  min) 

Buoy  heave  and  cable  roll  (10  min). 

Winds  were  15  to  20  knots,  northeast,  with  an  estimated,  sea  state  3  at  the 
buoy.  Results  of  a  simple  "quick  look"  analysis  are  shown  in  figures  34  through 
37. 

One- hundred  samples  of  each  record  were  digitizeu  and  analyzed  on  the 
GSA-440  time  sharing  computer.  Figure  34  shows  means,  variances,  and 
standard  deviations  for  each  parameter.  In  addition,  the  correlation  matrix  for 
simple  product-moment  correlations  is  shown.  Parameters  that  should  be 
coupled  appear  to  be  coupled,  and  parameters  that  should  be  decoupled  appear 
to  be  decoupled  For  example,  heave-surge,  heave-pitch,  and  surge-pitch  are 
coupled,  and  heave-roll,  surge-roll,  and  pitch-roll  are  decoupled.  Also,  cable 
angles  are  mildly  coupled  to  buoy  displacements  but  decoupled  from  buoy  angles. 
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Figure  34.  Spherical  Buoy  Motion 
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VARIANCE 
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DEVIATION 

0.0242  g2 

0.1555  g 

0.00978  g2' 

0.0991  g 

104. 2  deg2 

10. 20  deg 

55.3  deg2 

7.58  deg 

24. 0  deg2 

4. 90  deg 

16. 2  deg2 

4. 01  deg 

ROLL 

CABLE 

PITCH 

CABLE 

ROLL 

-0.0706 

0.2430 
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Figure  36.  Spherical  Buoy  Motion  Amplitude  Histograms  —  Pitch  and  Roll 
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One  hundred  amplitudes  from  each  parameter  record  were  digitized  and 
were  used  to  generate  parameter  amplitude  histograms.  These  probability 
distributions  are  shown  in  figures  35  and  37.  The  median  amplitudes  were 
computed  and  were  used  to  compute  Rayleigh  distributions,  which  are  plotted 
over  the  histograms.  The  Rayleigh  distributions  of  the  form 


pW)  =  £  ji  e 


_  7T  H* 

r7F 


where 


p(H)  is  the  probability  of  parameter  H 
H  is  the  parameter 
H  is  the  mean  value  of  the  parameter, 
were  found  to  match  the  histograms  quite  well.  The  fact  that  the  amplitude 
probability  distributions  all  appear  to  fit  a  Rayleigh  distribution  indicates  a 
linear  transform  from  wave  height  to  buoy  response. 

If  the  functions  relating  the  mean  buoy  motion  amplitude  parameters  to 
sea  state  were  known,  the  probability  distribution  for  any  parameter  in  any  sea 
state  can  be  computed  from  the  Rayleigh  distribution. 

From  the  statistical  accelerometer  data,  the  mean  heave  acceleration 

80 

amplitude  can  be  plotted  versus  mean  wave  height  (using  Vine  and  Volkman's 
relations  for  sea  state  and  mean  wave  height).  This  curve,  shown  in  figure  38, 
can  be  approximated  in  the  region  H  =  1  to  10  ft  with  the  linear  function; 

Hv  =  o.u  ■+  o.  ooze  H  . 

Assuming  a  linear  transform  for  the  other  motion  parameters  and  averaging 
pitch  and  roll  amplitude  means,  we  can  write  a  set  of  linear  equations  for  the 


«sgwg 
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mean  motion  amplitudes.  In  addition,  these  equations  can  be  substituted  into 
the  Rayleigh  equation  to  find  the  amplitude  probability  distribution  for  any 
parameter.  The  linear  equations  and  their  amplitude  distributions  are  as 
follows: 


For  1  <  H  <  10  ft,  the  mean  heave  acceleration  amplitude  is 

Hv  =  0.11  +  O.oozeH, 

The  amplitude  probability  distribution  is 

p(Hv)=  £ , 

The  mean  surge  acceleration  amplitude  is 


Su  “  0, 06 XX  +  0.  00537  H  , 

and  the  probability  distribution  is 


p(3u) 

The  mean  sway  acceleration  amplitude  is 


Sw  -  0.06S1+  0-0053 7H, 

and  the  probability  distribution  is 

ptew)=  f &e  f  * 

The  mean  pitch  angle  amplitude  is 


Pt  ~  5.51  +  0.H3G  H  } 

and  the  probability  distribution  is 

P(Pe)=£#e-^ 

The  mean  roll  angle  amplitude  is 

RL  =  s.?1?  +  aw  H  , 

and  the  probability  distribution  is 


* 


r. 
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;(RL)  «  I  Jt 


7T  ftt 

5tp 


piA<-;=  •r-e=ze 

The  mean  cable  pitch  angle  amplitude  is 

CjH  =  3.HI  +  0. 2<7V  H  , 


and  the  probability  distribution  is 


p(Cpt)=me 


7T  CP6C 

“T  w 


ZCF? 

The  mean  cable  roll  angle  amplitude  is 

CM  *  3.‘//  +  0.2‘?‘<  H  > 

and  the  probability  distribution  is  „ 

7r  CRL2 

p(c«o-  f  $,  e  »■. 


The  preceding  empirical  equations  will  serve  as  a  first-order  approxi¬ 


mation  to  the  buoy  motion  parameters  for  the  spherical  buoy  and  can  be  used 
for  design  purposes. 


4.2.2  Torroidal  Buoy  at  Station  BRAVO 

An  8-ft-diameter,  torroidal  oceanographic  buoy  was  outfitted  with  buoy 
motion  sensing  instrumentation  and  telemetry  and  was  installed  in  Block  Island 
Sound  during  May- June  1970.  The  buoy  (figure  39)  was  installed  by  the  USCGC 
MARIPOSA  on  30  April  1970  at  station  BRAVO  in  120  ft  of  water.  A  telemetry 
receiving  station  was  established  in  the  generator  building  at  the  Watch  Hill 
Lighthouse,  Watch  Hill,  Rhode  Island  (figures  40  and  41).  A  schematic  of  the 
instrumentation  arrangement  is  shown  in  figures  42  and  43.  A  detailed  descrip¬ 
tion  of  the  instrumentation,  circuitry,  and  calibrations  is  described  in 
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The  buoy  contained  >e  following  motion  sensing  instrumentation: 

Heave  accelerometer 
Surge  accelerometer 
Sway  accelerometer 
Pitch  angle  potentiometer 
Roll  angle  potentiometer. 

In  addition,  a  self-recording  current  meter  was  attached  to  the  mooring  cable 
a'  a  depth  of  60  ft.  Wind  speeds  were  recorded  on  a  paper  tape  recorder  at  the 
Watch  Hill  Lighthouse.  Wind  directions  were  logged  every  4  hr  by  the  duty 
personnel  at  the  Lighthouse.  Also,  the  ESSA  weather  reports  on  163.5  kHz 
were  monitored  during  each  data  run,  and  reported  conditions  at  Coast  Guard 
Stations  bordering  Block  Island  Sound  were  logged  on  the  data  sheets. 

The  buoy  transmitted  data  every  12  hr  for  a  ^-hr  time  period.  Each  data 
transmission  was  preceded  by  a  calibration  sequence  consisting  of  two  voltage 
levels  from  the  potentiometer  sensors.  The  FM  signals  from  the  buoy  con¬ 
taining  the  five  mixed  frequencies  were  demodulated  and  recorded  on  an  FM 
tape  recorder.  The  composite  signal  and  a  12.5-kHz  phase-lock  signal  were 
also  recorded  on  two  AM  channels. 

During  the  51  days  the  buoy  was  on  station,  55  data  transmissions  were 
recorded.  The  only  major  problem  encountered  was  the  failure  of  the  opera¬ 
tional  amplifier  that  mixed  the  five  frequencies  from  the  voltage  controlled 
oscillators  in  the  buoy.  This  component  was  replaced  and  tests  continued. 
Structurally,  the  only  failure  was  the  loss  of  a  cotter  pin  on  a  shackle  holding 
one  of  the  three  chain  bridle  legs  under  the  buoy.  The  loss  of  this  pin  allowed 
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the  shackle  to  come  undone  in  such  a  manner  that  the  buoy  bridle  had  only  two 
active  legs.  The  buoy  was  recovered  on  19  June  1970  by  the  USCGC  REDWOOD. 

Because  of  the  large  amount  of  recorded  data,  a  simple  "quick  look" 
analysis  of  the  data  was  performed  by  playing  back  the  recorded  data  for  each 
of  the  five  motion  parameters  on  a  "memoscope. "  A  2-min  sample  of  each 
parameter  for  each  run  was  traced  out.  The  mean  width  of  the  band  trace-out 
was  measured,  and  the  value  of  the  parameter  double  amplitude  (in  volts)  was 
logged.  This  value  was  assumed  to  be  the  "significant"  amplitude,  i.e. ,  the 
mean  of  the  hig>  it  one-third  of  the  amplitudes.  It  was  found  that  the  data  on 
the  pitch  channel  were  too  noisy  to  be  used  in  this  fashion;  thus,  pitch  was  not 
included  in  this  analysis.  The  logged  values  of  buoy  heave,  surge,  and  sway 
acceleration  amplitudes  along  with  buoy  roll  angle  amplitude  are  shown  in 
table  4,  The  environmental  conditions  are  also  shown  for  each  run. 

The  statistics  of  the  measured  parameters  are  shown  in  figure  44.  The 
matrix  of  simple  product- moment  correlation  coefficients  is  also  shown.  These 
statistics  were  computed  on  the  GSA-440  time  sharing  computer.  Inspection  of 
the  correlation  matrix  indicates  that  all  the  significant  buoy  motions  are  well 
correlated.  In  addition,  they  correlate  well  with  wii.d  speed  and,  to  a  lesser 
extent,  with  computed  wind  wa^e  height.  A  similar  analysis  was  conducted  with 
current  speed.  The  elements  of  that  correlation  matrix  are  as  follows: 

Heave  Sway  Surge  Roll 
Current  0.1537  -0.00572  0.1040  -0.0237. 

Thus,  it  appears  that  buoy  motion  amplitudes  are  not  correlated  with  current 
speed.  Buoy  motions  are  not  significantly  affected  by  current  speed. 
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Figure  44.  Torroidal  Buoy  Motion  Statistics  for  6-week  Period 


The  GSA-440  computer  was  also  used  to  fit  curves  to  the  observed 
parameter  significant  amplitudes.  Six  curves  were  fitted  using  least- squares 
methods. 

In  general,  the  linear  plot  had  the  best  index  of  determination  when  fitted 
to  the  data.  Figures  45  through  50  show  the  observed  buoy  motion  significant 
amplitudes  plotted  versus  wind  speed,  computed  significant  wave  height,  and 
buoy  heave  acceleration.  The  least-mean-squares  linear  curve  is  shown  on 
each  plot.  The  least-mean-squares  linear  functions  were  transformed  to 
engineering  units  and  are  summarized  in  figure  51.  If  the  buoy  motion  ampli¬ 
tudes  are  assumed  to  be  distributed  by  the  Rayleigh  distribution,  the  mean 
amplitude  is  62g  percent  of  the  significant  amplitude.  The  empirical  equations 

shown  in  figure  51  were  again  transformed  in  such  a  manner  that  they  are  in 
> 

terms  of  mean  amplitudes.  They  are  shown  in  figure  52. 

The  cumulative  probability  distributions  for  wind,  current,  observed  and 
computed  wave  height,  and  buoy  motion  amplitudes  are  shown  in  figures  53 
through  55  for  the  time  period  that  the  buoy  was  on  station. 

A  complete  spectral  analysis  was  performed  on  two  runs  of  buoy  motion 
data  taken  on  10  and  11  June  1970.  A  Fast  Fourier  Transform  method  that  is 
programmed  and  is  available  on  the  NUSC  UNIVAC  3108  computer  was  used  to 
compute  power  spectra  and  cross  correlations.  The  data  were  digitized  at  a 
rate  of  64  samples  per  second  on  automatic  data  processing  equipment  by  the 
Data  Analysis  Branch  at  NUSC.  A  set  of  buoy  motion  spectra  for  the  second 
data  run,  which  was  taken  at  2030  EDST,  on  11  June  1970,  is  shown  in  figures  56 
through  61.  The  wind  was  at  10  knots  from  the  southwest  and  had  been  blowing 
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Figure  46.  Buoy  Motions 


versus  Wind  Speed  -  Surge  and  Boll 
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1.  ON  WIND  {Win  knots) 

HEAVE  ACCELERATION 

SWAY-SURGE  ACCELER¬ 
ATION. 

PITCH-ROLL  ANGLE 


Hvl/3  =  °*0064  +  0-0163*W  (g's) 

BUl/g  =  Sw1/3  =0.1056+  0. 003825. W  (g's) 

Pil/3  =  Rl1/3  ---  12. 83°  +  0. 1373* W  (deg) 


2. 


ON  SIGNIFICANT  WAVE  HEIGHT  (Hl/3  in  ft) 


HEAVE  ACCELERATION 


HVi/g  =  °* 168  +  °* 01792  *Hi/3  <g's) 


SWAY-SUR GE  ACCELER-  Su  .  =  Sw.  ,  =  0. 120  +  0. 01015.H,  .  (g's) 

AT  ION  1/3  1/3  1/3 


PITCH-ROLL  ANGLE 


Pil/3=R1l/3  =  12*36  +  1-°-Hl/3  <d6g) 


3.  ON  HEAVE  ACCELERATION  AMPLITUDE  (Hv^3  in  g's) 

SWAY-SURGE  ACCELER-  Su  .  =  Svv  ,  =  0. 0268  +  0. 557. Hv  (g's) 

AT  ION  1/3  1/3  1/3 


fCH-ROLL  ANGLE 


Pil/3  =  R1l/3  ’  7,82  +  31-35*Hv1/3  <deS) 


Figure  51.  Empirical  Equations  cor  Significant  Buoy  Motion  Amplitudes 
(Based  on  Linear,  Least-Mean-Squares  Equations) 
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1.  ON  WIND  (Win  knots) 

HEAVE  ACCELERATION  Hv  =  0.004  +  0.0102®W  (g's) 

SWAY-SURGE  ACCELERATION  Su  =  Sw  =  0. 066  +  0. 002394  ®  W  (g’s) 

PITCH-ROLL  ANGLE  Pi  =  Rl  =  8. 03  +  0. 086  ®W  (deg) 


2.  ON  MEAN  WAVE  HEIGHT  (H  in  ft) 

HEAVE  ACCELERATION  Hv  =  0. 105  +  0. 01792  ®H  (g's) 

SWAY-SURGE  ACCELERATION  Sw  =  Su  =  0. 075  +  0. 01015  ®H  (g's) 

PITCH-ROLL  ANGLE  Pi  =  Rl  =  7.  72°  +  1.0®  H  (deg) 

(NOTE  THAT  THE  PROBABILITY  DISTRIBUTIONS  ARE  THE  SAME 
AS  ABOVE.) 

3.  ON  MEAN  HEAVE  ACCELERATION  AMPLITUDE  (Hv  in  g's) 

SWAY-SURGE  ACCELERATION  Sw  =  Su  =  0.01674  +  0.557® Hv  (g's) 

PITCH-ROLL  ANGLE  Pi  =  Rl  =  4.  89  +  31.  35®  Hv  (deg) 


Figure  52.  Empirical  Equations  for  Mean  Buoy  Motion  Amplitudes 
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steady  at  that  strength  for  1  hr.  The  wind  wave  spectrum  was  computed  by 
using  a  two-parameter  Bretschneider  spectrum.  Figure  56  shows  the  computed 
wind  wave  power  spectrum.  Figures  57  through  61  show  the  buoy  motion  power 
spectra.  Figure  57,  the  buoy  heave  acceleration  power  spectra,  shows  a  peak 
in  the  energy  around  0.5  rad/sec  (T  =  12.5  sec).  There  is  also  a  peak  around 
1.5  rad/sec,  which  corresponds  to  the  peak  of  the  computed  wind  wave  power 
spectrum. 

The  peak  at  the  lower  frequency  may  be  caused  by  swell;  a  1— ft  swell 
from  thq  south  was  observed  visually  during  the  data  run.  The  same  type  of 
energy  distribution  is  seen  in  the  other  spectra.  Buoy  roll  angle  and  sway 
acceleration  spectra  indicate  relatively  more  energy  at  wind  wave  frequencies; 
thus,  it  is  concluded  that  the  buoy  was  oriented  in  such  a  manner  that  the  sway- 
heave  plane  was  close  to, a  southwest-northeast  orientation. 

Buoy  motion  parameters  were  cross-correlated  with  buoy  heave  accelera¬ 
tion  and  are  shown  in  figures  62  through  65.  These  plots  indicate  the  existence 
of  two  modes  with  rather  strong  coupling  between  motion  parameters.  If  a 
linear  system  is  assumed,  the  input  and  output  spectra  can  be  related  by  the 
transfer  function: 

5Aa(W)  =  H  (*>)  •  SHz(w)  } 

where 

the  buoy  motion  power  spectrum 
is  the  wave  height  power  spectrum 

ffM  is  the  transfer  function  . 
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Wave  spectra  determined  by  Williams  for  similar  conditions  in  Block 
Island  Sound  indicate  that  swell  appears  in  the  wave  spectrum  as  a  horizontal 
line  to  the  left  of  the  wind  wave  peak  and  is  about  5  dB  down  from  the  peak. 

Figure  56  was  modified  to  include  swell.  The  buoy  motion  power  spectra  were 
smoothed  and  the  difference  in  decibels  between  them  and  the  wave  spectrum 
(including  swell)  were  plotted  as  the  square  of  the  absolute  value  of  the  transfer 
function  in  decibels.  Figures  66  through  70  show  the  transfer  functions  and 
indicate  the  existence  of  two  modes.  The  relative  magnitudes  of  the  peaks 
indicate  .that  the  modes  are  a  heave  mode  (low  frequency)  and  a  roll  mode  (high 
frequency).  Observations  made  by  the  writer  while  servicing  the  buoy  at  sea 
indicate  that  the  sway-surge  mode  has  the  lowest  frequency,  the  heave  mode 
has  the  next  highest  frequency,  and  the  pitch-roll  mode  has  the  highest  frequency. 

These  spectra  indicate  that  the  heave,  surge,  and  sway  motions  of  the 
buoy  are  primarily  excited  by  the  sea  swell,  whereas  buoy  pitch  and  roll  are 
excited  by  the  higher  frequency  wind  waves. 

Analysis  of  measured  buoy  motion  data  for  the  3^- ft  spherical  buev  and 
for  the  8-ft  torroidal  buoy  have  yielded  sets  of  empirical  equations  that  can  be 
used  to  predict  mean  buoy  motions  and  amplitude  distributions  for  various  sea 
states.  Buoy  motion  amplitude  distributions  appear  to  follow  a  Rayleigh 
distribution,  which  indicates  that  the  dynamic  system  is  linear  or  near-linear. 

Spectra  for  the  torroidal  buoy  indicate  the  existence  of  a  heave  mode  and 
a  pitch-roll  mode  within  the  range  of  wave  frequencies.  In  general,  both  buoys 


are  hard-coupled  to  the  sea  surface  in  heave.  The  torroidal  buoy  appears  to 
respond  in  pitch- roll  motions  to  a  greater  extent  than  does  the  spherical  buoy. 
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*4^  4 . 3  Simulation  and  Comparison  of 

Buoy  System  Dynamics 

4.3. 1  Spherical  Buoy  at  Station  D 

Steady- state  configurations  of  the  spherical  buoy  system,  installed  at 
Station  D,  are  shown  in  figure  71  for  various  uniform  currents.  Very  little  of 
the  3/4- in.  DiLock  chain  is  picked  up  off  the  bottom  even  under  the  worst  current 
conditions  —  1  knot.  Steady -state  tensions  never  go  over  100  lb,  and  the  watch 
circle  radius  can  vary  from  38  ft  at  0  knot  to  62  ft  at  1  knot.  These  configura¬ 
tions  were  computed  nn  the  GSA  time-sharing  computer  using  the  steady-state 
buoy  system  configuration  program  shown  in  appendix  B.  The  Savonius  rotor 
current  meter  used  in  the  spherical  buoy  system  failed  3  days  after  emplantment, 
thus;  cable  angles  can  not  be  correlated  with  current  strength  for  this  buoy 
system.  However,  the  computed  configurations  are  used  as  initial  conditions 

for  the  buoy  system  dynamics  simulation  when  current  strength  and  direction 

83 

are  computed  from  the  Coast  and  Geodetic  Survey  current  tables. 

The  lumped- mass  dynamic  eouations  of  motion  for  the  cable  and  the 
equations  of  motion  for  the  spherical  buoy  were  programmed  and  solved  numer¬ 
ically  in  the  time  domain  using  a  fourth-order, Iiunge-Kutta  numerical  integra¬ 
tion  scheme.  Five  mooring  line  mass  elements  were  used  —  three  for  the 
cable  and  two  for  the  chain.  The  buoy  was  allowed  six  degrees  of  freedom 
(heave,  surge,  sway,  yaw,  pitch,  and  roll),  and  each  mooring  line  mass 
element  was  allowed  three  degrees  of  freedom  (x.  ,  y  ,  and  z).  The  program 
shown  in  appendix  B  can  accept  wind  and  current  vectors  coming  from  any 
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Figure  71.  Computed  Steady-State  Configurations  of  the  Spherical  Buoy  System 


214 


«!_?;  direction;  however,  the  waves  are  assumed  to  be  two-dimensional  and  are 

constrained  to  come  in  on  the  z  axis  of  the  inertial  coordinate  system.  A 
complete  listing  of  the  input  data  needed  to  describe  the  spherical  buoy  system 
is  shown  in  appendix  D.  Mean  wave  heights  and  periods  were  computed  from 
reference  82  for  various  sea  state  conditions.  The  spherical  buoy  system 
dynamics  model  was  excited  with  the  ten-component  random  wave  model  based 
upon  a  two-parameter  Bretschneider  spectrum  having  the  mean  wave  height  and 
period  for  each  sea  state.  The  average  wind  strength  and  direction  and  the 
average  current  strength  over  the  time  period  in  which  the  statistical  acceler¬ 
ometers  were  in  operation  were  also  used  to  force  the  model.  The  wave 
amplitudes  for  each  component  were  allowed  to  build  linearly  over  two  compo¬ 
nent  wave  periods.  This  procedure  minimized  transient  motions.  The  solution 
was  allowed  to  proceed  in  time  as  transients  decayed.  Finally,  buoy  heave 
accelerations  were  sampled  over  a  time  period,  and  the  mean  heave  acceleration 
amplitude  was  computed. 

These  computed  amplitudes  are  shown  in  figure  72  and  are  compared  with 
the  amplitudes  derived  from  the  statistical  accelerometer  data.  The  computer 
model  overestimates  the  heave  acceleration  amplitudes  at  the  lower  sea  states 
and  agrees  quite  well  at  the  higher  sea  states.  Since  the  observed  sea  states  are 
based  upon  the  Block  Island  resident  engineer's  visual  observations,  a  plus  or 
minus  one  sea  state  error  band  is  shown  in  figure  72  for  the  observed  data.  It  is 
concluded  that  the  computer  model  provides  a  conservative  estimate  of  buoy 
accelerations  for  design  purposes. 


MEAN  HEAVE  ACCELERATION  AMPLITUDE  (g’s) 
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On  16  March  1970.  strip  chart  recordings  of  the  buoy  motion  instrumenta¬ 
tion  were  made.  The  computed  Bretsohneider  wind  wave  spectrum  for  that 


location,  date,  and  time  based  upon  the  observed  wind  speed,  direction,  and 
duration  is  shown  in  figure  73.  This  spectrum  was  used  to  determine  compo¬ 
nent  amplitudes  for  the  random  sea  forcing  function.  The  observed  wind  and 
the  current  as  computed  from  the  tidal  current  tables  were  also  used  to  force 
the  model  on  the  UNIVAC  1 108.  Computed  buoy  heave,  sway, and  surge  accel¬ 
erations,  pitch,  roll,  cable  pit-’h,  and  cable  roll  were  sampled  in  the  same 
manner  as  the  observed  data  and  were  analyzed  on  the  GSA  time-sharing 
-  computer. 


The  reduced  spherical  buoy  motion  parameters  based  upon  computed 
motions  and  observed  motions  are  shown  on  figures  74  and  75,  respectively. 
Since  the  spherical  buoy  was  not  equipped  with  a  yaw  sensing  device,  the  motions 
designated  surge  and  sway  and  pilch  and  roll  are  not  known  relative  to  the  z 
axis,  along  which  the  waves  are  traveling.  Visual  observations  made  from  the 
beach  while  the  data  were  being  recorded  indicated  that  the  buoy  was  aligned  in 
one  direction  with  little  or  no  yaw  motion.  A  comparison  of  figures  34  and  75 
indicates  that  the  observed  sway  may  really  have  been  surge  and  that  pitch  and 
roll  should  be  interchanged.  Figure  75  reflects  these  changes. 

A  comparison  of  figures  74  and  75  indicates  that  the  computer  model 

t , 

predict?  buoy  heave  accelerations  with  good  engineering  accuracy  (-9.5  percent 


error  for  heave  acceleration  standard  deviation),  underestimates  buoy  sway 
accelerations  (-71  percent  error  for  sway  acceleration  standard  deviation),  and 
overestimates  buoy  pilch  and  roll  motion  (+5*1  percent  error  for  pilch  angle  and 
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Figure  74.  Simulated  Spherical  Buoy  Motion  Parameter  Statistics 
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Figure  75.  Observed  Spherical  Buoy  Motion  Parameter  Statistics 


*74  percent  error  for  roll  angle  standard  deviations).  Mooring  line  pitch  and 
roll  angle  standard  deviations  are  underestimated  by  —33.  6  percent  and 
-19.0  percent,  respectively. 

The  simple  product-moment  correlation  matrix  for  the  observed  data 
indicates  moderate  coupling  for  heave-sway,  heave- roll,  sway- roll,  and  pitch- 
cable  pitch.  Because  of  the  axial  symmetry  of  the  buoy,  it  is  expected  that 
heave-pitch  and  surge-pitch  would  also  bo  coupled.  The  simulated  buoy  motion 
correlation  matrix  indicates  coupling  between  heave-sway,  sway-surge,  sway- 
pitch,  swav-roll,  and  pitch-roll.  The  modal  coupling  indicated  bv  the  two 
correlation  matrices  do  not  agree.  This  poor  agreement  is  probably  due  to  the 
rather  short  sample  time  used  to  compute  the  product  moment  correlations.  An 
error  estimate  based  upon  an  assumed  bandpass  while  noise  sea  spectrum  with 
a  bandwidth  ol  1  11/  indicates  normalized  errors  of  about  22  percent  in  the  stand¬ 
ard  deviations  and  about  18  percent  in  the  product-moment  correlations.  The 
comparison  of  simulated  I.u<>v  motions  w  th  observed  buoy  motions  for  this  case 


indicates  reasonable  ngivonu  nt  for  parameter  standard  deviations  but  poorer 


agreement  lor  product  moment  correlations.  In  view  of  the  limited  amount  of 


data  and  its  rclaliwm  poor  'pialih.  this  particular  comparison  will  not  be 


extended. 


1.8.2  Torroiilal  Buoy  at  Station  BRAVO 

The  motion  computer  model  with  lumped- mass  cable  elements  was 
inoddied  and  input  values  were  changed  in  order  to  simulate  motions  of  the 
torroidal  bu<>\  BBAVO.  A  subroutine  to  compute  the  buo\ant  forces  and  moments 
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by  polar  integration  around  the  torroid  as  a  function  of  its  draft  and  lilt  angle 
was  incorporated  in  the  program  (appendix  E).  A  subroutine  to  compute  simple 
statistical  properties  (mean,  variance,  and  standard  deviation)  of  the  input  wave 
height  and  output  buoy  motions  (heave,  surge,  and  sway  accelerations  ami  pitch 
and  roll  angles)  was  also  incorporated  in  the  program. 

Winds  were  assumed  to  act  from  the  southwest  (the  predominant  wind 
direction  during  May  and  June),  and  an  average  ebb  current  of  0. 7  knot  was 
assumed  (uniform  in  depth,  setting  to  the  east).  Thus,  the  buoy  system  coor¬ 
dinate  system  has  the  z  axis  pointing  southwest  and  the  y  axis  pointing  south¬ 
east  (figure  76).  The  S-M-B  method  was  used  to  compute  the  mean  wave  height 
and  period  for  winds  of  5,  10,  15,  20,  and  25  knots  for  the  southwest  winds 
with  an  assumed  duration  of  4  hr.  The  computed  mean  wind  wave  heights  and 
periods  for  station  BRAVO  are  shown  in  figure  77.  Resulting  Bretschneider 
spectra  for  these  conditions  are  shown  in  figure  78.  These  spectra  were  used 
to  compute  random  wa\e  component  amplitudes  and  frequencies  which, in  turn, 
forced  the  buoy  motion  computer  model. 

Initial  runs  with  a  five-element  cable  model  were  found  to  be  very  time- 
consuming  in  machine  time  since  the  relatively  low  mass  and  high  elastic 
modulus  of  the  cable  elements  required  that  numerical  integration  step  sizes  on 
the  order  of  0.001  sec  be  used  for  numerical  stability.  A  step  size  of  0.0005  sec 
was  used  for  accuracy.  TIk*  program  was  rewritten,  and  inputs  were  recom¬ 
puted  for  a  three-element  mooring  line  model.  Caole  masses  were  concentrated 
at  the  current  meter  in  the  middle  of  the  cable,  and  the  lengths  of  heavv  anchor 
chain  were  broken  up  to  form  the  other  two  mooring  line  lumped  masses.  This 
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Figure  78.  Computed  Wind  Wave  Spectra 
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procedure  greatly  reduced  the  highest  natural  frequency  of  the  system  and  the 
step  size  could  be  increased  to  0. 01  sec  before  numerical  instability  occurred. 

A  step  size  of  0. 005  sec,  which  allowed  an  order  of  magnitude  increase  in  compu¬ 
tational  speed, was  used.  No  significant  difference  was  seen  in  computed  buoy 
motions  when  using  the  three-lump  or  five-lump  cable  model. 

The  run  procedure  is  as  follows.  First,  the  model  is  acted  upon  by  the 
mean  wind  and  current  components  (no  waves)  and  allowed  to  converge  to  its 
steady-state  configuration.  Then,  the  random  wave  components  are  introduced 
and  buoy  motions  sampled  at  every  integration  for  60  sec.  Each  run  required 
about  20  min  of  computer  time  on  the  UNIVAC  1108.  The  ratio  of  computer 
time  to  solution  time  was  9. 15:1. 

The  results  of  four  runs  with  mean  wind  speeds  of  5,  10,  15,  and  20  knots 
are  shown  in  figures  79,  80,  and  81.  The  least-mean-squares  plots  of  the 
observed  data,  as  shown  in  figures  45  through  50  and  as  summarized  in 
figure  51,  are  also  shown  in  figures  79,  80,  and  81.  Since  the  yaw  orientation 
of  the  buoy  was  not  known  during  the  at-sea  measurements,  the  coefficients  of 
the  surge  and  sway  empirical  functions  and  of  the  pitch  and  roll  empirical 
functions  w’ere  averaged.  Also,  since  the  winds  and  waves  in  the  computer  model 
are  acting  along  the  z  axis,  only  heave,  surge,  and  pitch  motions  are  compared 
with  observed  data.  The  s  verage  error  indicated  in  the  simulation  over  the 
range  of  wind  speeds  considered  are  as  follows: 

1.  On  Mean  Wind  Speed 

Mean  heave  acceleration  amplitude  +14.45  percent 
Mean  surge  acceleration  amplitude  -63.6  percent 


SURGE  ACCELERATION  AMPLITUDE  ( 
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Mean  pitch  angle  amplitude  -42.4  percent 

2.  On  Mean  Wave  Height* 

Mean  heave  acceleration  amplitude  +2.22  percent 
Mean  surge  acceleration  amplitude  -60.15  percent 

Mean  pitch  angle  amplitude  -41.43  percent 

3.  On  Mean  Heave  Acceleration  Amplitude 

Mean  surge  acceleration  amplitude  -60.4  percent 
Mean  pitch  angle  amplitude  -42.3  percent. 

Again,  assuming  bandpass-limited  white  noise  spectra  of  about  1-Hz  bandwidth, 
note  that  the  error  in  the  standard  deviations  of  the  simulated  buoy  motion  is  about 
+13  percent.  Also,  known  and  estimated  errors  in  the  sensors  and  instrumenta¬ 
tion  (including  possible  observer  error  in  reading  the  memoscope)  indicate  an 
overall  error  of  +20. 13  percent  in  the  observed  data. 

The  differences  in  observed  and  computed  mean  heave  acceleration  ampli¬ 
tudes  fall  within  these  error  bands  and  thus  indicate  that  the  computer  model 
offers  reasonable  accuracy  for  this  motion.  However,  even  if  these  error  bands 
are  taken  into  account,  the  model  is  systematically  underestimating  surge  and 
pitch  motions.  Since  surge  and  sway  motions  are  underestimated  with  both  the 
spherical  buoy  and  with  the  torroidal  buoy,  it  is  suspected  that  the  transverse 


hydrodynamic  mass  and  damping  used  in  the  model  may  be  in  error.  Recent 


*Note  that  wave  heights  for  the  observed  data  were  based  upon  visual 
observations  while  wave  heights  for  the  simulated  buoy  motions  were  computed 
by  the  S-M-B  method. 
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communications  with  E.  Geller  and  R.  Canada,  who  are  affiliated  with  the 
National  Data  Buoy  Project,  have  indicated  that  the  results  of  their  model  buoy 
tests  in  towing  basins  do  not  agree  with  present  theory  because  surge  and  sway 
hydrodynamic  forces  are  less  than  predicted  analytically.  This  communication 
confirms  the  above  suspicion,  but  judgment  is  reserved  until  more  test  tank 
data  are  published. 

The  environmental  conditions  as  measured  during  the  data  run  of  11  June 

1970  at  2030  EDST  were  used  as  input  to  the  computer  model,  and  buoy  motions 

were  computed.  Power  spectra  for  wave  height,  buoy  heave  accelerations, 

buoy  sway  acceleration,  buoy  surge  acceleration,  buoy  pitch  motions,  and  buoy 

84—86 

roll  motions  were  computed  by  a  Fast  Fourier  Transform  method.  A  total 

of  1,024  samples,  sampled  at  0.06-sec  intervals  for  each  parameter,  were 

transformed.  The  samples  were  smoothed  by  averaging  over  a  0.3-sec  interval 

to  prevent  aliasing  in  the  spectra.  Eight  ensemble  averages  were  used,  and  they 

resulted  in  64  statistical  degrees  of  freedom  and  a  standard  error  of  17.66  per- 

87 

cent.  The  frequency  resolution  is  0. 0163  Hz.  Williams  states  that  a 
resolution  on  the  order  of  0.02  IIz  is  adequate  to  define  ocean  swell  spectral 
peaks  in  his  study  of  ocean  wave  spectra  in  Block  Island  Sound. 

The  input  wave  spectrum  is  shown  in  figure  S2  along  with  the  frequencies 
of  the  components  that  form  the  "random"  wave  forcing  function.  Figure  83 
shows  the  computed  spectrum  level  of  the  "random"  wave  model  compared  with 
the  theoretical  spectrum  level.  As  expected,  the  "random"  wave  model  exhibits 
characteristics  of  narrow-band  white  noise,  which  is  reasonable  since  it  is 
made  up  of  a  finite  number  of  sine  wave  components  close  to  one  another  in 
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Figure  83,  Simulated  Wave  Height  Spectral  Levels 
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the  frequency  domain. 

The  standard  deviations  for  the  input  wave  height  and  the  output  buoy 
motions  are  as  follows: 


Wave  height 

0.4807  ft 

Heave  acceleration 

1.9321  ft/sec' 

Sway  acceleration 

0.4217  ft/sec 

Surge  acceleration 

0.4974  ft/sec 

Pitch  angle 

2. 98  deg 

Roll  angle 

4.437  deg. 

Output  spectral  levels  are  shown  in  figures  84  through  88.  In  general,  the 
spectra  exhibit  more  deterministic  properties  than  the  buoy  motion  spectra 
from  the  at-sea  data.  More  smoothing  was  done  on  the  at-sea  spectra  than  on 
the  simulated  spectra,  but  the  at-sea  spectra  do  not  appear  to  have  the  narrow- 
band  characteristics  of  the  simulated  spectra,  especially  in  heave,  surge,  and 
pitch.  This  indicates  that  the  model  may  be  inadequately  damped.  The  natural 
frequency  of  the  buoy  in  heave  when  at  middraft  is  computed  to  be  0.36  Hz.  The 
heave,  sui'ge,  and  pitch  spectra  exhibit  peaks  at  about  this  frequency.  Since 
surge  and  pitch  are  decoupled  from  heave  hydrodynamically,  the  coupling  must 
be  effected  through  the  cable  tensions  acting  on  the  buoy. 

The  second  peak  in  the  heave  spectrum  is  located  at  about  the  same 
frequency  as  the  peak  of  the  water  particle  acceleration  spectra.  The  forcing 
function  component  due  to  the  heave  hydrodynamic  inertia  force  caused  by  water 
particle  acceleration  is  driving  the  buoy  in  this  frequency  range.  The  surge 
spectrum  also  indicates  that  this  forcing  mechanism  is  active  in  that  mode. 
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FREQUENCY  (Hi) 


Figure  88.  Simulated  Roll  Spectral  Levels 


Both  pitch  and  roll  spectra  show  peaks  at  this  frequency.  Pitch  and  roll 
motions  have  computed  natural  frequencies  at  about  0.7  Hz  and  would  tend  to 
respond  to  this  forcing  mechanism  through  coupling  from  surge  and  sway 
motions.  Since  the  water  particle  motions  are  constrained  to  the  x  -  z  plane, 
one  would  not  expect  a  response  out  of  plane.  However,  the  current  vector  has 
an  out-of-plane  component;  thus,  the  mooring  line  tension  force  will  have  an  out- 
of-plane  component.  Heave  and  surge  motions  will  couple  into  sway  and  roll 
motions  because  of  the  mooring  line  tension.  This  coupling  also  illustrates 
another  unusual  feature  of  this  dynamic  system  in  that  the  response  of  the 
system  is  strongly  dependent  on  the  mean  values  of  the  system  element  spatial 
locations  —  a  strong  argument  against  linearized,  decoupled  models  of  buoy 
systems. 

Both  heave  and  surge  show  peaks  at  higher  frequencies  (1.05  and  1.4  Hz). 
These  peaks  are  probably  due  to  natural  frequencies  for  the  lumped-mass  cable 
model.  Sway  and  roll  spectra  have  nearly  the  same  shape  and  thus  indicate  a 
strong  coupling  in  these  modes.  Surge  and  sway  (also,  pitch  and  roll  through 
coupling)  indicate  some  response  at  very  low  frequencies.  This  is  due  to  the 
lowest  natural  frequency  of  the  buoy  system  in  a  horizontal  mode.  If  the 
analogy  is  made  with  a  pendulum  having  a  length  order  of  magnitude  with  the 
water  depth,  the  natural  frequency  of  this  motion  would  be  quite  low. 

The  heave  spectrum  computed  from  the  experimental  data  (figure  57) 
indicates  a  peak  at  0.4  Hz,  which  agrees  well  with  the  peak  in  the  simulated 
heave  spectrum  due  to  the  buoys  natural  frequency  in  heave.  Surge,  sway, 
pitch,  and  roll  spectra  all  have  spikes  at  0.6  Hz.  This  is  in  agreement  with 


the  peaks  in  the  simulated  sway,  surge,  pitch,  and  roll  spectra.  All  the 
spectra,  especially  in  roll,  indicate  energy  at  frequencies  out  to  2  Hz.  This 
energy  may  be  due  to  vibration  in  the  mooring  line  coupling  into  the  buoy. 

This  comparison  of  simulated  buoy  motions  with  observed  buoy  motions 
for  the  torroid  buoy  at  station  BRAVO  on  11  June  1970  indicates  that  the  com¬ 
puter  model  is  distributing  energy  in  its  response  in  about  the  same  way  as  the 
actual  buoy  system.  However,  the  model  is  apparently  underdamped  and  is 
filtering  out  some  energy  between  natural  frequencies.  Again,  more  test  tank 
data  are  needed  for  buoy  hulls  to  determine  their  hydrodynamic  characteristics. 

The  simulations  of  buoy  system  dynamics  for  the  spherical  and  torroidal 
buoys  moored  in  Block  Island  Sound  have  used  the  lumped- mass  model  of  cable 
dynamics.  In  the  course  of  this  research,  it  was  found  that  if  mean  tensions  in 
the  cable  are  very  low  or  if  a  number  of  force  discontinuities  are  present 
along  the  cable  (both  conditions  common  to  shallow  water  moorings),  the  finite- 
difference  method  is  usually  not  suitable.  With  very  low  tensions,  the  cable 
equations  can  go  ultrahyperbolic,  which  could  cause  the  numerical  method  to  break 
down.  When  many  force  discontinuities  are  in  the  line,  a  large  number  of  nodes 
are  needed  in  the  cable  segments  between  the  discontinuities  in  addition  to  equation 
of  motion  for  each  discontinuity.  The  resulting  computational  time  becomes 
prohibitive. 

4.3.2  WIIQI  Mooring  No.  238 

In  order  to  validate  the  buoy  dynamics  simulatation  using  the  finite- 
difference  cable  model,  mooring  line  tension  data  taken  with  WHOI  mooring 
No.  238  (reference  37)  is  compared  with  simulated  tension  data  for  the  same 
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buoy  system.  WHOI  mooring  No.  238  is  essentially  the  same  as  WHOI  mooring 
No.  279  (figure  27)  and  was  moored  at  the  same  location.  Tensions  were  meas¬ 
ured  just  below  the  torroidal  surface  buoy  and  were  telemetered  ashore.  Com¬ 
puted  inputs  for  the  simulation  of  WHOI  mooring  No.  279  are  shown  in  appendix  < 
D,  and  the  buoy  system  dynamics  program  incorporating  the  finite-difference 
cable  model  is  shown  in  appendix  B.  Since  the  WHOI  data  are  shown  against 
wind  speed,  the  primary  input  is  wind  speed.  The  fetch  was  assumed  to  be 
100  miles  for  this  location  in  the  North  Atlantic,  and  wind  durations  were  taken 
to  be  24  hr.  The  Webster  current  profile,  with  a  surface  current  of  1.5  ft/sec, 
was  assumed  to  act  in  the  direction  of  the  wind  for  all  cases. 

Strains  and  cable  angles  at  each  node  must  be  read  in  as  initial  conditions. 

The  steady-state  buoy  system  configuration  program  was  used  to  compute  these 

parameters  for  each  wind  speed  and  the  given  current  profile.  The  initial  cable 

angles  are  shown  in  figures  89  and  90.  The  standard  deviation  for  the  cable 

37 

tension  just  below  the  buoy  was  computed  for  each  run.  Millard  presents  a 
"scattergram"  of  "dynamic  tension  amplitude"  versus  2-hr  mean  wind  speed. 
Tensions  wore  recorded  on  a  Rustrak  recorder,  which  has  a  very  slow  chart 
speed  —  on  the  order  of  centimeters  per  day.  Thus,  this  "amplitude"  was  read 
by  measuring  the  breadth  of  a  very  thick  line.  This  actually  represents  the 
tension  difference  between  the  highest  and  lowest  tension  that  occurred  in  a  2-hr 
period,  it  is  assumed  that  these  "amplitudes"  correspond  to  the  highest  l/10th 
wave  heights  found  in  wind- wave  height  distributions.  The  WHOI  data  of 
"dynamic  tension  amplitude"  versus  wind  speed  is  shown  in  figure  91. 
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Figure  89.  Initial  Strains  for  WHOI  Mooring  No.  238 
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Figure  90.  Initial  Angles 


for  WHOI  Mooring  No.  238 


245 


4  r 


Tension  standard  deviations  for  each  run  were  first  converted  to  mean 
amplitudes  by  assuming  sinusoidal  variations,  converted  to  highest  l/10th 
amplitudes  by  assuming  a  Rayleigh  distribution,  and  finally  converted  to  double 
amplitudes.  The  resulting  expression  is 

Twm~(l.S7)(z.o3)(a.o)-Tso  . 

The  simulated  "dynamic  tension  amplitudes"  are  plotted  in  figure  91. 
Errors  for  the  simulation  were  not  computed  in  this  case  since  the  meaning  of 
the  W.HOI  "dynamic  tension  amplitude"  is  not  clear.  However,  the  simulation 
computes  dynamic  tensions  that  are  order  of  magnitude  and  that  increase  with 
increasing  wind  speed.  A  better  set  of  buoy  motion  and  tension  measurements 
are  needed  for  deep  sea  buoy  systems  in  order  to  fully  validate  the  buoy  motion 


dynamic  simulation  with  the  finite-difference  cable  model. 


V.  SUMMARY 


5.1  Restatement  of  the  Problem 

The  object  of  this  investigation  was  to  evaluate  the  forces  acting  on  the 
components  of  a  simple  buoy  system  exposed  to  the  oceanic  environment  and  to 
develop  a  numerical  model  of  buoy  system  dynamics  to  simulate  buoy  system 
response.  Axisymmetrie  buoy  hull  shapes  were  considered  in  general,  and 
hydrostatic  and  hydrodynamic  forces  and  moments  on  oblate  spheroids  were 
studied.  The  set  of  integro-differential  equations  of  motion  for  the  buoy  were 
reduced  to  a  set  oi  ordinary  differential  equations  with  nonconstant  coefficients 
by  using  frequency  dependent  hydrodynamic  force  coefficients  published  in  the 
literature. 

A  "quasi- random  "wind  wave  model  was  developed  to  simulate  the  motions 
ot  the  water  masses  in  the  immediate  vicinity  of  the  buoy.  Wind  wave  properties 
were  computed  with  the  S-M-B  method  from  the  mean  wind  speed,  fetch,  and 
duration.  Borgman's  energy  partitioning  method  was  applied  to  a  two-component 
Bretschneider  spectrum,  and  sinusoidal  wave  component  amplitudes  and  fre¬ 
quencies  were  computed.  The  random  phase  components  were  summed  to 
compute  instantaneous  water  particle  motions. 

Major  assumptions  made  in  the  investigation  of  buoy  dynamics  include  the 
following: 
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1.  Dissipative  forces  can  ho  separated  into  those  due  to  surface  wave 
generation  and  those  due  to  viscous  drag. 

2.  Infinitesimal  buoy  motions  exist;  this  assumption  was  made  by 
22  23 

John  ’  in  the  derivation  of  the  hydrodynamic  force  and  moment  integrals 

42 

and  subsequently  by  Kim,  who  evaluated  the  integrals  for  oblate  sphenoids. 

3.  The  Ilaskind  hypothesis  is  valid;  i.e. ,  it  is  assumed  that  the  general 
problem  of  an  object  moving  on  a  free  surface  in  response  to  gravity  waves  on 
that  free  surface  can  be  linearized  to  the  extent  that  the  velocity  potential  for 
that  motion  is  the  sum  of 

a.  the  linearized  velocity  potential  of  the  gravity  waves  alone, 

b.  the  velocity  potential  of  the  object  moving  on  the  free  surface 
with  no  waves,  and 

c.  the  velocity  potential  due  to  the  waves  generated  by  the  motion 

of  the  body. 

The  Ilaskind  hypothesis  was  used  in  order  to  apply  Kim's  coefficients  for  oblate 
spheroids,  which  were  derived  for  objects  oscillating  on  a  free  surface,  to  the 
case  of  objects  oscillating  on  a  free  surface  with  gravity  waves. 

4.  The  St.  Denis- Pierson  hypothesis  is  valid,  i.e.,  it  was  assumed  that 
the  sea  can  be  represented  as  the  linear  sum  of  elementary  waves  of  random 
phase. 

5.  Body  dimensions  are  small  compared  with  a  wavelength.  This  assump¬ 
tion  was  made  in  orde*  to  use  the  computed  water  particle  motions  as  the 
motions  of  the  mass  of  water  in  the  immediate  vicinity  of  the  buoy  hull. 
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Equations  of  motion  for  an  elastic  cable  capable  of  supporting  b'-'h  axial 
and  transverse  waves  simultaneously  were  written.  A  unique  finite-difference 
numerical  technique,  an  extension  of  Hartree's  method  for  hyperbolic  partial 
differential  equations,  was  developed  for  the  solution  of  sets  of  coupled  hyper¬ 
bolic  equations.  A  lemped-mass  cable  model  was  also  developed.  Major  assump¬ 
tions  made  in  the  development  of  the  cable  dynamics  model  are  that  the  cable  is 
homogeneous  and  perfectly  flexible  and  that  hysterisis  damping  is  negligible 
compared  with  viscous  damping.  The  equations  of  motion  for  the  buoy  and  cable 
were  programmed  and  solved  numerically  on  a  UNIVAC  1108  digital  computer  in 
the  time  domain.  Two  types  of  buoy  hulls  were  considered  —  a  spherical  buoy 
and  a  torroidal  buoy.  Both  shallow  water  and  deep  water  moorings  were  simu¬ 
lated  using  the  lumped- mass  and  finite-difference  cable  models.  Simulated  data 
were  compared  with  observed  data  in  two  steady-state  cases  and  in  three 
dynamics  cases. 

5 . 2  Conclusions 

Steady-state  buoy  system  configurations  were  simulated  using  the  method 
described  in  chapter  III  for  a  shallow  water  buov  system  and  for  a  deep  water 
buoy  system.  Comparison  of  simulated  tc  observed  configuration  parameters 
(current  meter  tilt  angles  for  the  shallow  water  mooring  and  mean  tensions  for 
the  deep  water  mooring)  indicate  good  agreement  bctw-cen  the  computer  model 
configurations  and  the  configurations  of  the  actual  buoy  systems.  Although  this 
comparison  of  a  few  bits  of  data  from  two  buoy  systems  does  not  constitute  a 
full  validation  of  the  steady-state  buoy  system  configuration  model  (angles  and 
tensions  all  along  the  mooring  line  for  many  current  profiles  and  winds  and  for 
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many  different  buoy  systems  should  he  simulated  and  correlated  for  a  full 
validation),  it  does  indicate  that  the  simulation  errors  are  on  the  order  of 
5  percent,  which  is  adequate  for  engineering  analysis. 

Buoy  system  dynamics  were  simulated  for  three  cases  (two  shallow  water 
moorings  where  the  lumped-mass  cal;*,  nodel  was  used  and  one  deep  water  moor¬ 
ing  where  the  finite-difference  cable  model  was  used),  and  simulated  buoy  motion 
parameters  were  compared  with  observed  parameters.  In  general,  (he  lumped- 
mass  cable  model  is  more  adaptable  to  shallow  water  buoy  systems,  which,  if 
slack-moored,  tend  to  have  very  low  tensions  at  the  bottom  and  cable  angles  that 
approach  T172  and,  if  laut-moored,  tend  to  have  very  high  dynamic  strain  levels, 
which  can  cause  slack  cable.  The  lumped-mass  cable  model  can  handle  these 
cases  easily  (if  tensions  on  an  element  go  to  zero,  the  element  simply  free  falls 
through  the  water),  whereas  the  finite- difference  cable  model  breaks  down  with 
slack  cable  conditions  since  the  equations  of  motion  are  ultrahyperbolic  and 
possess  an  infinite  numbev  of  equally  valid  solutions.  Also,  the  lumped-mass 
model  is  more  adaptable  when  the  mooring  line  eonlains  a  number  of  mass  or 
force  discontinuities,  i.c.,  instruments,  subsurface  buoys,  sentinels,  etc. 

The  comparison  for  both  shallow  water  cases,  one  a  spherical  buoy  and 
the  other  a  Richardson  lorroid,  indicate  that  the  surge  and  sway  hydrodynamic 
forces  were  underestimated  in  the  simulation.  Also,  the  comparison  of  observed 
to  sin  dated  shallow  water  torroidal  buov  motion  spectral  levels  indict  tes  that 
the  simulation  is  underdamped  in  all  modes.  There  is  good  agreement  in 
simulated  and  observed  heave  motions,  the  most  important  motion  parameter 


for  the  buoy  system  designer.  In  general,  the  simulation  appears  to  predict 
heave  motions  within +15  percent  and  the  other  motions  within  +50  percent. 

Since  the  environmental  conditions  were  not  monitored  at  either  buoy  during  the 
motion  measurements  but  were  inferred  from  wind  speeds  measured  ashore, 
from  computed  tidal  currents  based  upon  previous  current  measurements,  and 
from  visual  observations  of  sea  conditions,  it  is  impossible  to  draw  conclusions 
on  the  validity  of  the  model,  except  that  it  computes  buoy  system  motions  that  are 
order  of  magnitude  with  observed  motions. 

The  comparison  of  simulated  deep  sea  buoy  system  dynamics  using  the 
finite-difference  model  is  inconclusive  because  of  uncertainties  in  monitoring 
the  environment  and  in  the  statistical  meaning  of  the  tension  data  collected. 

5 . .?  Suggestions  for  Further  Study 

Future  research  in  the  area  of  buoy  system  dynamics  should  involve  the 
omplanlment  and  the  fitting  of  motion  sensing  instrumentation  to  a  wide  spectrum 
of  oceanic  buoy  system  types  —  both  shallow  and  deep  water.  The  environment 
at  each  site  should  be  adequately  monitored  (winds,  waves,  and  currents),  and 
a  complete  set  of  buoy  system  motions  should  be  recorded.  The  measurements 
should  include  angles  and  tensions  along  the  mooring  line  as  well. 

Another  key  area  that  should  be  investigated  involves  the  hydrodynamic 
forces  acting  on  a  body  on  the  free  surface  of  a  fluid  when  the  free  surface  is 
subjected  to  random  gravity  waves.  John's  analysis  assumed  infinitesimal  body 
motions  in  rrder  to  linearize  the  free  surface  boundary  conditions.  Kim's 
analysis  assumed  sinusoidal  body  motions  in  order  to  evaluate  the  hydrodynamic 
fer’e  and  moment  integrals.  An  investigation  of  the  validity  of  these  assumptions 
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should  be  made  in  order  to  better  understand  the  nature  of  the  body-fluid  inter¬ 
action.  Also,  test  tank  data  are  needed  to  validate  theoretical  force  and 
moment  coefficients. 

Extension  of  the  analysis  presented  in  this  dissertation  to  nonaxisymmet- 
ric  buoy  hulls  to  include  cross-coupled  hydrodynamic  forces  and  moments  in 
the  other  modes  would  be  of  significance  as  a  more  general  study  of  buoy 
system  dynamics. 

Finally,  more  research  on  the  behavior  of  wire  ropes  and  synthetic  lines 
that  are  used  in  the  ocean  environment  is  needed  for  the  prediction  of  buoy 


system  performance. 
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Appendix  A 


ANGULAR  STABILITY  OF  AXISYMMETRIC  BUOYS 

It  is  known  from  visual  observation  and  from  collected  motion  data  that 
buoys  do  not  become  unstable  in  the  sense  of  Liapunov;  i.e. ,  buoy  motions  will 
not  build  to  infinity.  However,  buoys  will  undergo  large  excursions  and  can  be 
upset.  For  example,  a  buoy  that  undergoes  gross  heaving  motions  in  ihe  sea  so 
that  it  is  alternately  awash  and  then  rises  out  of  water  to  fall  over  on  its  side 
would  be  deemed  unstable  by  the  casual  observer;  however,  the  buoy  would 
actually  be  stable  since  its  motions  are  bounded  and  do  not  tend  to  infinity  with 
time.  Of  course,  if  this  buoy  housed  meteorologic  or  oceanographic  instru¬ 
mentation,  it  would  be  of  little  value  because  of  its  wild  motions. 

The  "stability"  of  this  type  of  motion  is  best  described  in  the  pitch  (or 
roll)  phase  plane,  i.e. ,  a  plot  of  pitch  angular  velocity  versus  pitch.  However, 
pitching  (and  rolling)  motions  are  heavily  influenced  by  the  draft  of  the  buoy, 
i.e. ,  the  location  of  the  center  of  buoyancy  relative  to  the  center  of  gravity. 
Thus,  the  heave  motions  and  resulting  buoy  draft  must  be  considered.  Because 
the  righting  moment  depends  on  the  location  of  the  center  of  buoyancy  relative 
to  the  center  of  gravity,  the  pitch  (or  roll)  equations  ox  motion  are  similar  to 
the  nonlinear  equations  of  motion  for  a  pendulum. 

A  wide  spectrum  of  axisymmetric  buoy  shapes  will  be  considered,  ranging 
from  a  spar  buoy  with  a  high  draft-to-beam  ratio  to  a  discus  buoy  with  a  low 
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draft-to-beam  ratio  (figure  A-l).  First,  consider  the  spar  buoy.  In  the  absence 
of  external  forces  (mooring  line  tensions,  etc.),  static  stability  is  maintained 
only  if  the  center  of  gravity  is  below  the  center  of  buoyancy: 

Lcq  <  W/z.  . 

This  is  obvious  from  the  undamped  pitch  equation  of  motion  for  the  spar  buoy: 

Ijt  J+S§tFB*H(£- Lcer )  smfi  =  o  _ 

Also,  if  the  buoy  is  tilted  at  extreme  angles,  we  see  that 

jS  >TT  or  JQ  <  -TT  t 

The  spar  buoy  will  not  return  to  its  initial  equilibrium  position  but  will  undergo 
a  complete  resolution  and  seek  equilibrium  at  either  27T  or  -27 T.  A  typical 
pitch  phase  plane  for  the  spar  buoy  is  shown  on  figure  A-2  for  two  possible 
conditions; 

1*  J-'Cfr  ^  H/e 

2-  l~CGr  >  W/z.  , 

Note  that  if  the  buoy  is  initially  unstable,  L  >  H/2  ,  it  is  stable  at  +  TT; 
i.e. ,  it  is  stable  upside  down. 

A  cable  attached  to  the  bottom  of  the  spar  buoy  will  tend  to  stabilize  the 
buoy.  Adding  a  restoring  moment  due  to  the  cable  tension  in  the  equation  of 
motion  yields 

Ifij$  +^5  H  (h/2  “  ^-CG-)  SINJ3 

+  K  Lee,  (1-cosJB)  lC(k  SiN^  —  O  ) 

where  K  is  a  cable  spring  constant.  For  a  stable  system,  the  sum  of  the  tw'o 

restoring  moment  terms  must  be  positive.  Solving  for  L  ,  we  see  that 

00 


Figure  A-l.  Axisymmetrie  Buoy  Shapes 


SPAR  BUOY  PITCH 
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L__  >  H/2 
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u  =  f^/c’-yc^r , 

where 

j>girB*H 

^  ~  t+Ktt-cosjs)  . 

In  general,  a  spar  buoy  is  designed  to  minimize  heave  response.  Thus, 
the  shallowest  buoy  draft  would  be  the  static  draft  of  the  buoy  minus  one- half 
the  largest  wave  height: 


^mih  H  o  g  \  MAX 

For  a  stable  buoy,  we  see  that 


Hmih  _  H0-  z'K  MAa 


Lcs  <  ~  2 

The  spherical  buoy  shape  is  subject  to  the  same  types  of  moments,  but 

the  study  is  complicated  by  the  fact  that  the  center  of  buoyancy  will  deviate  from 

the  vertical  axis  of  the  buoy  as  the  buoy  pitches.  If  the  center  of  gravity  lies 

below  the  center  of  buoyancy  (5/8  II  from  the  bottom),  the  buoy  is  positively 

stable  and  will  have  a  phase  plane  representation  similar  to  that  shown  in 

figure  A-2  for  the  spar  buoy.  The  upper  phase  plane  applies  if  L  <5/8  H, 

and  the  lower  phase  plane  applies  if  L  >  H  .  However,  the  buoy  is  neutrally 

LG 

stable  if  5/8  H  <  L  <  H  ,  since  the  center  of  buoyancy  is  always  directly 
below  the  center  of  gravity.  In  this  situation  there  is  no  definite  stable  position, 
because  there  is  no  restoring  moment.  Again,  the  tension  of  a  mooring  line 
attached  below  the  center  of  buoyancy  will  stabilize  the  buoy. 

Any  buoy  with  a  draft- to- beam  ratio  less  than  0. 5  and  a  hull  height-to- 
beam  ratio  less  than  1  can  be  bistable;  i.e. ,  it  can  be  stable  right  side  up  or 


upside  down.  Torroid  and  discus  buoys  are  of  this  class.  If  the  center  of 
gravity  is  located  far  below  the  buoy  hull  (achieved  by  means  of  ballast  weights 
on  a  boom  or  tripod),  the  phase  plane  is  similar  to  the  upper  curve  for  the  spar 
buoy.  If  the  center  of  gravity  is  far  above  the  buoy  hull  (due  to  heavy  instru¬ 
ments  or  equipment),  the  lower  curve  for  the  spar  buoy  phase  plane  would 
apply  —  the  buoy  is  stable  when  upside  down.  If  the  center  of  gravity  is  near 
the  geometric  center  of  the  buoy  hull,  the  buoy  is  equally  stable  right  side  up  or 
upside  down.  The  phase  plane  for  this  situation  is  shown  on  figure  A-3.  Inspec¬ 
tion  of  figure  A-3  indicates  that  the  width  of  the  stable  (in  the  sense  that  the 
buoy  is  right-side-up)  region  in  the  phase  plane  can  vary  from  21Tto  0  depend¬ 
ing  on  the  vertical  location  of  the  center  of  gravity. 

This  simple  discussion  of  buoy  stability  did  not  consider  other  modes  of 
possible  unstable  motion  since  they  have  never  been  observed  to  offer  serious 
problems.  Cross-coupled  moments  due  to  the  hydrodynamic  forces  acting  on 
the  buoy  and  the  horizontal  tension  components  were  neglected  in  this  discussion. 
Even  with  those  restrictions,  a  few  design  guides  are  apparent.  The  buoy  should 
be  designed  with  the  center  of  gravity  below  the  center  of  buoyancy  if  at  all 
possible.  Also,  the  mooring  line  attachment  point  should  be  as  low  as  possible 
to  offer  the  greatest  righting  moment  if  the  buoy  does  capsize. 


GEOMETRIC  CENTER 


PITCH  PHASE  PLANES  FOR 
TORROID  OR  DISCUS  BUOYS 


Appendix  B 


COMPUTER  PROGRAMS  FOR  THE  MODELS 

Steady- State  Buoy  System  Configurations 

This  program  computes  the  three-dimensional  spatial  configuration, 
tensions,  and  strains  of  either  an  elastic  or  inelastic  buoy  system  mooring 
line.  Winds  from  any  compass  direction  can  act  on  the  buoy,  and  currents  that 
vary  in  strength  and  direction  as  a  function  of  depth  can  act  on  the  bucy  and 
mooring  line.  The  mooring  line  can  be  composed  of  segments  having  different 
properties  (weight  in  water,  mass,  drag,  elasticity,  etc.).  Point  mass  discon¬ 
tinuities  (to  simulate  current  meters,  hydrophones,  etc.)  can  also  be  accounted 
for. 

The  basic  cable  equations  (equations  (230)  to  (233)) and  a  discussion  of  the 
development  of  this  program  are  included  in  chapter  III  of  the  main  text.  The 
logic  employed  in  this  program  is  shown  in  figure  B-l,  which  generally  illus¬ 
trate:;  the  computational  operations. 

>rhe  input  data  are  as  follows: 

; .  Buoy  major  diameter  (ft) 

2.  Buoy  minor  diameter  (ft) 

a.  Vertical  diameter  for  an  oblate  spheroid 

b.  Hole  diameter  for  a  torroid 
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Figure  B-l.  Steady-State  Buoy  System  Configuration  Flow  Chart 
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c.  0  for  a  cylinder 


d.  0  for  a  sphere 

2.  Buoy  weight  (lb) 

4.  Maximum  draft  of  hull  (ft) 

5.  Free  draft  with  no  mooring  line  (ft) 

2 

6.  Buoy  windage  (ft  ) 

7.  Wind  drag  coefficient 

8.  Cable  diameter  (in.) 

9.  Cable  weight  in  water  per  unit  length  (lb/ft) 

.  2 

10.  Effective  cable  modulus  of  elasticity  (lb/in.  ) 

11.  Unstretched  cable  length  (ft) 

12.  Surface  current  (knots) 

13.  Water  depth  (ft) 

14.  Wind  speed  component  in  y  direction  (ft/sec) 

15.  Wind  speed  components  in  z  direction  (ft/sec). 

The  integration  step  size  "B"  is  normally  set  at  l/100th  of  the  total  cable 
length.  However,  if  the  cable  properties  are  changed  from  a  lightweight  line  to 
a  very  heavy  line  (for  example,  anchor  chain),  the  step  size  should  be  changed 
in  inverse  proportion  to  the  in-water  weights.  Also,  the  normal  and  tangential 
drag  coefficients  ("DRGON”  and  "DRGC'f")  should  be  changed  accordingly  if  the 
mooring  line  section  is  not  circular,  for  example,  hair-faired  cable  and  chain. 
Occasionally,  the  solution  will  not  converge  into  the  depth  error  band.  This 


occurs  if  the  first  draft  increment  is  too  large  and  the  computed  x  dimension  of 


the  mooring  line  "overshoots"  the  depth.  With  each  overshoot,  the  buoy  draft 
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increment  is  halved;  thus,  convergence  may  be  very  slow.  The  initial  buoy 
draft  increment  should  be  halved  if  this  occurs . 
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Buoy  System  Dynamics  for  the 
Spherical  Buoy 

This  program  simulates  the  motions  of  a  spherical  buoy  exposed  to  winds, 
currents,  and  a  random  sea.  The  mooring  line  is  simulated  as  five,  elastically 
connected,  lumped  masses.  Cable  weights,  drag  forces,  and  inertia  forces  are 
concentrated  at  each  mass.  The  buoy  is  allowed  six  degrees  of  freedom,  and 
each  mooring  line  mass  element  is  allowed  three  translational  degrees  of  free¬ 
dom.  With  a  total  of  21  degrees  of  freedom,  42  first-order  differential  equations 
are  integrated  simultaneously  in  the  time  domain. 

The  equations  of  motion  for  the  buoy  and  the  development  of  the  forces 
acting  on  the  buoy  are  presented  in  chapter  III  of  the  main  text.  The  lumped- 
mass  cable  equations  are  also  shown  in  that  chapter. 

Major  computational  procedures  in  this  program  are  shown  in  figure  B-2. 

The  input  data  are  as  follows : 

1.  Buoy  hull  radius  (ft) 

2.  Height  of  the  center  of  gravity  of  the  buoy  above  the  mooring  line 
connection  point  on  the  buoy  (ft) 

3.  Height  of  the  mooring  line  connection  point  below  the  buoy  hull  (ft) 

4.  Buoy  weight  (lb) 

2 

5.  Buoy  structural  and  floodwater  mass  (lb-sec  /ft) 

2 

6.  Yaw  mass  moment  of  inertia  (lb-sec  /ft) 

2 

7.  Pitch  mass  moment  of  inertia  (lb-sec  /ft) 

2 

8.  Roll  mass  moment  of  inertia  (lb-sec  /ft) 

9.  Effective  buoy  wind  drag  coefficient  (dimensionless) 
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10.  Effective  buoy  wind  lift  cc  efficient  (dimensionless) 

2 

11.  Buoy  windage  (profile  area)  (ft  ) 

2 

12.  Buoy  plan  area  (ft  ) 

13.  Height  of  the  wind  center  of  pressure  above  the  center  of  gravity  of 
the  buoy  (ft) 

14.  Mean  wave  height  (ft) 

15.  Mean  wave  period  (sec) 

16.  The  unstretched  cable  lengths  between  mass  elements  (ft)  (6  required) 

17.  Upper  mooring  line  segment  diameter  (ft) 

18.  Upper  mooring  line  segment  weight  in  water  per  unit  length  (Ib/ft) 

2 

19.  Upper  mooring  line  segment  mass  per  unit  length  (lb-sec  /ft) 

20.  Lower  mooring  line  segment  diameter  (ft) 

21.  Lower  mooring  line  segment  weight  in  water  per  unit  length  (lb/ft) 

2  . 

22.  Lower  mooring  line  segment  mass  per  unit  length  (lb-sec  /ft) 

23.  Surface  current  y  component  (ft/sec) 

24.  Surface  current  z  component  (ft/sec) 

25.  Wind  speed  y  component  (ft/sec) 

26.  Wind  speed  z  component  (ft/sec) 

27.  Initial  buoy  displacements  x,  y,  and  z  (ft) 

28.  Initial  cable  element  displacements  (including  the  anchor  clump) 
x,  y,  and  z  for  each  of  6  sets  (ft). 

The  user  of  this  program  should  first  estimate  the  highest  natural  fre¬ 
quency  in  the  system.  In  general,  the  upper  mooring  line  segment  will  be 
lightest  and  the  highest  natural  frequencies  are  in  the  axial  mode  along  the 
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cable.  Using  the  same  variable  names  as  those  in  the  program  we  can  estimate 
the  highest  natural  frequency: 

v?  _  i  iz . *  ea.  it 

'hi  ~  ZTT^j  OCSMi-CiOCZ)  ) 

where 

EA  is  the  cable  modulus*  (lb) 

DCSM1  is  the  mass  per  unit  length  of  the  upper  mooring  line  segment 
(lb-sec2/ft2) 

CLO(2)  is  the  unstretched  cable  length  between  the  first  and  second 
cable  mass  elements  (ft). 

It  the  .:able  lengths  were  very  short  or  the  cable  modulus  were  very  high  in  the 
lower  mooring  line  segment,  an  estimate  should  be  made  of  its  highest  natural 
frequency  also. 

For  numerical  stability,  the  integration  step  size  should  be  about  l/2Cth 
the  shortest  period  present.  Thus,  the  step  size  is  computed: 

R  ~o.o?( l/fj  . 

In  this  program,  the  wave  component  amplitudes,  frequencies,  and  phase 
angles  are  computed  externally  and  are  listed  in  the  body  of  the  program  (MOMEG, 
AMP,  and  PHS).  The  water  depth  (DEEP)  is  also  listed,  and  the  x  coordinate 
of  the  mooring  anchor  elump  should  be  set  equal  to  the  water  depth.  Tb*  total 
time  of  the  simulation  is  controlled  by  a  logical  "IF"  statement  (statement 


+The  cable  modulus  is  the  product  of  the  cable  material  elastic  modulus 
and  the  actual  cable  cross-sectional  area. 
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'4jP  no-  491),  which  shifts  to  the  "STOP"  control  if  the  maximum  simulated  time  is 

exceeded. 

In  this  program,  buoy  motion  output  data  are  not  printed  for  the  first 
6  sec  of  simulation  as  initial  transients  decay.  After  this  time,  the  following 
buoy  system  outputs  are  printed  every  800  time  steps: 

Simulation  time  (sec) 

2 

Water  particle  vertical  acceleration  (ft/sec  ) 

2 

Buoy  heave  acceleration  (ft/sec  ) 

2 

Buoy  sway  acceleration  (ft/sec  ) 

2 

Buoy  surge  acceleration  (ft/sec  ) 

Buoy  pitch  angle  (deg) 

Buoy  roll  angle  (deg) . 
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Buoy  System  Dynamics  for  the  Torroidal 
Buoy  at  Station  BRAVO 

The  previous  program,  developed  for  the  spherical  buoy,  was  modified  to 
simulate  the  torroidal  buoy  dynamics.  As  indicated  in  figure  B-3,  the  basic 
computational  procedures  remain  the  same  as  for  the  spherical  buoy.  However, 
the  wave  component  amplitudes,  frequencies,  and  phase  angles  are  computed  in 
a  subroutine  (RWAVE)  in  the  program  and  do  not  have  to  be  listed  in  the  program. 
The  buoyant  forces  and  moments  for  the  torroidal  buoy  are  also  computed  in  a 
subroutine  (TORBU)  using  an  integration  method  developed  in  appendix  E. 

Finally,  since  the  output  motions  were  to  be  displayed  as  spectrum,  subroutines 
using  Fast  Fourier  Transform  (FFT)  methods  were  employed  to  compute  the  spectra. 

Buoy  hull  hydrodynamic  force  functions  in  the  program  were  modified 
according  to  Kim's  data  for  a  half-beam  to  draft  ratio  of  3.2:1.  Terms  are 
included  to  account  for  the  hydrodynamic  mass  and  drag  of  the  three-leg  chain 
bridle  under  the  buoy,  and  the  effective  hydrodynamic  centers  are  modified  to 
account  for  the  bridle.  The  computed  inputs  for  the  torroidal  buoy  are  shown 
in  appendix  D. 

A  three-element,  lumped-mass  cable  model  was  used  in  this  simulation 
since  the  actual  masses  in  the  system  were  concentrated  at  three  places  along 
the  mooring  line  (the  current  meter,  the  sentinel,  and  halfway  down  the  3/4-in. 
chain).  Initial  runs  using  a  five- lump  cable  model  were  compared  with  runs 
using  a  three- lump  cable  model,  and  no  significant  difference  in  buoy  motions 
was  noted.  However,  there  was  an  ordei  jf  magnitude  increase  in  computational 
speed  since  the  two,  high- natural- frequency  cable  lumps  were  summed  into  the 
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current  meter  lump. 

Input  data  for  this  program  are  as  follows: 

1.  Buoy  hull  radius  (ft) 

2.  Height  of  the  center  of  gravity  of  the  buoy  above  the  mooring  line 
connection  point  on  the  buoy  (ft) 

3.  Height  of  the  mooring  line  connection  point  below  the  buoy  hull  (ft) 

4.  Buoy  weight  (lb) 

2  . 

5.  Buoy  structural  and  floodwater  mass  (lb-sec  /ft) 

2 

G.  Yaw  mass  moment  of  inertia  (lb-sec  /ft) 

2  . 

7.  Pitch  mass  moment  of  inertia  (lb-sec  /ft) 

2 

8.  Roll  mass  moment  of  inertia  (lb-sec  /ft) 

9.  Effective  buoy  wind  drag  coefficient  (dimensionless) 

10.  Effective  buoy  wind  lift  coefficient  (dimensionless) 

2 

11.  Buoy  windage  (profile  area)  (ft  ) 

2 

12.  Buoy  plan  area  (ft  ) 

13.  Height  of  the  wind  center  of  pressure  above  the  center  of  gravity  of  the 
buoy  (ft) 

14.  Mean  wind  speed  causing  the  wind  waves  (ft/sec) 

15.  Wind  duration  (hr) 

1G.  Surface  current  y  component  (ft/sec) 

17.  Surface  current  z  component  (ft/sec) 

18.  Wind  y  component  (ft/sec) 

19.  Wind  z  component  (ft/sec) 

20.  Initial  displacements  of  the  buoy  x,  y,  and  z  (ft) 
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21.  Initial  cable  displacements  (including  the  anchor  clump)  x,  y,  and  z 
for  each  of  four  sets  (ft) 

Again,  the  user  should  estimate  the  highest  natural  frequency  in  the 
system  and  adjust  the  integration  step  size  as  required.  Computed  buoy  motions 
are  sampled  every  12th  time  step,  and  1024,  or  2*^,  samples  are  stored. 

Simple  statistical  estimates  (mean,  variance,  and  standard  deviation)  of  the 
output  motions  are  computed.  Each  data  point  is  smoothed  to  reduce  aliasing 
by  averaging  with  the  four  data  points  closest  to  it  in  time.  Spectra  are  com¬ 
puted  using  FFT  subroutines  and  are  smoothed  for  plotting.  Statements  557 
through  595  plot  the  spectral  levels  in  dB.  A  Stromberg  Datagraphics,  Inc. 
integrated  graphics  system  peripheral  to  the  Naval  Underwater  Systems  Center 
UNIVAC  1108  computer  was  used  for  this  procedure. 
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ENCTSF  UNIVAC  1106  FORTRAN  V  COHpTCaYION.  IT *DTaW55TIC*  TC5S*0C (Tf 


00000*  377777777777 


Buoy  System  Dynamics  for  the  Torroidai 
Buoy  Used  in  WHOI  Mooring  No.  238 

This  program  (figure  B-4)  is  basically  the  same  as  the  previous  program 
except  that  the  mooring  line  forces  and  dynamics  are  simulated  with  the  finite- 
difference  method  described  in  chapter  III  of  the  main  text.  Subroutine  "MOOR" 
(figure  B-5)  takes  the  spatial  values  of  the  six  variables  (strain,  two  angles, 
and  the  velocity  components)  describing  the  cable  motions  and  updates  them  for 
the  next  time  step.  The  tension  and  angles  at  the  top  of  the  cable  are  then  used 
to  compute  the  mooring  line  forces  acting  on  the  buoy. 

For  numerical  stability,  the  value  of  the  tensile  wave  characteristic 
should  never  exceed  the  quotient  of  the  nodal  spacing,  H,  and  the  time  step,  K, 
in  the  subroutine: 

Chi  <  H/k, 

The  list  of  input  values  is  the  same  as  that  for  the  previous  program.  In  this 
particular  program,  the  mooring  line  is  composed  of  two  segments,  and  toe 
cable  weights,  masses,  etc.  are  changed  at  a  cable  length  of  4800  ft.  There 
are  20  nodes  spaced  400  ft  apart  to  simulate  dynamics  of  an  8,000-ft  mooring 
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Figure  B-5.  Finite- Difference  Cable  Dynamics  Simulation  Flow  Chart 
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Appendix  C 


DERIVATION  OF  THE  CABLE  CHARACTERISTICS 


In  order  to  rewrite  the  cable  equations  in  their  "normal"  form,  the 
characteristic  roots  of  the  cable  equations  must  be  determined.  From  equation 
(168)  of  the  main  text,  we  see  that  A 14  +BV+C  ~o  % 
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Multiplying  through  by  the  inverse  of  the  A  matrix  yields 

A' A Vt  +A'6t4+A"C  =o  , 

or 

ut  +a"bvs+a'c  **  0 . 

To  find  the  inverse  of  the  A  matrix,  partition  the  A  matrix  into  four  3-by-3 


matrices: 


JX V  /UW  cos (p 
_jxV  -yUlS  SIN  (p 
o  -^(Vcasip-VsMt) 


-I  0  0 

0  -(1+6)  0 


0  0  ~(l+6)C0S(f> 


S 


Let  A  be  partitioned  in  the  same  wav. 


K,  K2' 

K,  K,J  . 

then 


d  -^ilk  kJ  [r  o 

e  0  J  LK3  kiJ  ~  0  1 
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DK,  -><  K3  -  1 

£>K  -yU  K,  =0 

^  0  1=1  0  £  Mor.^rfl  qultir  . 

EK*+0=I  .-.  K2  =  E'j 

Substituting  into  the  first  equation,  we  see  that 

K3  =  jil, 

F„om  the  second  equation,  we  see  Lhal 

=  i,  DK  =  1  pg" 
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The  A  matrix  is  written 
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Summarizing,  the  six  characteristics  are 
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COMPUTED  INPUT  DATA  FOR  THE  SIMULATIONS 

This  appendix  provides  computed  input  data  for  the  following: 

1.  Steady-state  configurations  of  the  8- ft  torroidal  buoy  and  current 
meter  array  at  Block  Island  Sound  station  BRAVO. 

2.  Steady-state  configurations  of  the  8-ft  torroidal  buoy  and  cable  used 
for  WHOI  Mooring  No.  279. 

3.  Dynamics  of  the  3-2— ft  spherical  buoy  and  cable  at  station  DELTA  off 
New  Harbor,  Block  Island. 

4.  Dynamics  of  the  8-ft  torroidal  buoy  and  cable  at  station  BRAVO  in 
Block  Island  Sound. 

5.  Dynamics  of  the  8-ft  torroidal  buoy  and  cable  used  for  WHOI  Mooring 
No.  238. 


Inputs  for  the  Simulation  of  Steady- State 
Configurations  of  the  8-ft  Torroid  and 
Current  Meter  Array  at  BRAVO 

The  components  used  in  the  buoy  system  are  shown  in  figure  21  and  are 
described  in  chapter  IV  of  the  main  text.  The  input  data  are  as  follows: 

1.  Buoy  diameter  BD1  S.O  ft 

2.  Torroid  section  diameter  BD2  2.5  ft 

3.  Buoy  weight  WB  1200.0  lb 

4.  Maximum  bull  draft  IIM 

301 


2.5  ft 
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At 


5.  Free  draft  (when  displacement  equals  buoy 
weight)  HFREE 

6.  Buoy  windage  WAREA 

7.  Wind  drag  coefficient  WCD 

8.  Cable  diameter  DIA: 

Upper  line  -  5/8-in.  polypropylene 
Lower  line  -  3/8-in.  wire  rope 

9.  Cable  weight  in  water  per  unit  length 
WTC: 

Upper  line  - 
Lower  line  - 

10.  Cable  modulus  of  elasticity  EC: 

Upper  line  - 
Lower  line  - 

11.  Mooring  line  length  SM 

12.  Current  speeds  at  three  depths: 

CUR1,  CUR2,  and  CUR3 

13.  Current  directions  at  three  depths: 

DIR  1,  DIR  2,  and  DIR  3 

14.  Water  depth  DEEP 

15.  Wind  speed  v  component 
1G.  Wind  speed  x  component 


0. 7  it 

LL  84  ft2 
0.971 

0.625  in. 
0.375  in. 


-0.02  lb/ft 
0.2  lb/ft 


1.67  x  105  lb/in.  “ 
1.20  x  107  lb/in.2 


235.0  ft 

1.0  to  1.  77  knots 

0-360  deg 
120*.  0  ft 
0. 0  ft/sec 
0.0  ft/sec. 


The  cable  properties  DIA,  WTC,  and  EC  are  changed  at  a  cable  length  of 
;30.0  ft  from  the  values  for  the  5/8-in.  polypropylene  to  the  values  for  3/8- in. 
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wire  rope.  Changes  in  tension  and  angles  across  the  current  meters  were 
found  by  inserting  the  following  procedures  in  the  program: 

1.  Use  and  IF  statement  to  locate  the  current  meters  along  the  mooring 

line. 

2.  Use  the  computed  tension  and  angles  to  compute  the  force  components 
acting  on  the  top  of  the  current  meter. 

3.  Solve  the  statics  equations  for  the  current  meter  by  using  the  current 
meter  in-water  weight  and  computed  drag  force  components  to  find  the  force 
components  acting  on  the  bottom  of  the  current  meter. 

4.  Compute  the  current  meter  tilt  angle. 

5.  Transform  back  to  cable  coordinates  to  find  the  new  cable  tension  and 
angles. 

6.  Increase  the  cable  length  by  3  ft,  i.e. ,  the  length  of  the  current  meter. 

7.  Continue  integration  down  the  cable. 

The  following  data  were  required  to  accomplish  the  above: 

Current  meter  weight  in  water  WCM  67  lb 

Current  meter  drag  coefficient  CDCM  0.59 

2 

Current  meter  frontal  area  ARCM  2.285  ft  . 

Current  strengths  and  directions  as  functions  of  depth  were  computed  as 
follows: 

1.  Average  strengths  and  directions  from  current  meters  1  and  2. 

2.  For  depths  of  0  to  70  ft,  set  the  upper  layer  strength  and  o-  < 


equal  to  these  values. 
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3.  Test  current  meter  3  strength  against  mean  strength  of  meters  1  and 
2.  If  greater,  set  the  lower  layer  current  equal  to  the  mean  strength  in  the 
upper  layer. 

4.  Current  strength  and  direction  in  the  lower  layer  (70  to  120  ft)  set  equal 
to  current  meter  3  strength  and  direction. 

Integration  step  size  (B)  was  set  equal  to  1.0  ft,  and  the  limiting  depth 
error  bandwidth  DDP  was  set  equal  to  1  ft. 

Inputs  for  the  Simulation  of  Steady- State 
Configurations  of  the  8-ft  Torroid  and 
Cable  Used  for  WHOI  Mooring  No.  279 

The  components  used  in  this  mooring  are  shown  in  figure  27  and  are 
described  in  chapter  IV  of  the  main  text.  Buoy  dimensions  are  the  same  but  the 
buoy  weight  is  increased  to  account  for  the  instruments  in  the  buoy,  the  chain 
bridle,  and  instruments  directly  beneath  the  buoy.  The  input  data  are  listed  as 


follows: 

1. 

Buoy  diameter  BD1 

8. 0  ft 

2. 

Torroid  section  diameter  BD2 

2.5  ft 

3. 

Buoy  weight  WB 

2100.0  lb 

4. 

Maximum  hull  draft  HM 

2.5  ft 

5. 

Free  draft  IIFREE 

1.04  ft 

6. 

Buoy  windage  WARE  A 

19.84  ft! 

7. 

Buoy  wind  drag  coefficient 

0.971 

8. 

Cable  diameter  DIA: 

Upper  cable  -  1/4-in.  GAC 

0.  25  in 
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Lower  cable  -  5/8- in.  plaited  nylon  0.625  in. 

9.  Cable  weight  in  sea  water  per  unit  length  WTC: 

Upper  cable  -  0.090  lb/ft 

Lower  cable  -  0.0105  lb/ft 

t 

10.  Cable  modulus  of  elasticity  EC: 

7  2 

Upper  cable  -  1.682  x  10  Ib/in. 

5  2 

Lower  cable  -  3.52  x  10  lb/in.  0  <  T<  1000  lb 

6,  79  x  105  Ib/in.  2  1000  lb  <  T  <  2000  lb 
1.  041  x  106  lb/in. 2  2000  lb  <  T 


11. 

Mooring  line  length  SM 

8000.0  ft 

12. 

Surface  current  CUR 

0-1.46  knots 

13. 

Water  depth  DEEP 

8800.0  ft 

14. 

Wind  speed  y  component 

0.0  ft/sec 

15. 

Wind  speed  z  component 

0.0  ft/sec. 

Cable  properties  were  changed  at  a  length  of  4800  ft  from  the  buoy  and  the 
changes  in  tension  and  angles  were  computed  across  the  instruments  in  the  line 
as  before.  The  input  data  for  the  current  meter  and  tensiometers  are  listed  as 
follows: 

Current  meter: 

Weight  in  water  120.0  lb 

Drag  coefficient  1.4 

2.92  ft2 


Frontal  area 
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Tensiometer: 

Weight  in  water 
Drag  coefficient 
Frontal  area 


50.0  lb 
1.4 

0.875  ft^ 


The  drag  force  acting  on  the  buoy  was  increased  to  account  for  the  chain 

bridle  and  the  instruments  directly  below  the  buoy.  An  effective  area-drag 

2 

coefficient  product  of  22. 18  ft  was  added  to  the  area-drag  coefficient  product 
of  the  buoy  to  account  for  the  instruments  and  chains. 

Inputs  for  the  Simulation  of  Buoy  System 
For  the  Spherical  Buoy 

The  spherical  buoy  dimensions  were  mea.-ured  snd  the  buoy  and  its 
instrumentation  were  weighed.  The  center  of  gravity  and  mass  moments  of 
inertia  were  calculated  from  the  known  weights  and  dimensions.  The  cable  and 
chain  were  weighed  in  air  and  the  mass  and  weight  in  sea  water  of  each  were 
computed.  Mean  wind  wave  heights  and  periods  and  the  frequencies,  amplitudes, 
and  phases  were  computed  on  the  GSA  computer  by  using  program  RWAVE  (see 
subroutine  RWAVE  in  appendix  A)  for  a  given  wind  speed,  a  10-hr  duration,  and 
a  15-mile  fetch.  Input  data  are  listed  as  follows: 


1. 

Buoy  hull  radius  BR 

1.75  ft 

2. 

Center  of  gravity  height  XML 

0.208  ft 

3. 

Mooring  line  connection  height  XML 

0. 20b  ft 

4. 

Buoy  weight  WB 

440.0  lb 

5. 

Buoy  mass  MB 

13.60  lb-secVft 

6. 

Yaw  mass  moment  of  inertia  ALIN 

1 1. 967  lb-see2/ft 
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7. 

Pitch  mass  moment  of  inertia  BTIN 

17.46  lb-sec  /ft 

8. 

Roll  mass  moment  of  inertia  GMIN 

19.  875  Ib-sec2/ft 

9. 

Wind  drag  coefficient  WCD 

0.5 

10. 

Windage  (profile  area)  WAREA 

5.435  ft2 

n. 

Wind  lift  coefficient  WCL 

0.25 

12. 

Plan  area  WAREL 

4. 81  ft2 

13. 

Wind  center  of  pressure  height 

1.068  ft 

14. 

Mean  wave  height  WHTM 

0.5  to  4.05  ft 

15. 

Mean  wave  period  PER 

2. 2  to  6. 0  sec 

16. 

Unstretched  cable  lengths  CLO(I) 

18.  75  ft 

18. 75  ft 

18. 75  ft 

18.75  ft 
12.5  ft 
12.5  ft 

17.  Upper  cable  diameter  DIA1  0.0416  ft 

18.  Upper  cable  weight  in  sea  water  DWC1  0.35  Ib/ft 

2  2 

19.  Upper  cable  mass  per  unit  length  DCSMl  0.0124  lb-sec  /ft 

20.  Lower  cable  diameter  DLA2  0.125  ft 

21.  Lower  cable  weight  in  sea  water  DWC2  7.3  lb/ft 

2  2 

22.  Lower  cable  mass  per  unit  length  DCSM2  0.233  lb-sec  /ft 

23.  Surface  current  y  component  CYS  0.0  ft/sec 

24.  Surface  current  z  component  CZS  0.845  ft/sec 


25.  Wind  y  compor.ent 


-30.0  ft/sec 


•gwe *<j*rfs*Js**>. 
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26.  Wind  z  component 

-8.0  ft/sec 

4.^ 

27.  Initial  buoy  displacements  y(2) 

-0.21  ft 

y(4) 

-41.1  ft 

y(fi) 

62. 35  ft 

28.  Initial  cable  element  displacements  y (I)  15.04  ft 

27.54  ft 
38.  76  ft 

49.46  ft 
55.95  ft 
62. 00  ft 

y(I  +  2)  -29.79  ft 

-20.54  ft 
-13.06  ft 
-6.86  ft 
-3.31  ft 
0.00  ft 

y  (1  +  4)  55.86  ft 

46-06  ft 

34.47  ft 

20.48  ft 
10.43  ft 

0.00  ft. 

The  highest  natural  frequency  of  the  system  was  esti;..'*ted  by  using  th" 
method  described  in  appendix  A,  and  a  time  step  size  (11/  of  Ox  lr-  '  aoc  was 


used.  The  currents  were  uniform  over  the  62— ft  water  depth.  A  30— lb  force 
was  added  to  the  vertical  cable  force  component  acting  on  the  buoy  to  account 
for  the  instiument  package  just  below  the  buoy. 


Again,  buoy  dimensions  and  weights  were  measured  or  taken  from  the 
manufacturers' drawings,  and  the  properties  of  the  buoy  werf  calculated.  Thechain 
bridle  below  the  buoy  was  assumed  to  l>e  rigid,  and  its  mass  and  drag  were 
included  in  the  computations.  In  this  case,  each  cable  element  weight,  mass, 
and  drag  were  listed  directly  in  the  program  since  the  mooring  line  was  com¬ 
posed  of  many  elements  (cable  chains,  current  meter,  sentinel,  etc.).  Input 
data  arc  listed  as  follows: 


1.  Buoy  hull  radius  BTl 

2.  Center  of  gravity  height  XCG 

3.  Mooring  line  connection  height  XML 

4 .  Buoy  weight  WB 

5.  Buoy  mass  MB 

6.  Yaw  mass  moment  of  inertia  ALIN 

7.  Pitch  mass  moment  of  inertia  BTIN 

8.  Roll  mass  moment  of  inertia  GMIN 

9.  Wind  ag  coefficient 

10.  Wind  lift  coefficient 


4.0  ft 
13.064  ft 
11.782  ft 
2100.0  lb 
65.25  lb 

283.44  lb-sec2/ft 
445. 12  lb-sec2/ft 
445. 12  lb-sec2/ft 
0.971 
0.25 

2 

19.84  ft 


11.  Windage  (profile  area) 
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13. 

Wind  center  of  pressure  height 

6. 19  ft 

14. 

Wind  speed  causing  the  waves 

(11  June  1970) 

16. 9  ft/sec 

15. 

Wind  duration 

16. 

Surface  cur  rail  y  component 

(11  June  1970) 

-0. 21  ft/sec 

17. 

Surface  current  z  component 

(11  June  1970) 

-0.597  ft/sec 

18. 

Wind  y  component 

0.0  ft/sec 

IS. 

Wind  z  compcneni 

(11  June  1970) 

-16.9  ft/sec 

20. 

Initial  buoy  displacements  y(2) 

-0.25  ft 

y(4) 

-0.7  ft 

y(6) 

208.6  ft 

21. 

Initial  cable  element  displacements  y(I) 

64  18  ft 

110.45  ft 

120.91  ft 

120.9  ft 

y(I  +  2)  -0.8358  ft 

-0.7416  ft 
-0.6817  ft 
0.0  ft 

y(I  +  4)  200.84  ft 


193.35  ft 
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14 '1.95  ft 

0.0  ft. 

-3 

Jn  this  case,  the  time  step  size  was  taken  to  be  5x10  sec.  Because 
of  the  chain  bridle,  the  hydrodynamic  force  moment  arm  was  computed  to  be 
4. 16  ft  below  the  center  of  gra<  tty.  The  surge  and  sway  area-drag  coefficient 
products  include  the  effects  o.  the  steel  bracing  and  chain  bridle  under  the  buoy. 

Inputs  for  the  Simulation  of  Buoy  System 
Dynamics  for  the  Torroidal  Buoy  and 
Cables  Used  in  WIIOI  Mooring  No,  238 

Buoy  dimensions,  weights,  masses,  etc.  are  the  same  as  those  for  the 
torroidal  buoy  at  station  BRAVO.  The  wind  fetch  length  in  subroutine  It  WAVE 
was  changed  to  100  miles  to  better  simulate  deep-sea  wind  wave  conditions. 
Cable  dynamics  were  simulated  with  subroutine  MOOR,  and  the  six  dependent 
cable  properties  (tension,  two  angles,  and  three  velocity  components)  were 
calculated  at  20  points  along  the  cable.  Total  unstretched  mooring  line  length 
was  8000  ft:  4800  ft  of  1/4-in.  galvanized  steel  aircraft  cable  (polyolefin 
jacketed  to  3/8-in.  diameter)  and  3200  ft  of  5/8-in.  plaited  nylon. 

The  upper  cable  properties  arc  listed  as  follows: 


1. 

Diameter 

0.0312  ft 

2. 

Weight  per  unit  length  in  sea  water 

0.125  lb/ft 

3. 

Mass  per  unit  length 

0.0054  lb-secVft' 

1. 

Characteristic  velocity 

11,300.0  ft/sec 

5. 

Elastic  modulus-cross  section  product 

6. 87  x  10°  lb 
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The  lower  cable  properties  are  listed  as  follows: 


1. 

Diameter 

0.052  ft 

2. 

Weight  per  unit  length  in  sea  water 

0.0105  Ib/ft 

3. 

Mass  per  unit  length 

0.00455  lb-sec2/ft 

4. 

Characteristic  velocity 

2600.0  ft/sec 

5. 

Elnetie  modulus-cross  section  product 

3.06  x  104  lb. 

Cable  properties  were  changed  from  wire  rope  to  nylon  rope  at  4800  ft 
(1  =  14).  The  Webster  current  profile  was  used  to  compute  steady  drag  forces 
on  the  cable.  Numerical  stability  was  maintained  by  using  a  time  step  (d)  of 
0.02  sec.  The  H/K  quotient  is  equal  to  20,000  ft/sec,  which  is  larger  than  the 
tensile  wave  speed  in  the  upper  cable  (11,300  ft/sec).  Initial  strains  and  cable 
angles  were  computed  with  the  steady-state,  buoy  system  configuration  program 
with  the  same  current  profile  and  a  surface  current  of  1.5  knots.  Winds  of  10, 
20,  and  30  knots  were  used,  and  the  computed  steady-state  strains  and  cable 
angles  served  as  input  (initial  conditions)  '  3  dynamics  model. 
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BUOYANT  FORCES  AND  MOMENTS  FOR 
A  TORROIDAL  BUOY 

In  order  to  simulate  the  dynamics  of  the  torroidal  buoy,  a  method  to 
compute  the  buoyant  force  and  righting  moment  for  a  given  buoy  draft  and  tilt 
angle  was  developed.  This  computation  was  included  in  the  program  as  sub¬ 
routine  TORBU,  and  it  updated  the  buoyant  force  and  moment  for  each  integra¬ 
tion  time  step. 

Consider  a  torroid  with  major  radius  R  and  minor  radius  r  ,  partially 
immersed  in  a  fluid  with  a  mean  draft  H  and  a  tilt  angle  0  (figure  E-l  ,  The 
area  of  an  immersed  circular  segment  with  draft  h  is  given  by 
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Thus,  the  area  becomes 
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Define  the  section  draft,  H  ,  as  a  function  of  <p  ,  the  radial  angle  about  the 

s 

torroid  axis.  The  maximum  draft  in  the  direction  of  the  tilt  angle  0  (<P~  ’W/z) 


is 
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and  the  minimum  draft  is 
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For  any  angle  <f>  around  the  torroid, 


Hs  =  H  +  R  SlN<t>SlN9  -+•  (Z.q-  H)(l- sing)  , 
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The  last  term  represents  the  small  change  in  the  mean  draft  H  due  to  the  fact 
that  the  tilt  is  about  the  center  of  gravity  and  not  about  the  waterplane  center. 

By  using  the  draft  for  any  section  as  defined  above,  we  can  compute  the 
immersed  area  of  any  section  from  equation  (E-l)  or  (E-3).  The  immersed 
volume  is  found  by  integrating  around  the  torroid, 
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and  the  buoyant  force 
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where  is  the  weight  density  of  the  fluid.  The  centroid  of  the  submerged 
volume  V  is  computed  in  order  to  find  the  righting  moment: 


(£-7) 


The  righting  arm  is 
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and  the  righting  moment  is 

•M*  ^  ®  %  (E_9) 

The  above  equations  were  programmed  for  subroutine  TORBU,  and  the 
differential  volumes  were  summed  in  2-deg  steps  around  the  torroid.  Buoy 
draft  and  till  angle  are  the  inputs,  and  the  heave  buoyant  force,  tilt  righting 
moment,  cross-coupled  tilt-heave  buoyant  force,  and  cross-coupled  heave-tilt 
righting  moment  are  the  outputs.  Plots  of  the  displacement,  righting  moment, 
and  the  cross-coupled  force  and  moment  for  the  Richardson  torroid  used  at 
station  BRAVO  are  shown  in  figures  (E-2)  through  (E-5). 


